I. Collective Behavior, From Particles to Fields 


I.A Introduction 


The object of the first part of this course was to introduce the principles of statistical 
mechanics which provide a bridge between the fundamental laws of microscopic physics, 
and observed phenomena at macroscopic scales. 

Microscopic Physics is characterized by large numbers of degrees of freedom; for ex- 
ample, the set of positions and momenta {j;,q;}, of particles in a gas, configurations of 
spins {s;}, in a magnet, or occupation numbers {n,;}, in a grand canonical ensemble. The 
evolution of these degrees of freedom is governed by an underlying Hamiltonian H. 

Macroscopic Physics is usually described by a few equilibrium state variables such as 
pressure P, volume V, temperature 7, internal energy FE, entropy S, etc., which obey the 
laws of thermodynamics. 

Statistical Mechanics provides a probabilistic connection between the two realms. For 
example, in a canonical ensemble of temperature T’, each micro-state pu, of the system 
occurs with a probability p(w) = exp (—BH(u))/Z, where 3 = (kpT)~'. To insure 
that the total probability is normalized to unity, the partition function Z(T) must equal 
+ 1 &XP (—BH(n)). Thermodynamic information about the macroscopic state of the sys- 
tem is then extracted from the free energy F = —kpT |nZ. 

The above program can in fact be fully carried out only for a limited number of sim- 
ple systems; mostly describing non-interacting collections of particles where the partition 
function can be calculated exactly. Some effects of interactions can be included by per- 
turbative treatments around such exact solutions. However, even for the relatively simple 
case of an imperfect gas, the perturbative approach breaks down close to the condensa- 
tion point. One the other hand, it is precisely the multitude of new phases and prperties 
resulting from interactions that renders macroscopic physics interesting. In particular, we 
would like to address the following questions: 

(1) In the thermodynamic limit (N — oo), strong interactions lead to new phases of 
matter. We studied the simplest example of the liquid gas transition in some detail, 
but there are in fact many other interesting phases such as solids, liquid—crystals, 
magnets, superconductors, etc. How can we describe the emergence of such distinct 
macroscopic behavior from the interactions of the underlying particles? What are 


the thermodynamic variables that describe the macroscopic state of these phases; and 


it 


what are their identifying signatures in measurements of bulk response functions (heat 

capacity, susceptibility, etc.)? 

(2) What are the characteristic low energy excitations of the system? As in the case of 
phonons in solids or in superfluid helium, low energy excitations are typically collective 
modes, which involve the coordinated motions of many microscopic degrees of freedom 
(particles). These modes are easily excited by thermal fluctuations, and probed by 
scattering experiments. 

The underlying microscopic Hamiltonian for the interactions of particles is usually 
quite complicated, making an ab initio particulate approach to the problem intractable. 
However, there are many common features in the macroscopic behavior of many such sys- 
tems that can still be fruitfully studied by the methods of statistical mechanics. Although, 
the interactions between particles are very different at the microscopic scale, one may hope 
that averaging over sufficiently many particles leads to a simpler description. (In the same 
sense that the central limit theorem ensures that the sum over many random variables 
has a simple Gaussian distribution.) This expectation is indeed justified in many cases 
where the collective behavior of the interacting system becomes more simple at long wave- 
lengths and long times. (This is sometimes called hydrodynamic limit by analogy to the 
Navier-Stokes equations for a fluid of particles.) The averaged variables appropriate to 
these length and time scales are no longer the discrete set of particle degrees of freedom, 
but slowly varying continuous fields. For example, the velocity field that appears in the 
Navier-Stokes equations is quite distinct from the velocities of the individual particles in 
the fluid. Hence the productive method for the study of collective behavior in interacting 
systems is the Statistical Mechanics of Fields. Thus the program of this course is, 

e Goal: To learn to describe and classify states of matter, their collective properties, and 

the mechanisms for transforming from one phase to another. 

e Tools: Methods of classical field theories; use of symmetries, treatment of nonlinearities 

by perturbation theory, and the renormalization group (RG) method. 

e Scope: To provide sufficient familiarity with the material so that you can follow the 

current literature on such subjects as phase transitions, growth phenomena, polymers, 


superconductors, etc. 


I.B_ Phonons and Elasticity 


The theory of elasticity represents one of the simplest examples of a field theory. We 


shall demonstrate how certain properties of an elastic medium can be obtained, either by 
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the complicated method of starting from first principles, or by the much simpler means of 
appealing to symmetries of the problem. As such, it represents a prototype of how much 
can be learned from a phenomenological approach. The actual example has little to do 
with the topics that will be covered in this course, but it fully illustrates the methodology 
that will be employed. The problem of computing the low temperature heat capacity of a 
solid can be studied by both ab initio and phenomenological methods. 

(i) Ab initio (particulate) approach: Calculating the heat capacity of a solid material 
from first principles is rather complicated. We breifly sketch some of the steps: 

e The ab initio starting point is the Schrodinger equation for electrons and ions which can 
only be treated approximately, say by a density functional formalism. Instead, we start 
with a many body potential energy for the ionic coordinates V(qi, d2,---, dn), which may 
itself be the outcome of such a quantum mechanical treatment. 

e Ideal lattice positions at zero temperature are obtained by minimizing V, typically 
forming a lattice ¢*(€,m,n) = [¢a + mb + nc] = ¢%, where 7 = {€,m,n} is a triplet of 
integers, and G, 6, and é@ are unit vectors. 

e Small fluctuations about the ideal positions (due to finite temperature or quantum 
effects) are included by setting g7 = ¢%+ u(r’). The cost of deformations in the potential 


energy is given by 


Uol®) ug(*’) + O(u?), (1.1) 


where the indices a and (3 denote spatial components. (Note that the first derivative of 
Y vanishes at the equilibrium position.) The full Hamiltonian for small deformations is 
obtained by adding the kinetic energy >) >, Dal? )?/2m to eq.(I.1), where pa(7) is the 
momentum conjugate to uq(7). 

e The next step is to find the normal modes of vibration (phonons) by diagonalizing the 
matrix of derivatives. Since the ground state configuration is a regular lattice, the elements 
of this matrix must satisfy various translation and rotation symmetries. For example, they 
can only depend on the difference between the position vectors of ions 7 and 7°”, i.e. 


0?V 
——_—_ = K,4(F — 7’). 1.2 
are al ) (1.2) 


This translational symmetry allows us to at least partially diagonalize the Hamiltonian by 


using the Fourier modes, 
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tal®) = Fottelk ). (1.3) 
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(Only wavevectors k inside the first Brillouin zone contribute to the sum.) The Hamilto- 
nian then reads 
Lise | poli) 2 ine 
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While the precise form of the Fourier transformed matrix Kap(k ) is determined by the 
microscopic interactions, it has to respect the underlying symmetries of the crystallographic 
point group. Let us assume that diagonalizing this 3 x 3 matrix yileds eigenvalues {k.(k) } 
The quadratic part of the Hamiltonian is now decomposed into a set of independent (non— 
interacting) harmonic oscillators. 


e The final step is to quantize each oscillator, leading to 
4 | 
Hav +S twelE) (mal) +5). (15) 


where wa(k) = \/ka(k)/m, and {ng(k)} are the set of occupation numbers. The average 


energy at a temperature T is given by 
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E(P) =V* + So hwa(k) ((ra()) + 5) ; (1.6) 
kya 


where we know from elementary statistical mechanics that the average occupation num- 


hwa 
kpT 


bers are given by (na(k)) = 1/(exp(##s;) — 1). Clearly E(T), and other macroscopic 
functions, have a complex behavior, dependent upon microscopic details through {k,(k) }. 
Are there any features of these functions (e.g. the functional dependence as T — 0) that 
are independent of microscopic features? The answer is positive, and illustrated with a 


one dimensional example. 


I.1. Displacements of a one-dimensional chain, and coarse-graining. 


Consider a chain of particles, constrained to move in one dimension. A most general 


quadratic potential energy for deformations {u,,}, around an average separation of a, is 


_ K K 
Y=YV fe in Un)? + > Di lins2 = Uy)? oe (1.7) 
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The decomposition to normal modes is obtained via 


n/a : ; 
Un =i aK —iknay(f), where u(k) = Ca : (1.8) 


—t/a 20 


(Note the difference in normalizations from eq. (1.3).) The resulting potential energy is, 


esi la dk ak es Laat 
as yt @_§41\(e*® %_—] i(k+k’' )na : doles 1. 
Pao ee u(Bu(kl) +. (19) 
Using the identity 5 ve rene = 0(k + k’)2n/a, and noting that u(—k) = u*(k), we 


obtain 


n/a 
V=vVrt+ mS a [K1(2 — 2coska) + Ko(2—2cos2ka) + ---] |u(k)|?. (1.10) 
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The frequency of normal modes, given by w(k) = \/[2K1(1 — cos ka) + ---]|/m, is depicted 
below. In the limit k — 0, w(k) — v|k|, where the ‘sound velocity’ v equals a,/K /m (see 


below). 


The internal energy of these excitations, for a chain of N particles, is 


n/a 
BD=V+ Na | a a) 


_ n/a 2m exp (hw(k)/kpT) — 1)" (1.11) 
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1.2. The dispersion of phonon modes for a linear chain. 


As T — 0, only modes with hw(k) < kgpT are excited. Hence only the k — 0 part of the 
excitation spectrum is important and E(T) simplifies to 
°° dk ho|k| n 


E(D) =V" + Na =V*+Na 


SS kpT)?. 1.12 
_oo 27 exp(hv|k|/kpT) — 1 hiv | BT) oe 


Note: 
(1) The full spectrum of excitation energies is 
ao = K,(1—coska) + Ko(1 — cos2ka) +--- => +h as k— 0, (L113) 
Further neighbor interactions change the speed of sound, but not the form of the 
dispersion relation as k — 0. 

(2) The heat capacity C(T’) = dE/dT is proportional to T. This dependence is a universal 
property, i.e. not material specific, and independent of the choice of the interatomic 
interactions. 

(3) The T? dependence of energy comes from excitations with k — 0 (or A — oo), ice. 
from collective modes involving many particles. These are precisely the modes for 
which statistical considerations may be meaningful. 

(ii) Phenomenological (field) approach: We now outline a mesoscopic approach to 

the same problem, and show how it provides additional insights and is easily generalized 

to higher dimensions. Typical excitations at low temperatures have wavelengths A > 

A(T) & (hv/kT) >> a, where a is a lattice spacing. We can eliminate the unimportant 


short wavelength modes by an averaging process known as coarse graining. The idea is 


6 


consider a point x, and an interval a < dx < X(T) around it (Fig. 1.1). In dz all the 
displacements wu are approximately the same; and we can define an average deformation 
field u(x). The kinetic energy is now related to the density p = m/a via p f dxti(x)?/2. 
Note that u(x) is treated as a continuous function, but it certainly does not have any 
variations over distances smaller than the lattice spacing a. 

What is the most general potential energy functional V[u], for the chain? A priori, 
we don’t know much about the form of V|u], but we can construct it using the following 


general principles: 


Locality: In most situations, the interactions between particles are short range, al- 
lowing us to define a potential energy density ® at each point x, with Viu] = 
J dx®(u(x), du/dx,-- -). Naturally, by including all derivatives we can also describe long- 
range interactions. In this context, the term locality implies that the higher derivative 


terms are less significant. 


Translational symmetry: A uniform translation of the chain does not change its internal 
energy, and hence the energy density must satisfy the constraint ®[u(x2) + c] = ®[u(x)]. 
This implies that ® cannot depend directly on u(x), but only on its derivatives du/dz, 
d?u/dx?, 

Stability: Since the fluctuations are around the equilibrium solution, there can be no 
linear terms in u or its derivatives. (Stability further requires that the quadratic part of 
Y|u] must be positive definite. ) 


The most general potential consistent with these constraints is 


K (du\? L (@u\? Ou 
= ay (ene = pa M tee 1.14 
vil ic E (=) 2 (St) ++ pee @) % ee 
which after Fourier transformation gives 
dk 
Yul -/> [> ke+ ai | |u(k)|? 
2n | 2 
(1.15) 
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As k — 0, higher order gradient terms (such as L) become unimportant. Also, for small 
deformations we may neglect terms beyond second order in u (such as M). Adding the 


kinetic energy, we get a simple one dimensional field theory, with a Hamiltonian 


n= 8 fee[() +H) 


7 


This is a one-dimensional elastic (string) theory with material dependent constants p 
and v = /K/p. While the phenomenological approach cannot tell us the value of these 
parameters, it does show that the low energy excitations satisfy the dispersion relation 
Salk: 

We can now generalize the elastic theory of the string to arbitrary dimensions d: The 
discrete particle deformations wu, are coarse grained into a continuous deformation field 
u(Z). For an isotropic material the potential energy V|u], must be invariant under both 
rotations and translations ua(Z) + Ragug(@) + ca, where Rgg is a rotation matrix. A 
useful local quantity is the symmetric strain field, 


1 /0ua . dug 
Uap (#) = 3 te oF — ; (1.16) 


in terms of which the most general quadratic deformation Hamiltonian is 


Otig OUe , 2 pI 
H= S- jez c es Oe + ae ae + F aatla9 ‘ (1.17) 
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The elastic constants 44 and A are known as Lame coefficients. Summing over the repeated 
indices ensures that the result is rotationally invariant. This rotational invariance is more 


> 


transparent in the Fourier basis, u(k) = [ dtzet* Fq(z), since the Hamiltonian 


alk, pasos U ~ pty 2 

= | —— |-lu(k)|? + =k? |a(k)|? + —— (kak 1.18 
=f S [Srieye + Leia? +S Rate’) (1.18) 
manifestly includes only rotationally invariant quantities k - k, u- wu, and k- a. We can 
further decompose the Hamiltonian into two types of sound modes: longitudinal modes 
where k || ad, with vg = ,/(2u+ A)/p, and transverse modes with k L @, where v; = L/p. 


The internal energy is then given by 


d’k huck (d— 1)hurk 
ee / (2Q7)4¢ en —1 - exp(hu:k/kpT) — 1 Avert \ka 
1.19 


The specific heat vanishes as C x T%, as T — 0. 

Note: 

(1) All material dependent parameters end up in the coefficient A, while the scaling with 
T is universal. 

(2) The universal exponent originates from the (hydrodynamic) modes with k > 0. The 


high frequency (short wavelength) modes come into play only at high temperatures. 
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(3) The scaling exponent depends on dimensionality and the range of interactions. (Long— 
range Coulomb interactions lead to a different result.) 

(4) Experimental observation of a power law alerts us to the physics. For example, in 
superfluid helium, the observation of C « T° (as opposed to C « T 3/2 expected for 
an ideal Bose gas), immediately implies phonon-—like excitations as noted by Landau. 

x There are many other well known examples demonstrating the universality and impor- 

tance of power laws. For example, consider a cloud of tracers moving in some unspecified 

medium. The scaling of some characteristic dimension x with time ¢t can alert us to the 
possible dynamics that governs the motion of the particles. Three simple possibilities are: 

(1) Diffusion, in which case x x V Dt. 

(2) Dissipate transport, where x x vt. 

(3) Free forced motion, where x x gt?/2, as in a gravitational field. 

The Navier-Stokes equation for fluid flow is yet another example. We can use these 

examples to construct general guidelines for setting up and analysing phenomenological 

field theories. Some of steps in the procedure are: 

(1) Input for construction of the coarse grained Hamiltonian comes from symmetry, range 
of interactions, and dimensionality. 

(2) Unlike the above example, in general nonlinearities cannot be ignored in the resulting 
effective field theory. We shall learn how to treat such nonlinearities by the methods 
of perturbation theory and the renormalization group. 

(3) Output of the analysis is expressed in terms of universal exponents, and other func- 


tional dependances that can be directly compared with experiments. 


I.C Phase Transitions 


The most spectacular consequence of interactions among particles is the appearance 
of new phases of matter whose collective behavior bears little resemblance to that of a few 
particles. How do the particles then transform from one macroscopic state to a completely 
different one. From a formal perspective, all macroscopic properties can be deduced from 
the free energy or the partition function. Since phase transitions typically involve dramatic 
changes in various response functions they must correspond to singularities in the free 
energy. The canonical partition function for a finite collection of particles is always an 
analytical function. Hence phase transitions, and their associated non—analyticities, are 
only obtained for infinitely many particles, i.e. in the thermodynamic limit, N — oo. The 
study of phase transitions is thus related to finding the origin of various singularities in 
the free energy and characterizing them. 

The classical example of a phase transition is the condensation of a gas into a liquid. 
Some important features of the liquid—gas condensation transition are: 

(1) In the temperature/pressure plane, (7, P), the phase transition occurs along a line 
that terminates at a critical point (T., Pe). 

(2) In the volume/pressure plane, (P,v = V/N), the transition appears as a coexistence 
interval, corresponding to a mixture of gas and liquids of densities p, = 1/v,, and 
pi = 1/v, at temperatures T < Ty. 

(3) Due to the termination of the coexistence line, it is possible to go from the gas phase to 
the liquid phase continuously (without a phase transition) by going around the critical 
point. Thus there are no fundamental differences between liquid and gas phases. 

From a mathematical perspective, the free energy of the system is an analytical func- 

tion in the (P,T) plane, except for some form of branch cut along the phase boundary. 

Observations in the vicinity of the critical point further indicate that: 

(4) The difference between the densities of coexisting liquid and gas phases vanishes on 
approaching T%, i.e. Pliquid > Pgas, a8 T — Tr. 

(5) The pressure versus volume isotherms become progressively more flat on approaching 
Tc from the high temperature side. This implies that the isothermal compressibility, 
Kr = — OV/OP|,, /V, diverges as T > T;. 

(6) The fluid appears “milky” close to criticality. This phenomenon, known as critical 
opalescence, suggests collective fluctuations in the gas at long enough wavelengths to 
scatter visible light. These fluctuations must necessarily involve many particles, and 


a coarse graining procedure may thus be appropriate to their description. 
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A related, but possibly less familiar, phase transition occurs between paramagnetic 
and ferromagnetic phases of certain substances such as iron or nickel. These materials 
become spontaneously magnetized below a Curie temperature 7,. There is a discontinuity 
in magnetization of the substance as the magnetic field h, goes through zero for T < Ty. 
The phase diagram in the (h,7) plane, and the magnetization isotherms M(h), have 
much in common with their counterparts in the condensation problem. In both cases a 
line of discontinuous transitions terminates at a critical point, and the isotherms exhibit 
singular behavior in the vicinity of this point. The phase diagram of the magnet is simpler 
in appearance, because the symmetry h ++ —h ensures that the critical point occurs at 
ha Me= 0: 


I.D Critical Behavior 


The singular behavior in the vicinity of a critical point is characterized by a set of 
critical exponents. These exponents describe the non—analyticity of various thermodynamic 
functions. The most commonly encountered exponents are listed below: 

e The Order Parameter: By definition, there is more than one equilibrium phase on 
a coexistence line. The order parameter is a thermodynamic function that is different 
in each phase, and hence can be used to distinguish between them. For a magnet, the 
magnetization F 

mL) = V lim M(h,T), 


serves as the order parameter. In zero field, m vanishes for a paramagnet and is non—zero 


in a ferromagnetic, i.e. 


0 for ST 


Ee foe Peet (1.20) 


m(T,h=0) x | 


where t = (JT. — T)/T. is the reduced temperature. The singular behavior of the order 
parameter along the coexistence line is therefore indicated by the exponent (3. The singular 
behavior of m along the critical isotherm is indicated by another exponent 6, defined 
through 

mT =Te,h) x hi, (1.21) 


The two phases along the liquid—gas coexistence line are differentiated by their density, 
and the density difference p — p., where p, is the critical density, serves the role of the 


order parameter. 


11 


e Response Functions: The critical system is quite sensitive to external perturbations, 
as typified by the infinite compressibility at the liquid—gas critical point. The divergence 
in the response of the order parameter to a field conjugate to it is indicated by an exponent 


y. For example, in a magnet, 


xe A= 0) ole 5 (1.22) 


where in principle two exponents y; and y_ are necessary to describe the divergences on 
the two sides of the phase transition. Actually in almost all cases, the same singularity 
governs both sides and y, = y_ = y. The heat capacity is the thermal response function, 


and its singularities at zero field are described by the exponent aq, i.e. 


Cu(T,h =0) ox |t\-*. (1.23) 


e Long-range Correlations: Since the response functions are related to equilibrium 
fluctuations, their divergence in fact implies that fluctuations are correlated over long dis- 
tances. We shall prove this statement by considering the magnetic susceptibility. Starting 
from the (Gibbs) partition function in a field h, Z(h) = tr{exp|—BHo + GhM]}, the mag- 
netization can be computed as (M) = 01InZ/O(Gh) = tr{M exp|—BHpo + BhM]}/Z. The 
susceptibility is then related to the variance of magnetization by 


y= OM =3 ize [M_? exp (—GHo + BhM)| > = tr [M exp (—GHo + samy} } 
Z Z (1.24) 


=" ((M?) = (M)”) 


The overall magnetization is obtained by adding contributions from different parts of 


the system, i.e. 


M= fe rm(r). (1.25) 


(For the time being we treat the magnetization as a scalar quantity.) Substituting the 


above into eq.(I.24) gives 
kpT x = [ bra" (m@ymir") — (m(¥)) (m(F"))) . (1.26) 


Translational symmetry of a homogeneous system implies that (m(7’)) = m is a constant, 
while (m(7’)m(r’)) = G(r —7') depends only on the separation. We can express the result 


in terms of the connected correlation function, defined as 
(m(F)m(F")), = ((m(#) — (m(F))) (MF) — (m(F")))) = GF —F)— m7, (1.27) 
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Integrating over the center of mass coordinates in eq.((I.26)) results in a factor of volume 


V, and the susceptibility is given by 
X= av f avr (m(F’)m(0)),. (1.28) 


The connected correlation function is a measure of how the local fluctuations in one 
part of the system effect those of another part. Typically such influences occur over a 
characteristic distance €, called the correlation length. (It can be shown rigorously that 
this function must decay to zero at large separations; in many cases G,(7") = (m(7)m(0))., 
decays as exp(—|7’|/&) at separations |7’| > €.) Let g denote a typical value of the correla- 
tion function for |7'| < €. It then follows from eq.(I.28) that kgTy/V < g&?; and x — ox, 
necessarily implies € — oo. This divergence of the correlation length also explains the ob- 
servation of critical opalescence. The correlation function can be measured by scattering 


probes and its divergence 


(Th = 0) x fel, (1.29) 


is controlled by exponents vp =v_ =v. 
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Il. THE LANDAU-GINZBURG APPROACH 


Il.A Introduction 


We noted in the previous section that the singular behavior of thermodynamic func- 
tions at a critical point (the termination of a coexistence line) can be characterized by 
a set of critical exponents {a, 3,y,---}. Experimental observations indicate that these 
exponents are quite universal, i.e. independent of the material under investigation, and 
to some extent, of the nature of the phase transition. For example, the vanishing of the 
coexistence boundary in the condensation of C'O2 has the same singular behavior as that 
of the phase separation of protein solutions into dilute and dense components. This univer- 
sality of behavior needs to be explained. We also noted that the divergence of the response 
functions, as well as direct observations of fluctuations via scattering studies, indicate that 
the correlations become long ranged in the vicinity of the critical point. Such correlated 
fluctuations involve many particles (€ >> a, where a is a typical inter-particle spacing), and 
a coarse graining approach, in the spirit of the theory of elasticity, may be appropriate to 
their description. Here we shall construct such a statistical field theory. 

We shall frame the discussion in the language of a magnetic system, whose symmetries 
are more transparent, although the results are of more general applicability. Consider a 
metal, say iron, close to its Curie point. The microscopic origin of magnetism is quantum 
mechanical, involving such elements as itinerant electrons, their spin, and the exclusion 
principle. Clearly a microscopic approach is rather complicated, and material dependent. 
Such a theory is necessary to find out which elements are likely to produce ferromagnetism. 
However, given that there is such behavior, the microscopic theory is not necessarily useful 
to describe its disappearance as a result of thermal fluctuations. This is because the 
(quantum) statistical mechanics of the collection of interacting electrons is excessively 
complicated. The important degrees of freedom close to the Curie point, whose statistical 
mechanics is responsible for the phase transition, are long wavelength collective excitations 
of spins (much like the long wavelength phonons that dominate the heat capacity at low 
temperatures). We can thus coarse grain the magnet to a scale much larger than the 
lattice spacing, and define a magnetization field m™(x), which represents the average of the 
elemental spins in the vicinity of a point x. It is important to emphasize that while x is 


treated as a continuous variable, the functions m do not exhibit any variations at distances 
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of the order of the lattice spacing, i.e. their Fourier transforms involve only wave-numbers 
less than some upper cutoff A ~ 1/a. 

In describing other types of phase transitions, the role of (x) is played by the appro- 
priate order parameter density. It is then useful to examine a generalized magnetization 


for n-component spins, existing in d-dimensional space, i.e. 
= d he é 
X = (21, %2,...,ta) ER (space) , m= (m1,mMg,...,Mn) € R” (spin). 


Some specific problems covered in this framework are: 
n = 1 describes liquid—gas transitions, binary mixtures, as well as uniaxial magnets; 
n = 2 applies to superfluidity, superconductivity, and planar magnets; 
n = 3 corresponds to classical magnets. 
While most physical situations occur in three-dimensional space (d = 3), there are also 
important phenomena on surfaces (d = 2), and in wires (d = 1). Relativistic field theory 
is described by a similar structure, but in d= 4. 

As in the case of a deformed solid, we construct a local effective Hamiltonian GH([m] = 
{ d*x®[m(x)], on the basis of appropriate symmetries. We shall assume that the material 
is uniform and isotropic, so that all locations and directions in x space are equivalent. In 
the absence of an external magnetic field, all directions for magnetization are equivalent, 
and hence H[R,m(x)] = Him(x)], where R, is a rotation in the n-dimensional space. 
Some of the terms consistent with these symmetries that can appear in the expansion of 


®[m7(x)] are 


m?(x) = m7i(x) - (x) = S| mj(x)mj(x) , m(x)= (m?(x))° , m(x) , ; 
n d 
(Vm) =S > S damdam, , (Vm)? , m(vVm)? , 


Including a small magnetic field h, that breaks the rotational symmetry, the lowest order 


terms in the expansion of ® lead to, 
a, |t, 2 4 K 2 Lb. 
BCS: [aes 5m (x) +um (x) + = (Vm) +---—h-m(x)], (II.1) 


which is known as the Landau-Ginzburg Hamiltonian. (The magnetic field also generates 
terms proportional to m?77- h, which are of higher order, and actually less important than 


the um* term.) 
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Eq.(II.1) is constructed on the basis of symmetry alone, and depends on a set of phe- 
nomenological parameters {t,u, K,---}. These parameters are non-universal functions of 
microscopic interactions, as well as external parameters such as temperature and pressure. 
It is essential to fully appreciate the latter point, which is usually the source of much con- 
fusion. The probability for a particular configuration of the field is given by the Boltzmann 
weight exp{—GH|m7i(x)|}. This does not imply that all terms in the exponent are propor- 
tional to (kgT)~'. Such dependence holds only for the true microscopic Hamiltonian. 
The Landau—Ginzburg Hamiltonian is more correctly an effective free energy obtained by 
integrating over (coarse graining) the microscopic degrees of freedom, while constraining 
their average to ™m(x). It is precisely because of the difficulty of carrying out such a first 
principles program that we postulate the form of the resulting effective free energy on the 
basis of symmetries alone. The price paid is that the phenomenological parameters have 
an unknown functional dependence on the original microscopic parameters, as well as such 
on external constraints as temperature (since we have to account for the entropy of the 


short distance fluctuations lost in the coarse graining process). 


II.B Saddle Point Approximation, and Mean—Field Theory 


The original problem has been simplified considerably by focusing only on the coarse 
grained magnetization field described by the Landau-Ginzburg Hamiltonian in eq.(II.1). 
Various thermodynamic functions (and their singular behavior) can now be obtained from 


the associated partition function 
ae i: Diti(x) exp{—BH[im1(x)]}. (11.2) 


Since the degrees of freedom appearing in the Hamiltonian are functions of x, the symbol 
J Dii(x) refers to a functional integral. In actuality, the functional integral should be 
regarded as a limit of discrete integrals. After discretizing the coordinate x into a lattice 


of N points i, at a distance a from each other, 


N 
Seta | 


(There are some mathematical concerns regarding the existence of functional integrals. 
The problems are associated with having too many degrees of freedom at short distances, 


allowing rather badly behaved functions. These issues need not concern us since we know 
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that the underlying problem has a well defined lattice spacing that restricts the short 
distance behavior.) 

It is still difficult to calculate the Landau—Ginzburg partition function. As a first step, 
we perform a saddle point approximation in which the integral in eq.(II.2) is replaced by 
the maximum value of the integrand. The natural tendency of interactions in a magnet 
is to keep the magnetizations vectors parallel, and hence we expect the parameter K in 
eq.(II.1) to be positive. The configuration of m™ that maximizes the integrand is then 
uniform, and 


BF =-InZ=V min{V(m)}m. (11.3) 


The uniform magnetization occurs for ™(x) = mh, which minimizes 
m +umt +---— Adm. (II.4) 


In the vicinity of the critical point, m is a small quantity, and it is justified to keep 
only the lowest powers in the expansion of V(m). (We can later check self consistently that 
the terms left out are indeed small corrections.) The behavior of U(m) depends strongly 
on the sign of the parameter t. 

(1) For t > 0, we can ignore the quartic term, and the minimum occurs for ™ ~* h/t. The 
vanishing of magnetization as h—-0 signals paramagnetic behavior. The susceptibility 

x = 1/t, diverges as t > 0. 

(2) For t < 0, a quartic term with a positive value of u is required to insure stability 

(i.e. a finite magnetization). The function Y(m) now has degenerate minima, at 

a non-zero value of m. There is thus a spontaneous magnetization, even at h = 

0 indicating ferromagnetic behavior. The direction of the m is determined by the 

systems preparation, and can be realigned by an external field h. 

Thus a saddle point evaluation of the Landau—Ginzburg partition function results in 
paramagnetic behavior for t > 0, and ferromagnetic behavior for t < 0. Hence we can map 


the phase diagram of the Landau—Ginzburg Hamiltonian to that of a magnet by setting 


(fo \=or HP) tO aT": 


(II.5) 
u(T,-++) =u+u(T — To) + O(F — Te)’, 


where a and u are unknown positive constants, dependent upon material properties. The 
basic idea is that the phenomenological parameters are functions of temperature that can 


be expanded in a Taylor series in T — T,.. The minimal conditions needed to reproduce the 
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experimental phase diagram are contained in eqs.(II.5). It is of course possible that some 
other terms in the expansion, e.g. a or u are also zero. However, these are non-generic 
situations which can presumably be removed by changing some other system parameter. 


We now examine the singular behaviors predicted by eqs.(II.3) and (II.4). 


e Magnetization: In zero field, from OV/Om = tm + 4um3? = m(t + 4um?) = 0, and we 


obtain 


(II.6) 


0 for t > 0, 
= fort < 0. 


We thus find a universal exponent 3 = 1/2, while the amplitude is material dependent. 


e Heat Capacity: The free-energy is given by (h = 0) 


(are mec fort, 
y= min U(m) = { ae eG. (II.7) 
Since t = a(T — T.) +--+; 0/0T x O/0t, and 
7 f O? (BF 0 for t > 0, 
C= -Taa X —a5 (=) = { + fort <0. IES) 


We observe a discontinuity, rather than a divergence, in the heat capacity. If we insist 

upon describing the singularity by a power law, we have to choose the exponent a = 0. 

e Susceptibility: In the presence of h, we expect 77 = m(h)h, and from 0V/dm = 0, we 

obtain tm + 4um? = h. Hence 
_1_ Oh 


So ep n? = { 
Xe dm |, —o oe um 


t for t > 0, 


II. 
—2t for t < 0. (11.9) 


Thus the singularity in susceptibility is describable by y+ ~ A+|t|~7*, with 7, = 7_ = 


1. Although the amplitudes A+ are material dependent, their ratio is predicted to be 


universal, with A;/A_ = 2. (As we shall see later, what we have calculated so far is a 
longitudinal susceptibility. There is also a transverse susceptibility that is always infinite 
below Ti.) 
e Equation of State: On the critical isotherm t = 0, the magnetization behaves as 
m = (h/4u)/3, ice. 

m(t=0,h)~hi/?, with 6 =3. (II.10) 
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II.C Spontaneous Symmetry Breaking and Goldstone Modes 


For zero field, h= 0, although the microscopic Hamiltonian has full rotational sym- 
metry, the low-temperature phase does not. As a specific direction in n-space is selected 
for the net magnetization M , there is a spontaneously broken symmetry, and a correspond- 
ing long-range order is established in the system. The original symmetry is still present 
globally, in the sense that if all local magnetizations (x), are rotated together (i.e. if 
m(x) > RmM(x)), there is no change in energy. Such a rotation transforms one ordered 
state into an equivalent one. If a uniform rotation costs no energy, by continuity we expect 
a rotation that is slowly varying in space (e.g. (x) + R(x)m(x), where R(x) only has 
long wavelength variations) to cost very little energy. Such low energy excitations are called 
Goldstone modes. They are present in any system with a broken continuous symmetry. 
There are no Goldstone modes when a discrete symmetry is broken, since it is impossible 
to produce slowly varying rotations from one state to an equivalent one. Phonons are an 
example of Goldstone modes, corresponding to the breaking of translation and rotation 
symmetries by a crystal structure. 

Let us explore the origin and consequences of Goldstone modes in the context of 
superfluidity. In analogy to Bose condensation, the superfluid phase has a macroscopic 


occupation of a single quantum ground state. The order parameter, 
h(x) = de + ive = |b(x) |e”, (II.11) 


is the ground state component (overlap) of the actual wavefunction in the vicinity of x. 
The phase of the wavefunction is not an observable quantity and should not appear in any 
physically measurable probability. For example, the effective coarse grained Hamiltonian 


can be obtained as an expansion, 
a, |A 2, !) 12 4 
BH = | d°x rae +5 +ulpl +--+]. (II.12) 


Clearly, eq.(II.12) is equivalent to the Landau-Ginzburg Hamiltonian with n = 2 (% = 
(wx, Wy)). The superfluid transition is signaled by the onset of a finite value of w for t < 0. 
Minimizing the Hamiltonian fixes the magnitude of ~, but not its phase 8. Now consider 
a state with a slowly varying phase, i.e. with q(x) = we’? Inserting this form in the 


Hamiltonian yields an energy 
K d 2 
BH = BHo + 5 a x( VO), (11.13) 
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where K = Ky?. Taking advantage of translational symmetry, Eq.(IJ.13) can be decom- 
posed into independent modes (in a region of volume V) by setting 0(x) = )/, ety /V/V, 
as 


K 
BH = BHo+ > > 4° |9(a)/?. (11.14) 
q 


Clearly the long wavelength Goldstone modes cost little energy and are easily excited by 
thermal fluctuations. 
Assuming that the amplitude of the order parameter is indeed uniform, the probability 


of a particular configuration is given by, 
ee ore 2 
P[O(x)] « exp =r d°x(V0)*| . (11.15) 
Alternatively, in terms of the Fourier components, 
K 2 2 
P[A(a)] x exp |—> D @lA(a)I?| x [] p(Ga). (11.16) 
q q 


Each mode Og, is an independent random variable with a Gaussian distribution of zero 


mean, and with 1 


Og —al 
(Oq9q') = oe (11.17) 


' Note that the Fourier transform of a real field (x), is complex 0g = 0g, + i0q,3- 
However, the number of fields is not doubled, due to the constraint of 6g = 04 = 9q,." — 


i0q,5. A Gaussian translational invariant weight has the generic form 


P [{Oq}] « Heo A) 646-4] = [Tow |- 2A (02 5 + 02 S) 


q>0 


While the first product is over all q, the second is restricted to half of the space. There 


are clearly no cross correlations for differing q, and the Gaussian variances are 


1 
(25) = (#28) = 52 


from which we can immediately construct 


(Ba8-za) = (8a.n) * (88.3) = Her 
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From eq.(II.17) we can calculate the correlations in the phase 0(x) in real space. Clearly 
(0(x)) =0 by symmetry, while 

1 iq:xtiq’-x’ 1 

(O(x)0(x')) = 5 OE ee aCe a 7 S>——. (11.18) 


q,q’ q 


In the continuum limit, the sum can be replaced by an integral ($7, > V f d%q/(27)%), 


and ie. cade 
eq et = Ca(x — x’) 
IN\es a SEN II.1 
(oxo) = f Gea — (1.19) 
The function, 
d?q e' ax 
es = IL.2 
cu == | ayia al 


is the Coulomb potential due to a unit charge at the origin in a d—dimensional space, since 


it is the solution to 


d 2 
V?Cu(x) = / on 7 ea = §4(x). (11.21) 


We can easily find a solution by using Gauss’ theorem, 
/ d’aV°Ca= ¢ dS-VCq 


For a spherically symmetric solution, VCqg = (dCq/dx)#, and the above equation simplifies 


to a0 
1= Sau, (11.22) 
where 
See es 11.23 
is the total solid angle (area of unit sphere) in d dimensions. Hence 
dC a 1 go 4 
peel — C a II.24 
ae Sat! ’ a(x) (2 = d)Sa + Co, ( ) 
where co is a constant of integration. 
The long distance behavior of Cq(x) changes dramatically at d= 2, as 
Co GF 2 
2-4 
Jim Calz) = 4 @—aSq d<2 | (11.25) 
In(z) 
= 2 
20 : 


The constant of integration can obtained by looking at 
([8(x) — 0(x')]?) = 2(8(x)*) — 2(6(x)0(x’)), (II.26) 


which goes to zero as x > x’. Hence, 


or 2 (|x — x’|?-4 — a2-4) 
([A(x) — O(x*)]°) = —KQ-dSa” (1.27) 
where a is of the order of the lattice spacing. 
For d > 2, the phase fluctuations are finite, while they become asymptotically large 
for d < 2. Since the phase is bounded by 27, this implies that long range order in the 
phase is destroyed. This result becomes more apparent by examining the effect of phase 


fluctuations on the two point correlation function 
((x)p*(0)) = G2 (e899 ON) (11.28) 


(Since amplitude fluctuations are ignored, we are in fact looking at a transverse correlation 


function.) We shall prove later on that for any collection of Gaussian distributed variables, 


(exp(a6)) = exp (Se) | 


Taking this result for granted, we obtain 


= it 7 2-4 — 92-4 
(w(x)b*(0)) = ¥? exp — 5 Gx) = H(0))) =~ exp || , (II.29) 
and asymptotically 
tandateas a apr ford >2 _ 
Jim WejuO) = rd>? (11.30) 


The saddle point approximation to the order parameter w, was obtained by ignoring fluctu- 
ations. The above result indicates that inclusion of phase fluctuations leads to a reduction 
of order in d > 2, and its complete destruction in d < 2. 

The above example typifies a more general result known as the Mermin—Wagner theo- 
rem. The theorem states that there is no spontaneous breaking of a continuous symmetry 
in systems with short range interactions in dimensions d < 2. Some corollaries to this 
theorem are: 

(1) The borderline dimensionality of two, known as the lower critical dimension, has to 
be treated carefully. As we shall demonstrate later on in the course, there is in fact 

a phase transition for the two dimensional superfluid, although there is no true long 

range order. 

(2) There are no Goldstone modes when the broken symmetry is discrete (e.g. for n = 1). 

In such cases, long range order is possible down to the lower critical dimension of 

dg= tk. 
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II.D Scattering and Fluctuations 


In addition to bulk thermodynamic experiments, scattering measurements can be used 
to probe microscopic fluctuations at length scales of the order of the probe wavelength 4. 
In a typical set up, a beam of wavevector k; is incident upon the sample and the scattered 
intensity is measured at wavevector k, = k;+q. For elastic scattering, |k;| = |k,| =k, and 
q = |q| = 2ksin@, where @ is the angle between incident and scattered beams. Standard 
treatments of scattering start with the Fermi Golden Rule, and usually lead to a scattering 


amplitude of the form 
A(q) « (k, ® f|U|k; © 7) « o(q) [ aixe!** (0). (11.31) 


In the above expression, |i) and |f) refer to the initial and final states of the microscopic 
scattering element (the atom or ion), U/ is the scattering potential that can be decomposed 
as a sum due to the various scattering elements in the sample. The amplitude has a local 
form factor o(q) describing the scattering from an individual element. For our purposes, 
the more interesting global information is contained in p(q), the Fourier transform of the 
global density of scatterers p(x). The appropriate scattering density depends on the nature 
of the probe. Light scattering senses the actual atomic density, electron scattering measures 
the charge density, while neutron scattering is usually used to probe the magnetization 
density. Most such probes actually do not look at a snapshot of the system, but look at 


time averaged configurations. Thus the observed scattering intensity is 


S(q) x (|A(a)|*)  (lo(q)|’). (11.32) 


Here (e) indicates the thermal average of e, which can be used in place of the time average 
in most cases due to ergodicity. 

Eq.(I1.32) indicates that a uniform density only leads to forward scattering (q = 
0), while the long-wavelength fluctuations can be studied by working at small angles 
or with small k. If scattering is caused by the magnetization density, we can use the 
Landau-—Ginzburg Hamiltonian to compute its intensity. The probability of a particular 


configuration is given by 
K t 
P[m(x)] ox exp - | d¢x F(vmy le 5m + um! \ : (11.33) 
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As discussed earlier, the most probable configuration is uniform, with m(x) = mé1, where 
€, is a unit vector (7m is zero for t > 0, and equal to ,/—t/4u for t < 0). We can examine 
small fluctuations around such a configuration by setting 


n 


T(x) = [M+ be(x)]ér + D> ot,0(x)Eas (II.34) 

a=2 
where dy and ¢; refer to longitudinal and transverse fluctuations respectively. The latter 
can take place along any of the n—1 directions perpendicular to the average magnetization. 
After the substitution of eq.(II.34), the terms appearing in the Landau-Ginzburg 


Hamiltonian can be expanded to second order as 
(Vm)* = (Vge)? + (Ver)?, 
m =m? + 2mobe+ 7 + ¢, 
m* = m* + 4m de + 6762 + 2m79? + O(¢43, 43), 
resulting in a quadratic energy cost 


t+12um? 


BH =-InP =V (5m? + um) - [atx [F(v0) a 5 bj 


F / . Eve : (11.35) 


t+ 4um? 


| + O(¢7, $3). 


For uniform distortions, the longitudinal and transverse restoring potentials have “spring 


constants” given by, 


Ke. _5  0?U(m) t for t > 0 
= =t+12um eer yr milan ee (1.36) 
£L £ m 
and - 520 (m) 
7 m t fort > 0 
ey Ue eerie { 11.37 
a ie Og? |, 0 fort <0 ( ) 


(The physical significance of the length scales €¢ and &; will soon become apparent.) Note 
that there is no distinction between longitudinal and transverse components for the param- 
agnet (t > 0). For the ordered magnet in t < 0, there is no restoring force for the transverse 
fluctuations which correspond to the Goldstone modes discussed in the previous section. 
Following the change of variables to the Fourier modes, $(x) = )/, bge’t*/VV, the 


probability of a particular fluctuation configuration is given by 


P [{¢2,q3 $t,q}] « [Lex -$@? + &*)|beq *} - exp Se + & *\\big . (IL38) 
q 
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Clearly each mode behaves as a Gaussian random variable of zero mean, and the two point 


correlation functions are 
0a,69q,—a! 
K(q? + &2°) 


where the indices refer the longitudinal, or any of the transverse components. By 


(¢a,a%6,q') = (11.39) 


using a spin polarized source of neutrons, the relative orientations can be adjusted 
to probe either the longitudinal or the transverse correlations. The Lorentzian form, 
S(q) «x 1/(q? + €-?), usually provides an excellent fit to scattering line shapes away from 
the critical point. Eq.(II.39) indicates that in the ordered phase, longitudinal scattering 
still gives a Lorentzian form (on top of a delta-function at q = 0 due to the spontaneous 
magnetization), while transverse scattering always grows as 1/q?. The same power law 
decay is also predicted to hold at the critical point, t = 0. Actual experimental fits yield 
a power law of the form 


S(q,T =T.) « ; (11.40) 


qn 


with a small positive value of 1. 


II.E Correlation Functions and Susceptibilities 


We can also examine the extent of fluctuations in real space. The averages (¢a(x)) = 


(Ma(x) — Mq), are clearly zero, and the connected correlation function is 


Go,a(%) x’) =((Mo(X) — Ma) (™ma(x') — Ma) 


a ia-xtiaq’-x’ II.41 
=(alx)oa(x')) = 5 Sel (a.qh3.4") es 
qq’ 
Using eq.(II.39), we obtain 
ba.8 cia: (x—x’) ba,6 
Cc / — ’ ——r ? IT = / ; 
Goa(x, x) V > K(@ eee) K a(X x So) (II 42) 


where in the continuum limit, 


dt iq: 


Alternatively, Ig is the solution to the following differential equation 


; _ d?q ged = d¢q ea? teas: Sat T(x) 
Vv (a) = [ Sh = [ Sh |- sale = 64x) +P. ay 
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The solution is spherically symmetric, satisfying 


ta. a= lily Ta d 

— + —_— = 5+ : 11.45 

dx? my go dae. ~ “62 ee ( ) 

We can try out a solution that decays exponentially at large distances as 
Pj ees) 


geP 


(11.46) 


(We have anticipated the presence of a subleading power law.) The derivatives of Iq are 


given by 
dl 1 
ae —— a + — Ly 
dx 1 ae 
: (11.47) 
Ti (BOE De eet Dig 
dx2 x? aS ad 
For « # 0, the requirement that eq.(II.46) satisfies eq.(II.45) gives 
1 2 1 d—1 d—1 1 
DID 2) BN NO) o (11.48) 


Ae 
The choice of € as the decay length ensures that the constant terms in the above equation 
cancel. The exponent p is determined by requiring the next largest terms to cancel. For 
x < €, the 1/x? terms are the next most important; we must set p(p +1) = p(d —1), 
and p =d-— 2. This is the familiar exponent for Coulomb interactions, and indeed at this 


length scale the correlations don’t feel the presence of €, and decay as 
(ee F).. (11.49) 


(Note that a constant term can always be added to the solution to satisfy the limits 
appropriate to the correlation function under study.) At large distances x >> €, the 1/(x€) 
term dominates eq.(II.48), and its vanishing implies p = (d —1)/2. Matching to eq.(II.49) 


at x & € yields 
EB—-d)/2 
Iy(x) a CE exp (—2/€) (o5%6). (11.50) 
The length scale € is known as the correlation length. From eq.(II.42), we observe that 
transverse and longitudinal correlations behave differently. Close to the critical point, the 


longitudinal correlation length (eq.(II.36)) behaves as 


te i/K for t > 0 
f= (11.51) 
(—2t)-/2/./K fort <0 
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The singularities can be described by €4 ~ 9 B|t|~”+, where vz = 1/2 and B,/B_ = V2 
are universal, while &) « 1//K is not. The transverse correlation length (eq. (II.37)) equals 
€¢ for t > 0, and is infinite for all t < 0. 

Eq.(IL.49) implies that right at T.., correlations decay as 1/r¢~?. Actually, the decay 


exponent is usually indicated by 1/a¢~?+", where 7 is the same exponent introduced in 
eq.(II.40). Integrating the connected correlation functions results in bulk susceptibilities. 


For example, the divergence of the longitudinal susceptibility is also obtained from, 


d?xGe caer vee 11.52 
Xe x xG5(x) x pana * 8 ~ Pe Ae (11.52) 
0 


The universal exponents and amplitude ratios are again recovered from the above equation. 
For T < T,, there is no upper cut-off length for transverse correlations, and the divergence 


of the transverse susceptibility can be related to the system size L, as 
L dy 
Xt x [ aixceen x f =z & oe (11.53) 
o wv 


II.F Comparison to Experiments 
The true test of the validity of the approach outlined in previous sections is in its 
comparison with experiments. A table of exponents, and the appropriate experimental 


materials is provided below. 


Transition Type Material a B y V 
Ferromagnets (n= 3) Fe, Ni -0.1 0.4 1.3 

Superfluid (n = 2) Het 0 0.3 1.3 0.7 
Liquid-Gas (n = 1) CO2, Xe 0.1 0.3 12 0.7 
Ferroelectrics TGS 0 1/2 1 1/2 
and Superconductors 

Mean-field theory 0 1/2 1 1/2 


The exponents are actually known to much better accuracy than indicated in the 
table. The final row (mean-field theory) refers to the results obtained from the saddle point 
approximation. They agree only with the experiments on ferroelectric and superconducting 
materials. The disagreement between the exponents for different values of n suggests that 
the mean-field results are too universal, and leave out some essential dependence on n (and 
d). How do we account for these discrepancies? The starting point of the Landau-Ginzburg 
Hamiltonian is sufficiently general to be trustworthy. The difficulty is in the saddle point 
method used in the evaluation of its partition function, as will become apparent in the 


following sections. 
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II.G Gaussian Integrals 


In the previous section, the energy cost of fluctuations was calculated at quadratic 
order. These fluctuations also modify the saddle point free energy. Before calculating 
this modification, we take a short (but necessary) mathematical diversion on performing 
Gaussian integrals. 

The simplest Gaussian integral involves one variable ¢, 

T= [. dpe ze the = otk, (11.54) 


By taking derivatives of the above expression with respect to h, integrals involving powers 


of ¢@ are generated; e.g. 


d ne 2 Q2n7 12 h 
mee caer —>¢ the _ ——p@ 2K «.— 
Th [« pe ra a 
P 2 5 ; 12 (11.55) 
is. -Kdrtng _ /47 2 |S 
The: [a 6% K. Faroe 


If the integrand represents the probability density of the random variable ¢, the above 
integrals imply the moments (¢) = h/K, and (67) = h?/K?4+1/K. The corresponding 
cumulants are (¢). = (6) =h/K, and (67). = (¢7) — (¢)? =1/K. In fact all higher order 


cumulants of the Gaussian distribution are zero since 


y= (6), = exp [isn = = , (11.56) 


Now consider the following Gaussian integral involving N variables, 


In= [] d¢i exp ye gi: » hid; | . (1.57) 
~ ij 

It can be reduced to a product of N one dimensional integrals by diagonalizing the matrix 

K = K;,;. Since we need only consider symmetric matrices (K;,; = Kj,i), the eigenvalues 

are real, and the eigenvectors can be made orthonormal. Let us denote the eigenvectors 

and eigenvalues of K by g¢ and K, respectively, i.e. Kg = K,q. The vectors {q} form a 

new coordinate basis in the original N dimensional space. Any point in this space can be 


represented either by coordinates {¢;}, or { ,} with @ = >> 5 ba di- We can now change 
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the integration variables from {¢;} to {da}. The Jacobian associated with this unitary 


transformation is unity, and 


se ee Tapes ese NO f Dae h 

q=lv 
The final expression can be represented in terms of the original coordinates by using 
the inverse matrix K~1+, such that K~!K = 1. Since the determinant of the matrix is 


independent of the choice of basis, det K = [[, Aq, and 


at 

(27) Ki 

Lg = “ Ajh;|. II. 

N aK Dy : hih,; (11.59) 
iJ 

Regarding {¢;} as Gaussian random variable distributed with a joint probability distri- 

bution function proportional to the integrand of eq.(II.57), the joint characteristic function 


is given by 


= hes Kk7!} 
(« M8) — exp iS Kj hak; — 5° a ae (11.60) 


4,9 4, 
Moments of the distribution are obtained from derivatives of the characteristic function 


with respect to kj, and cumulants from derivatives of its logarithm. Hence, eq.(II.60) 


implies 
(Pie = yh 
j (11.61) 
(bibj)e = Key 
Another useful form of eq.(II.60) is 
(exp(A)) = exp (ay. + 5(4").| (11.62) 


where A = )}, a;; is any linear combination of Gaussian distributed variables. We used 
this result earlier in computing the order parameter correlations in the presence of phase 
fluctuations in a superfluid. 

Gaussian functional integrals are a limiting case of the above many variable integrals. 


Consider the points 2 as the sites of a d—dimensional lattice and let the spacing go to zero. 


29 


In the continuum limit, {¢;} go over to a function ¢(x), and the matrix Kj; is replaced by 


a kernel K(x,x’). The natural generalization of eq.(II.59) is 


[Docs |- f ataate AS oenace')+ f atarcoots 


es (11.63) 
x (det K)~!/? exp [/ tincdtx’ EE nosin(x’) ; 
where the inverse kernel K~!(x, x’) satisfies 
[atx x!)K 7! (x'!, x”) = 64(x — x"). (11.64) 


The notation D¢(x) is used to denote the functional integral. There is a constant of 
proportionality, (27)%/?, left out of eq.(II.63). Although formally infinite in the continuum 
limit of N — oo, it does not effect the averages that are obtained as derivatives of such 


integrals. In particular, for Gaussian distributed functions, eq.(II.61) generalizes to 


(6(x))e = i dx! K(x, x/)h (x!) 
(b(x)O(x')\e = Kx, x’) 


(11.65) 


In dealing with small fluctuations to the Landau—Ginzburg Hamiltonian, we encoun- 


tered the quadratic form 
[ atxvoy? + 62/22] = f atxatxo(x')34ex — x)(-V? +€7)600, —_(IL66) 
which implies the kernel 
K(x,x’) = K64(x — x')(-V? + €-”). (11.67) 
Following eq.(II.64), the inverse kernel satisfies 
K f atxts(x = x"\(—V? 46°) K(x" — x’) = 54 (x! =x); (II.68) 
which implies the differential equation 
K(-V2 + €-*)K 71 (x) = 64(x). (11.69) 


Comparing with eq.(II.44) implies K~'(x) = (¢(x)¢(0)) = —Ia(x)/K, as obtained before 
by a less direct method. 
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Il.H Fluctuation Corrections to the Saddle Point 


We can now examine how fluctuations around the saddle point solution modify the free 
energy, and other macroscopic properties. Starting with eq.(II.35), the partition function 


including small fluctuations is 


Z = exp [-v (Si + um) [ Pees exp i-5 jetx v6)? + | \ 


Jonwael & fetefowrsg]} 


Each of the Gaussian kernels is diagonalized by the Fourier transforms 


(11.70) 


ia) = / dx exp (~iq- x) 6(0)/VV, 


and with corresponding eigenvalues K(q) = K(q?+€~). The resulting determinant of K 


is a product of such eigenvalues, and hence 


d 
Indet K = “In K(q) = vf ot In[K(q? + €~?)]. (11.71) 


The free energy resulting from eq.(II.70) is then given by 


InZ Ae) 1 dt —1 dé 
peo ae / on n[K(q? + €2)) + — / Gat In[K (q? + &*)]. 


V 2 9 2 
(11.72) 


(Note that there are n—1 transverse components.) Using the dependance of the correlation 


lengths on reduced temperature, the singular part of the heat capacity is obtained as 


dé 1 
043 [oo for t > 0 


oa 2 J (27)4 (Kq? +t)? 
Ge suete 2 ae @ a ) | 1.73) 
q 
aap cs Eee 
Bu lesar casa: oa 


The correction terms are proportional to 


1 d¢q iF 


The integral has dimensions of (length)*~?, and changes behavior at d = 4. For d > 4 
the integral diverges at large q, and is dominated by the upper cutoff A ~ 1/a, where a is 


dl 


the lattice spacing. For d < 4, the integral is convergent in both limits. It can be made 


dimensionless by rescaling q by €~', and is hence proportional to €4~¢. Therefore 


i canoe ford>A4 
(11.75) 


Crr=> : 
Etre ford<4 


In dimensions d > 4, fluctuation corrections to the heat capacity add a constant 
term to the background on each side of the transition. However, the primary form of the 
singularity, a discontinuity in C, is not changed. For d < 4, the divergence of € « t~!/2, 
at the transition leads to a correction term from eq.(II.75) which is more important than 
the original discontinuity. Indeed, the correction term corresponds to an exponent @ = 
(4 —d)/2. However, this is only the first correction to the saddle point result. The diverge 
of Cr merely implies that the saddle point conclusions are no longer reliable in dimensions 
d < 4, called the upper critical dimension. Although we obtained this dimension by 
looking at the fluctuation corrections to the heat capacity, we would have reached the same 
conclusion in examining the singular part of any other quantity, such as magnetization or 
susceptibility. The contributions due to fluctuations always modify the leading singular 


behavior, and hence the critical exponents, in dimensions d < 4. 


II.I_ The Ginzburg Criterion 


We have thus established the importance of fluctuations, and identified them as the 
probable reason for the failure of the saddle point approximation to correctly describe 
the observed exponents. However, as noted in sec.II.G, there are some materials, such as 
superconductors, in which the experimental results are well fitted to the singular forms 
predicted by this approximation. Can we quantify why fluctuations are less important in 
superconductors than in other phase transitions? 

Eq.(II.75) indicates that fluctuation corrections become important due to the diver- 
gence of the correlation length. Within the saddle point approximation, the correlation 
length diverges as € & &9|t|~!/?, where t = (T, — T)/T, is the reduced temperature, and 
&) = WK is a microscopic length scale. In principal, ) can be measured experimentally 
from fitting scattering line shapes. It has to approximately equal the size of the units that 
undergo ordering at the phase transition. For the liquid—gas transition, €) can be estimated 


1/3 where ve is the critical atomic volume. In superfluids, &) is approximately the 


as (Ue) 
thermal wavelength (7). Both these estimates are of the order of a few atomic spacings, 


1-10A. On the other hand, the underlying unit for superconductors is a Cooper pair. The 
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paired electrons are forced apart by their Coulomb repulsion, resulting in a relatively large 
separation of 9 + 10°A. 

The importance of fluctuations can be gauged by comparing the two terms in 
eq.(II.73); the saddle point discontinuity ACs p, « 1/u, and the correction term Cr. 


—(4-d)/2, 


Since K «x €%, the correction term is proportional to € ae Thus fluctuations are 


important provided, 
1 


(E4ACs pF 
The above requirement is known as the Ginzburg criterion. Naturally in d < 4, it is 


654s" > ACsp, => ltl<ter (11.76) 


satisfied sufficiently close to the critical point. However, the resolution of the experiment 
may not be good enough to get closer than the Ginzburg reduced temperature tg. If so, 
the apparent singularities at reduced temperatures t > tg may show saddle point behavior. 
It is this apparent discontinuity that then appears in eq.(II.76), and may be used to self— 
consistently estimate tg. Clearly, ACs.p, and €) can both be measured in dimensionless 
units; €9 in units of atomic size a, and ACs p, in units of Nkg. The latter is of the order of 
unity for most transitions, and thus tg & ° in d= 3. In cases where 9 is a few atomic 
spacings, a resolution of tg ~ 107! — 10~? will suffice. However, in superconductors with 
&) & 10%a, a resolution of tg < 107!8 is necessary to see any fluctuation effects. This 
is much beyond the ability of current apparatus. The newer ceramic high temperature 
superconductors have a much smaller coherence length of 9 ~ 10a, and they indeed show 
some effects of fluctuations. 

Again, it is worth emphasizing that a similar criterion could have been obtained by 
examining any other quantity. Fluctuations corrections become important in measurement 
of a quantity X fort < tg(X) ~ A(X)éy° 4, However, the coefficient A(X) may be 
different (by one or two orders of magnitude) for different quantities. So, it is in principle 
possible to observe saddle point behavior in one quantity, while fluctuations are important 
in another quantity measured at the same resolution. Of course, fluctuations will always 
become important at sufficiently high resolutions. 

A summary of the results obtained so far from the Landau—Ginzburg approach is as 
follows: 

e For dimensions d greater than an upper critical dimension of d,, = 4, the saddle point 
approximation is valid, and singular behavior at the critical point is described by exponents 
a= 0; 0 S17/2,7 = 14 = 1/2,%>. 

e For d less than a lower critical dimensions (dg = 2 for continuous symmetry, and d; = 1 
for discrete symmetry) fluctuations are strong enough to destroy the ordered phase. 

e In the intermediate dimensions, dy < d < d,,, fluctuations are strong enough to change 
the saddle point results, but not sufficiently important to completely destroy order. Un- 
fortunately, or happily, this is the case of interest to us in d= 3. 
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III. The Scaling Hypothesis 


IIIl.A The Homogeneity Assumption 


In the previous chapters, the singular behavior in the vicinity of a continuous transi- 
tion was characterized by a set of critical exponents {a, 3, y,6,v,7,---}. The saddle—point 
estimates of these exponents were found to be unreliable due to the importance of fluctua- 
tions. Since the various thermodynamic quantities are related, these exponents can not be 
independent of each other. The goal of this chapter is to discover the relationships between 
them, and to find the minimum number of independent exponents needed to describe the 
critical point. 

The non-analytical structure is a coexistence line for t < 0 and h = 0, that terminates 
at the critical point t = h = 0. The various exponents describe the leading singular be- 
havior of a thermodynamic quantity Q(t, h), in the vicinity of this point. A basic quantity 


in the canonical ensemble is the free energy, which in the saddle—point approximation is 


given by 
t? 
‘ —— for = 0-20 
fh) = min 5m tum*—hm| « ays : (III.1) 
ae 7173 for h Ea 0, t=0 
u 


The singularities in the free energy can in fact be described by a single homogeneous 


function’ in ¢ and h, as 
f(t,h) = |t)?gp (h/lel>) - (III.2) 


The function gr only depends on the combination x = h/|t|*, where A is known as the gap 
exponent. The asymptotic behavior of gf is easily obtained by comparing eqs.(III.1) and 
(III.2). The h = 0 limit is recovered if lim;-.9 g(x) ~ 1/u, while to get the proper power 
of h, we must set lim, g(x) ~ x4/3/u!/3. The latter implies f ~ |t|?h4/3/(u/3|t/44/9). 


* In general, a function f (a1, 22,---) is homogeneous if 
f (GPa, OF ea) = OFT ig Bay"), 


for any rescaling factor b. With the proper choice of b one of the arguments can be removed, 


leading to a scaling forms used in this section. 
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Since f can have no t dependence along t = 0, the gap exponent (corresponding to 
Eq.(III.1) has the value 
A= ; (III.3) 
2 
The assumption of homogeneity is that, on going beyond the saddle—point approxima- 
tion, the singular form of the free energy (and any other thermodynamic quantity) retains 
the homogeneous form 
feing.(t, 2) = |t|? “gy (h/It|>) - (III.4) 


The actual exponents a and A depend on the critical point being considered. The depen- 
dence on t is chosen to reproduce the heat capacity singularity at h = 0. The singular part 
of the energy is obtained from (say for t > 0) 


of 


ar (2—atl%gz (h/t) — Ant-9-4 9}, (n/lel4) 


Esing. ot 
ne (IIL.5) 
~ 12 — agg (h/|t|4) — ya if (h/|t|*) | = tC ge (h/|t\*). 


Thus the derivative of one homogeneous function is another. Similarly, the second deriva- 


tive takes the form (again for t > 0) 


0? 4 
Cris: oa at ~t "9c (h/|t|*) ; (III.6) 


Ot 


reproducing the scaling Cying, ~ (%, as h — 0. 


It may appear that we have the freedom to postulate a more general form 


Ce nyalil- "ge (hflel"*)., (III.7) 


with different functions and exponents for t > 0 and t < 0, that match at t = 0. However, 
this is ruled out by the condition that the free energy is analytic everywhere except on the 
coexistence line for h = 0 and t < 0, as proven as follows: Consider a point at t = 0 and 
finite h. By assumption, the function C’ is perfectly analytic in the vicinity of this point, 


expandable in a taylor series, 
C(t « h4) = A(h) + tB(h) + O(2?). (III.8) 


Furthermore, the same expansion should be obtained from both C, and C_. But eq. (III.7) 


leads to the expansions, 


Cx = [t\-4 [As (Ae) 3) (4) saa | (IIL.9) 
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where {p 


Q 


} are the leading powers in asymptotic expansions of g+ for large arguments, 


and {A+, B+} are the corresponding pre-factors. Matching to the taylor series in eq. (III.8) 


requires p,A 


Continuity at 


= —a, and qzA = —(1+ ay), and leads to 


Cy (t < AA) = AghO#/4t 4 Buh OC ee)/Ak |e] pee, (III.10) 


t = 0 now forces a4/A, = a_/A_, and (1+a,)/Ay = (l+a_)/A_, 


which in turn implies 


Despite using 


a+ =AaA_~_ =a 
(111.11) 
A, =A_=A 


|t] in the postulated scaling form, we can still ensure the analyticity of 


the function at t = O for finite h by appropriate choice of parameters, e.g. by setting 
B_ = —B, to match Eq.(III.10) to the analytic form in Eq.(III.8). Having established 


this result, we can be somewhat careless henceforth in replacing |t| in the scaling equations 


with t. Naturally these arguments apply to any quantity Q(t, h). 


Starting from the free energy in eq.(III.4), we can compute the singular parts of other 


quantities of interest: 


e The magnetization is obtained from 


m(t, h) ~ or patel gaa ay lel ys (III.12) 


In the limit x — 0, gm(x) is a constant, and 


m(t,h =0) ~ |t?-9-4, —> B=2-a—-A. (III.13) 


On the other hand, if x — 00, gm(x) ~ x?, and 


Pp 
m(t = 0,h) ~ |t|2-2-4 (=) (111.14) 


Since this limit is independent of t, we must have pA = 2—a-— A. Hence 


m(t,h =0)~ An@-e-AV/A, => = § =A/(2—-a—A)=A/f. (111.15) 


e Similarly, the susceptibility is computed as 


Om 
x(t, h) ~ Oh 


~ |t?-9 74g, (h/EA), => x(t,h = 0) ~ |e)? 8-74, > y= 2A-2+4+0. 
(III.16) 
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Thus, the consequences of the homogeneity assumption are: 

(1) The singular parts of all critical quantities Q(t,h), are homogeneous, with the same 
exponents above and below the transition. 

(2) Because of the interconnections via thermodynamic derivatives, the same gap expo- 
nent A, occurs for all such quantities. 

(3) All (bulk) critical exponents can be obtained from only two independent ones, e.g. a 
and A. 

(4) As a result of the above, there are a number of exponent identities. For example, 
eqs.(III.13), (III.15), and (III.16) imply 


a+28+7y=a+2(2-—a—A)4+(2A—2+a) =2 (Rushbrooke’s Identity), 
6-1l= — —-l= an = 5 (Widom’s Identity). 
(111.17) 
These identities can be checked against the following table of critical exponents. The 
first three rows are based on a number of theoretical estimates in d = 3; the last row comes 
from an exact solution in d = 2. The exponent identities are completely consistent with 


these values, as well as with all reliable experimental data. 


a B a7 6 V n 
n=1 0.11 0.32 1.24 4.9 0.63 0.04 
n=2 -0.01 0.35 1.32 4.7 0.67 0.04 
n=3 -0.11 0.36 1.39 4.9 0.70 0.04 
n=1 0 1/8 7/4 15 1 1/4 


III.B Divergence of the Correlation Length 


The homogeneity assumption relates to the free energy and quantities derived from 
it. It says nothing about the behavior of correlation functions. An important property of 
a critical point is the divergence of the correlation length, which is responsible for, and can 
be deduced from, diverging response functions. In order to obtain an identity involving 
the exponent v for the divergence of the correlation length, we replace the homogeneity 
assumption for the free energy, with the following two conditions: 


(1) The correlation length € has a homogeneous form, 


E(t, h) ~ |t|~"g (h/t!) . (II1.18) 
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(For t = 0, € diverges as h~”" with v;, = v/A.) 

(2) Close to criticality, the correlation length € is the most important length in the system, 
and is solely responsible for singular contributions to thermodynamic quantities. 
The second condition determines the singular part of the free energy. Since In Z(t, h) 


is extensive and dimensionless, it must take the form 


ae LN 
InZ = (=) X gst (=) Ga. 2s (III.19) 
a 


where g, and g, are non-singular functions of dimensionless parameters (a is an appropriate 
microscopic length). The singular part of the free energy comes from the first term, and 


behaves as 


InZ _ 7 
fied ht) re 1m lt| gp (h/|t|4) . (III.20) 


A simple interpretation of the above result is obtained by dividing the system into units 
of the size of the correlation length. Each unit is then regarded as an independent random 
variable, contributing a constant factor to the critical free energy. The number of units 
grows as (L/€), leading to eq.(III.19). 

The consequences of the above assumptions are: 
(1) The homogeneity of fsing.(t, h) emerges naturally. 


(2) We obtain the additional exponent relation 
2—a=dv — (Joshephson’s identity) . (III.21) 


Identities obtained from the generalized homogeneity assumption involve the space 
dimension d, and are known as hyperscaling relations. The relation between a and v is 
consistent with the exponents in the above table. However, it does not agree with the 
saddle point values, a = 0 and v = 1/2, which are valid for d > 4. Any theory of critical 
behavior must thus account for the validity of this relation in low dimension, and its 


breakdown for d > 4. 
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III.C Critical Correlation Functions and Self-Similarity 


One exponent that has not so far been accounted for is 7, describing the decay of 
correlation functions at criticality. Exactly at the critical point, the correlation length 
is infinite, and there is no other length scale (except sample size) to cut off the decay 
of correlation functions. Thus all correlations decay as a power of the separation. As 
discussed in the previous chapter, the magnetization correlations fall off as 


Cain (X) = (n(x) (0) = (mm)? ~ 1/[x|*-24, (111.22) 
Similarly, we can define an exponent 1 for the decay of energy—energy correlations as 
p(X) = (H(x)YH(0)) — (H)? ~ 1/|x|OP (11.23) 


Away from criticality, the power laws are cut off for distances |x| >> €. As the response 
functions can be obtained from integrating the connected correlation functions, there are 
additional exponent identities, such as (Fisher’s identity) 


§ day 
xX~ [aim ~ ij corn w EM lt] -VA-M yy = (2— .. (TIL-24) 
Similarly, for the heat capacity, 


g d 
Cw [axe ~ i; re wT ny eV), sy = (2-1 
(111.25) 
As before, two independent exponents are sufficient to describe all singular critical behavior. 
An important consequence of these scaling ideas is that the critical system has an 
additional dilation symmetry. Under a change of scale, the critical correlation functions 
behave as 

Gaia Ax) =A Gatieaile): (III.26) 
This implies a scale invariance or self-similarity: if a snapshot of the critical system 
is blown up by a factor of A, apart from a change of contrast (multiplication by X”), the 
resulting snapshot is statistically similar to the original one. Such statistical self-similarity 
is the hallmark of fractal geometry. As discussed by Mandelbrot, many naturally occurring 
forms (clouds, shore-lines, river basins, etc.) exhibit such behavior. The Landau—Ginzburg 
probability was constructed on the basis of local symmetries such as rotation invariance. If 
we could add to the list of constraints the requirement of dilation symmetry, the resulting 
probability would indeed describe the critical point. Unfortunately, it is not possible to 
directly see how such a requirement constrains the effective Hamiltonian. One notable 
exception is in d = 2, where dilation symmetry implies conformal invariance, and a lot 
of information can be obtained by constructing conformally invariant theories. We shall 
instead prescribe a less direct route of following the effects of the dilation operation on the 

effective energy; the renormalization group procedure. 
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IlI.D The Renormalization Group (Conceptual 


Success of the scaling theory in correctly predicting various exponent identities 
strongly supports the assumption that close to the critical point the correlation length 
€, is the only important length scale, and that microscopic length scales are irrelevant. 
The critical behavior is dominated by fluctuations that are self-similar up to the scale €. 
The self-similarity is of course only statistical, in that a magnetization configuration is 
generated with a weight W |m(x)] « exp{—GH|m(x)]}. Kadanoff suggested taking advan- 
tage of the self-similarity of the fluctuations to gradually eliminate the correlated degrees 
of freedom at length scales x < €, until one is left with the relatively simple, uncorrelated 
degrees of freedom at scale €. This is achieved through a procedure called the renor- 
malization group (RG), whose conceptual foundation is the three steps outlined in this 


section. 


(1) Coarse Grain: There is an implicit short distance length scale a, for allowed variations 
of (x) in the system. This is the lattice spacing for a model of spins, or the coarse graining 
scale that underlies the Landau-Ginzburg Hamiltonian. In a digital picture of the system, 
a corresponds to the pixel size. The first step of the RG is to decrease the resolution by 
changing this minimum scale to ba (b > 1). The coarse-grained magnetization is then 
given by 

7m; (x) = = : d¢x'm;(x’). (III.27) 

U8 Jcaueenerea atx 

(2) Rescale: Due to the change in resolution, the coarse grained ‘picture’ is grainier than 
the original. The original resolution of a can be restored by decreasing all length scales by 
a factor of b, i.e. by setting 


Xnew = aa (111.28) 


(3) Renormalize: The variations of fluctuations in the rescaled magnetization profile is 
in general different from the original, i.e. there is a difference in contrast between the 
pictures. This can be remedied by introducing a change of contrast by a factor ¢, through 


defining a renormalized magnetization 


1 


Bar [ : ia d¢x'm”(x’). (III.29) 
ell centered at bXnew 


Wide (new) 


By following these steps, for each configuration 77ioiq(x), we generate a renormalized 
configuration Mnew(x). Eq.(III.29) can be regarded as a mapping from one set of random 


variables to another, and can be used to construct the probability distribution, or weight 


Al 


W,|new(x)] = exp{—CHo|i7inew(x)|}. Kadanoff’s insight was that since on length scales 
less than €, the renormalized configurations are statistically similar to the original ones, 
they may be distributed by a Hamiltonian GH, that is also ‘close’ to the original. In 
particular, the original Hamiltonian becomes critical by tuning the two parameter t and 
h to zero, at which point the dominant configurations are similar to those of the rescaled 
system; the critical Hamiltonian is thus invariant under such rescaling. In the original 
problem, one moves away from criticality for finite t and h. Kadanoff’s assumption is that 
the corresponding new Hamiltonian is also described by non-zero tpew Or Anew- 

The assumption that the vicinity of the original and renormalized Hamiltonians to 
criticality is described by the two parameters t and h greatly simplifies the analysis. The 
effect of the RG transformation on the probability of configurations is now described by 
the two parameter mappings thew = to(toia, Pola) aNd Rnew = Ao(toia, hoa). The next 
assumption is that since the transformation only involves changes at the shortest length 
scales, it cannot cause any singularities. The renormalized parameters must be analytic 


functions of the original ones, and hence expandable as 


( ty(t,h) = A(b)t + B(b)ht+--- (III.30) 


hy(t,h) = C(b)t + D(b)h+--: 


Note that there are no constant terms in the above Taylor expansions. This expresses the 
condition that if GH is at its critical point (t = h = 0), then 3H, is also at criticality, and 
tnew = Anew = 0. Furthermore, due to rotational symmetry, under the combined transfor- 
mation (m(x) > —m(ax),ht> —h,t + t) the weight of a configuration is unchanged. As 
this symmetry is preserved by the RG, the coefficients B and C in the above expression 


must be zero, leading to the further simplifications 
(111.31) 


The remaining coefficients A(b) and D(b) depend on the (arbitrary) rescaling factor 
b, and trivially A(1) = D(1) = 1 for b = 1. Since the above transformations can be 
carried out in sequence, and the net effect of rescalings of b; and bz is a change of scale 
by 6,b2, the RG procedure is sometimes referred to as a semi-group. The term applies 
to the action of RG on the space of configurations: each magnetization profile is mapped 
uniquely to one at larger scale, but the inverse process is non-unique as some short scale 


information is lost in the coarse graining. (There is in fact no problem with inverting the 


42 


transformation in the space of the parameters of the Hamiltonian.) The dependence of 
A and D in eqs.(III.31) on b can be deduced from this group property. Since at b = 1, 
A= D=1, and t(bib2) & A(b,)A(b2)t & A(b1b2)t; we must have A(b) = bY, and similarly 
D(b) = 6", yielding 
{ =t, =oMt+--- 
(111.32) 
h’ = hyp = DY h+--- 
If GHoia is slightly away from criticality, it is described by a large but finite correlation 
length €oia. After the RG transformation, due to the rescaling in eq.(III.28), the new 
correlation length is smaller by a factor of b. Hence the renormalized Hamiltonian is less 
critical, and the RG procedure moves the parameters further away from the origin, i.e. the 
exponents y; and y, must be positive. 


e We can now explore some consequences of the assumptions leading to eq.(III.32). 


1. The free energy: The RG transformation is a many to one map of the original configu- 
rations to new ones. Since the weight of a new configuration, W'([m’]), is the sum of the 


weights W({m]), of old configurations, the partition function is preserved, i.e. 
L= [ Pew (in) = [ Pmt (nt) = Ae (I1I.33) 
Hence In Z = In Z’, and the corresponding free energies are related by 
VEG Ra Vie By (III.34) 
In d-dimensions, the rescaled volume is smaller by a factor of b7, and 
f(t, h) = b-7f (bt, bh), (111.35) 


where we have made use of the assumption that two free energies are obtained from 
the same Hamiltonian in which only the parameters t and h have changed according to 
eqs. (III.32). Eq.(III.35) describes a homogeneous function of t and h. This is made appar- 
ent by choosing a rescaling factor b such that bY is a constant, say unity, ie. b= t7!/¥%, 
leading to 

f(t,h) =t/ f(A, h/te/) = tg p(h/t4e/M*), (111.36) 


We have thus recovered the scaling form in eq.(III.4), and can identify the exponents as 
2-a=d/y : A= yn/Yt- (III.37) 
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2. Correlation length: All length scales are reduced by a factor of b during the RG 


transformation. This is also true of the correlation length, €’ = €/b, implying 
E(t, h) = bE (bt, DY h) = 7 1/ME(1, h/tr/H) EO”, (111.38) 


This identifies vy = 1/y:, and using eq.(III.37) the hyperscaling identity, 2 — a = dv, is 
recovered. 

3. Magnetization: From the homogenous form of the free energy (eq.(III.36)), we can 
obtain other bulk quantities such as magnetization. Alternatively, from the RG results for 
Z, V, and h, we may directly conclude 

10MmZ(t,h) _ 1 OlnZ'(t',h’) 


i ala eer ee 1 SY 


= bYr—4m(bY*t, bY h). (III.39) 


Choosing b = t~!/¥, we obtain 3 = (yp, — d)/y:, and A = yp/y as before. 
It is thus apparent that quite generally, the singular part of any quantity X has a 
homogeneous form 


X(t, h) = bY* X(bYt, bY*h). (III.40) 


For any conjugate pair of variables, contributing a term [ d’xF - X, to the Hamiltonian, 


the scaling dimensions are related by yx = yr — d, where F’ = bY" F under RG. 


IlIl.E The Renormalization Group (Formal 


In the previous section we noted that all critical properties can be obtained from the 
recursion relations in eqs.(III.32). Though conceptually appealing, it is not clear how such 
a procedure can be formally carried out. In particular, why should the forms of the two 
Hamiltonians be identical, and why are two parameters t and h sufficient to describe the 
transformation? In this section we outline a more formal procedure for identifying the 
effects of the dilation operation on the Hamiltonian. The various steps of the program are 


as follows: 


(1) Start with most general Hamiltonian allowed by symmetries. For example, in the 


presence of rotational symmetry, 


t K 
ou = | atx Em? ut uml 4 


(Vm)? + 5 (V2m)? fee]. (III.41) 


A particular system with such symmetry is therefore completely specified by a point in 


the (infinite-dimensional) parameter space S = (t,u,v,---,K,L,---). 
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(2) Apply the three steps of renormalization in configuration space: (i) Coarse grain by b; 


(ii) Rescale, x’ = x/b; and (iii) Renormalize, m’ = m/¢. This defines a change of variables, 


A) 


m'(x’) = Cpa 
c Cell of size b centered at bx’ 


d’xm(x). (111.42) 
Given the probabilities P[m(x)] « exp(—GH|m(x)]), for the original configurations, we can 
use the above change of variable to construct the corresponding probabilities P’[m’(x’)], 


for the new configurations. Naturally this is the most difficult step in the program. 


(3) Since rotational symmetry is preserved by the RG procedure, the rescaled Hamiltonian 
must also be described by a point in the parameter space of eq.(III.41), ice. 
BH'|m'(x')] = In P[m'(x’)] 
t! K! L’ III.43 
aa fot f atx’ En” + u'm'* + ym’ oh os ba Se = (Vm')’ a 3 (V?m')? SB on htel he ( ) 
The renormalized parameters are functions of the original ones, i.e. t/ = t,(t,u,...); u’ = 


up(t, u,...), etc., defining a mapping S’ = ¥,S in parameter space. 


(4) The operation #, describes the effects of dilation on the Hamiltonian of the system. 
Hamiltonians that describe statistically self-similar configurations must thus correspond 
to fixed points S*, such that #,S* = S*. Since the correlation length, a function of 
Hamiltonian parameters, is reduced by 6 under the RG operation (i.e. €(S) = b&(I,S)), 
the correlation length at a fixed point must be zero or infinity. Fixed points with €* = 
O describe independent fluctuations at each point and correspond to complete disorder 
(infinite temperature), or complete order (zero temperature). A fixed point with €* = co 


describes a critical point (T = T,). 


(5) Eqs.(III.32) represent a simplified case in which the parameter space is two dimensional. 
The point t = h = 0 isa fixed point, and the lowest order terms in these equations describe 
the behavior in its neighborhood. In general, we can study the stability of a fixed point by 
linearizing the recursion relations in its vicinity: Under RG, a point S*+0S is transformed 
to 


_ O86 


S* +65) = S*+(RE)apbSgt+---, where (RE)yg= | . 
03. | ¢- 


(III.44) 
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We now diagonalize the matrix (R/)ag to get the eigenvectors O;, and corresponding 


eigenvalues \(b);. Because of the group property’, 
REREO; = r(d)A(B');O; = REO; = Abb’) Oj. (111.45) 
Together with the condition \(1); = 1, the above equation implies 
A(b); = 0": (111.46) 


The vectors O; are called scaling directions associated with the fixed point S*, and 
y; are the corresponding anomalous dimensions. Any Hamiltonian in the vicinity of the 
fixed point is described by a point S = S*+4;9;O0;. The renormalized Hamiltonian has to 
interaction parameters S’ = S* + ¥j9;b¥‘O;. The following terminology is used to classify 
the eigen-operators: 
e If y; > 0, g; increases under scaling, and O; is a relevant operator. 
e If y; <0, g; decreases under scaling, and O; is an irrelevant operator. 
e If y; =0, g; is called a marginal operator, and higher order terms are necessary to track 
its behavior. 

The subspace spanned by the irrelevant operators is the basin of attraction of the fixed 
point S*. Since € always decreases under RG, and €(S*) = oo; then € is also infinite for 
any point on the basin of attraction of S*. For a general point in the vicinity of S*, the 


correlation length satisfies 


(91, 92,°++) = DE(bY g1, bY? ga, -- +). (111.47) 


For a sufficiently large 6, all the irrelevant operators scale to zero. The leading singularities 
of € are then determined by the remaining set of relevant operators. In particular if 
the operators are indexed in order of decreasing dimensions, we can choose b such that 


b¥1g, = 1. In this case, eq.(II1.47) implies 


€(91,.92-°°-) = 91 F(g2/g8/,---). (III.48) 


* The group property RERE = KE, = RERV, also implies that the linearized matrices 
for different b commute. It is thus possible to diagonalize them simultaneously, implying 


that the eigenvectors {O;} are independent of b. 
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We have thus obtained an exponent 11 = 1/y1, for the divergence of €, and a generalized 
set of gap exponents A, = Ya/y1, associated with gy. 

Let us imagine that the fixed point S* describes the critical point of the magnet in 
eq.(III.41) at zero magnetic field. As the temperature, or some other control parameter, 
is changed, the coefficients of the Hamiltonian are altered, and the point S moves along a 
trajectory in parameter space. Except for a single point (at the critical temperature) the 
magnet has a finite correlation length. This can be achieved if the trajectory taken by the 
point S intersects the basis of attraction of S* only at one point. To achieve this, the basin 
of attraction must have co-dimension one, i.e. the fixed point S* must have one and only 
one relevant operator. This provides an explanation of universality, in that the very many 
microscopic details of the system make up the huge space of irrelevant operators comprising 
the basin of attraction. In the presence of a magnetic field, two system parameters must 
be adjusted to reach the critical point ([. = T. and h = 0). Thus the magnetic field 
corresponds to an additional relevant operator at S*. Again, other ‘odd’ interactions, such 
as {m?,m°,---} should not lead to any other relevant operators. 

Although the formal procedure outlined in this section is quite rigorous, it suffers from 
some quite obvious shortcomings: How do we actually implement the RG transformations 
of step (2) analytically? There are an infinite number of interactions allowed by symmetry, 
and hence the space of parameters S, is inconveniently large. How do we know a priori that 
there are fixed points for the RG transformation; that Jt, can be linearized; that relevant 
operators are few, etc.? Following the initial formulation of RG by Kadanoff, there was a 
period of uncertainty until Wilson showed how these steps can be implemented (at least 


perturbatively) in the Landau-Ginzburg model. 


AT 


III.F The Gaussian Model (Direct solution) 


The RG approach will be applied to the Gaussian model in the next section. For 
the sake of later comparison, here we provide the direct solution of this problem. The 
Gaussian model is obtained by keeping only the quadratic terms in the Landau—Ginzburg 


expansion. The resulting partition function is 
= Oe cee 2, Le, \2 rat 
Z= | Dm(x)exp;— | d°*x gm + > (vin) + rad m)"+---—h-m| >. (II.49) 


Clearly the model is well defined only for t > 0, since there is no m* term to insure its 
stability for t < 0. The partition function still has a singularity at t = 0, and we can 
regard this as representing approaching a phase transition from the disordered side. 

The quadratic form is easily evaluated by the usual rules of Gaussian integration. The 
kernel is first diagonalized by the Fourier modes: The allowed q values are discretized in 
a finite system of size L, with spacing of 27/L. The largest q are limited by the lattice 
spacing, and confined to a Brillouin Zone whose shape is determined by the underlying 
lattice. We shall in fact use a slightly different normalization for the Fourier modes, and 


keep careful track of the volume factors, by setting 
m(q) = | d'xe'**ri(x) 


gies ae (III.50) 
fi(x) = nila) = f Senta) 


(We should really use a different symbol, such as m;(q) to indicate the Fourier modes. 
For the sake of brevity we use the same symbol, but explicitly include the argument q as 
the indicator of the Frourier transformed function.) The last transformation applies to the 
infinite size limit (Z — oo), and V is the system volume. 

In reexpressing the Hamiltonian in terms of Fourier modes, we encounter expressions 


such as 


e7ilata’)-x 2 


[ atxincx)? = [exe ail) -m(q’) = ye (111.51) 


The last expression follows from the vanishing of the integral over x unless q+ q/ = 0, in 


which case it equals the system volume. Similar manipulations leads to the Hamiltonian 


Koz + Lot +--- et 
on = Sy (FREE EE SY (ey? —fi-milq=0). (aL) 
q 
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With the choice of the normalization in eq.(III.50), the Jacobian of the transformation to 


Fourier modes is 1/\/V per mode, and the partition function equals 


24 Dot... Ay 
Z= [v-"? f arica) exp [A ima) +h-m(q = 0) 2» 1-53) 
qa 


The integral for q = 0 is 


Zy=Vv-r? i. din(0) exp sr Im(o)/ + i mio) - (22) ex eS . (111.54) 


After performing the integrations for q 4 0, we obtain 


Vi? 20 nye 
Ge — —___ —_____ : II. 
|r ll eee ee) oe 


The total number of modes N, equals the number of original lattice points. Apart from a 
constant factor resulting from (27)"%/?, the free energy is 
nZ on 


dt 
f(t) = =f, Tagya ln (t+ Ka + Lat ++.) — 5p. (II1.56) 


The integral in eq.(III.56) is over the Brillouin zone, which for a hypercubic lattice of 
spacing a, is a cube of side 27/a centered on the origin. However, we expect the singularities 
to originate from the long wavelength modes close to q = 0. The contributions from the 
vicinity of the Brillouin zone edge are clearly analytic, since the logarithm can be simply 
expanded in powers of ¢ for a finite q?. Thus to simplify the extraction of the singular 
behavior as t — 0, we approximate the shape of the Brillouin zone by a hypershpere of 
radius A = 7 /a. The spherical symmetry of the integrand then allows us to write 

" A 
foing.(t,h) = 5 Ka i: dqq’'In (t+ Kq? + Lq* +---) -— = (111.57) 


where Kg = Sq/(27)¢, and Sq is the d-dimensional solid angle. The leading dependence 


of the integral on t can be extracted after rescaling q by a factor of \/t/K, to 


n ey eae daa 2 4/42 h? 
fine ay S aKa (=) : dxx [In¢ + In (l+¢a + Lta*/K +---)| = ee 
(111.58) 


A9 


Ignoring analytic contributions in t, the leading singular dependence of the free energy can 


be written as 43 
fang (th) =~ [A+ reg] =P tap(h/t4) (IL59) 


Note that the higher order gradients, i.e. terms proportional to L,---, do not effect the 
singular behavior in eq.(III.59). On approaching the point h = 0 for t = 0*, the singular 


part of the free energy is described by a homogeneous scaling form, with exponents 
a, =2-—d/2 ; A=1/2+4d/4. (111.60) 


Since the ordered phase for t < 0 is not stable, the exponent @ is undefined. The suscep- 
tibility y « 0? f/Oh? « 1/t, however, diverges with the exponent y; = 1. 


III.G_ The Gaussian Model (Renormalization Group) 


The RG of the Gaussian model is most conveniently performed in terms of the Fourier 
modes. The goal is to evaluate the partition function 


A oad 2 44... > 
Dik [pn exp - f ee Ae) lm(q)|2 + ino) , (IL61) 


indirectly via the three steps of RG. Note that the Brillouin zone is approximated by a 
hypersphere of radius A. 
1. Coarse-Grain: Eliminating fluctuations at scales a < x < ba; is similar to removing 


Fourier modes with wavenumbers A/b < q < A. We thus break up the momenta into two 


subsets, 
{m(q)} = {4(a7)} © {m(q*)}, (III.62) 
and write 
Z= i Din(q<) i) Da(q> )e~ Pe l.4, (111.63) 
Since the two sets of modes are decoupled in the Gaussian model, the integration is trivial, 
and 
A 4d 
d 
Z ~exp -sv — nt ke eg te) x 
2) Jap (27) 
rae : ‘ (III.64) 
x, dig: ft Kg Lg pee \ 55, + 
[ Ditajexp -{ aq | ———3——_ ] lrn(a) ? +h m0) | 
0 (27) 2 
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2. Rescale: The partition function for the modes m(q<) is similar to the original, except 
that the upper cut-off has decreased to A/b, reflecting the coarsening in resolution. The 
rescaling x’ = x/b in real space, is equivalent to q’ = bq in momentum space, and restores 


the cut-off to its original value. The rescaled partition function is 


Z = eVohl®) x / Diii(q/) 
(III.65) 


A ada! —2,/2 —4 14 

d t+ kb Lb tee > x 

exp -{ oO ee sce a a ec Imn(q’)|? +R- (0)| . 
(Qa )4 2 

3. Renormalize: The final step of RG in real space is the renormalization of magnetization, 


mi! (x!) = m(x’)/¢. Alternatively, we can renormalize the Fourier modes according to 


m'(q!) = m(q’)/z, resulting in 


A dal —2_/2 —4 14 
d t+ Kb Lb ee > 
exp - f q —d 2 ( 2. q 5 q a ) |m!(q’)|? + zh - mm (0) 
(111.66) 
(Note that the factors ¢ and z, for rescalings of magnetization in real and Fourier space, 
are different.) 
Eq.(III.66) indicates the renormalized modes are distributed according to a Gaussian 


Hamiltonian, with renormalized parameters, 


t! = 27b-%¢ 

h' = zh 
Ks bre AR (III.67) 
dc ame’ 8 


The singular point t = h = 0 is mapped onto itself as expected. To make the fluctuations 
scale invariant at this point, we must ensure that the remaining Hamiltonian, stays fixed. 


This is achieved by the choice of z = b!+¢/2, which sets K’ = K, and makes the remaining 


parameters L,--- scale to zero. Away from criticality, the two relevant directions now scale 
as 
t' = bt ye = 2 
———— ; (111.68) 
h! = p't4/2p yn =14+d/2 
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Using the results of sec.(III.D), we can identify the critical exponents, 


y=1/y,=1/2, 

1 
ait EME eM e 
Ut 2 2 4 


a=2-—dvy=2-d/2, 


in agreement with the direct solution in the previous section. 


The fixed point Hamiltonian (t* = h* = L* =---=0), in real space is 
* K d 2 
BH = ri d°x(Vm)°*. (III.69) 


(The subscript a is placed on the integral as a reminder of the implicit short distance 
cut-off.) Under a simple rescaling x +> x’, and m™(x) > ¢m'(x'); K 6 K!' = b4-2¢7K. 
Scale invariance is achieved with the choice ¢ = b!~4¢/?. Forgetting the coarse-graining 


step, a general power of (x), added as a small perturbation to (GH)*, behaves as 
BH" + un / d?xm™ ~~ BH* +un,b4¢” / d?x'm!™, (III.70) 


suggesting that such perturbations scale as 


= opr), => «yn =n-d (5 = 1) : (III.71) 
The values y; = 1+ d/2, and y2 = 2, reproduce the exponents for y, and y; in eq.(III.68). 
The operators with higher powers are less relevant. The next most important operator in 
a system with spherical symmetry is y4 = 4 — d, which is irrelevant for d > 4 but relevant 
for d < 4; ye = 6 — 2d is relevant only for d < 3. Indeed the majority of operators are 


irrelevant at the Gaussian fixed point for d > 2. 
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IV. Perturbative Renormalization Group 


IV.A_ Expectation Values in the Gaussian Model 


Can we treat the Landau-Ginzburg Hamiltonian, as a perturbation to the Gaussian 


model? In particular, for zero magnetic field, we shall examine 
= d K ey ee dy, 4 
BH = BHot+U = | d°*x am a 5 (vm) + 5lV m)* +-+-|+u | d®xm*+---. (IV.1) 


The unperturbed Gaussian Hamiltonian can be decomposed into independent Fourier 


modes, as 


1 t+K@+Ll¢t+-:: 2 d4q t+K@+L¢+-:: 
> rm) ? = = 


RE = 2 (Qn)4 2 


The interaction mixes up the normal modes, and 


U=u f dtxm(xyt ++. 


d. jd. jd Jd 
i poe qid“qod*q3d M4 ix-(qi+q2+q3+q,) 


(Qn) 44 Ma(Q1)Ma(a2)mMe(a3)me(aa) 
(IV.3) 
where summation over a and (3 is implicit. The integral over x sets qi + q2 +q3+ qu = 0, 
and 
d*qid"qod“qz 
Uu= [Ena lanymalas)malas)ma(-a1 —qo—q3)+---. (IV.4) 
From the variance of the Gaussian weights, the two-point expectation values in a 


finite sized system with discretized modes are easily obtained as 


dq -—q 00,8 V 


? 


So —_—_.. IV.5 
t+Kq?+Lq*+-:- ee) 


(mo(q) ma(q‘))o = 


In the limit of infinite size, the spectrum becomes continuous, and eq.(IV.5) goes over to 


ba,8(27) 464 (q + q’) 


; IV.6 
t+K@et+Llg +. ae) 


(Mo(q) me(Q'))o = 
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The subscript 0 is used to indicate that the expectation values are taken with respect to 
the unperturbed (Gaussian) Hamiltonian. Expectation values involving any product of 


m’s can be obtained starting from the identity 


(ox > on ) = exp ys = (miM)o| > (IV.7) 
a 0 


Uj 


which is valid for any set of Gaussian distributed variables {m;}. (This is easily seen by 

‘completing the square.’) Expanding both sides of the equation in powers of {a;} leads to 

aia; 
2 


(mim;)o + 


aiajak Aj,AjAKA] 
6 24 


((mim;)o(MeMI)o + (MEME) o(M;M1)0 + (MEME) o(M;M21) 0) 


1+ a(mi)o + (mim;)o + (mum;mr)o + (mymj;MpMr)o ++°+ = 


aja; AjAjaKa] 


at 
7" 24 


(IV.8) 


Matching powers of {a;} on the two sides of the above equation gives 


£ 0 for 2 odd 
ee ms = (IV.9) 
w=1 0 


sum over all pairwise contractions for @ even 


This result is known as Wick’s theorem; and for example, 


(mumpmeEMi)o = (MEM;)o(MEMI)O + (MEME) o(MsMi)o + (MsME)O(M;M1)0 - 


IV.B Expectation values in Perturbation Theory 


In the presence of an interaction U, the expectation value of any operator O, is 
computed perturbatively as follows: 
_f[DmOcHo-U fp Dme PON —U+U?/2---] 
“Tpme to ~ ST pme Mot —u+ue/2—~ 
Zell = Uo + U?)o/2=--] 


Inverting the denominator by an expansion in powers of U gives 


(0) 
(IV.10) 


=(O)9 — ((OU)o — (O)o(U)o) + : ((OU7)9 — 2(OU)o(U)o + 2(O)o (U6 — (C)o(U*)o) 


+---= se ol (OU"), . 
n=0 


n! 


(IV.11) 
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The connected averages are defined as the combination of unperturbed expectation values 
appearing at various orders in the expansion. Their significance will become apparent in 
diagrammatic representations, and from the following example. 

Let us calculate the two point correlation function of the Landau-Ginzburg model 
to first order in the parameter wu. [In view of their expected irrelevance, we shall ignore 
higher order interactions, and also only keep the lowest order Gaussian terms.] Substituting 
eq.(IV.4) into eq.([V.11) yields 

da. ada. 7d 
(ma(a) ma(a!)) = (mala) mala!)ho — u [| SET Aae ae 

{(rma(a) maa’) (a1) maa) (as)™;(—a1 — a2 — as))o- ae 

(mMa(a) ™(Q'))o(mi(a1)m:i(G2)mj;(a3)™m;(—G1 — G2 — as))o] + O(u’). 
To calculate (OU) 9 we need the unperturbed expectation value of the product of six m’s. 
This can be evaluated using eq.(I[V.9) as the sum of all pair-wise contractions, 15 in all. 
Three contractions are obtained by first pairing m, to mg, and then the remaining four 
m’s in U. Clearly these contractions cancel exactly with corresponding ones in (O)9(U)o. 
The only surviving terms involve contractions that connect O to U. This cancellation 
persists at all orders, and (OU")6 contains only terms in which all n + 1 operators are 
connected by contractions. The remaining 12 pairings in (Ol/)o fall into two classes: 
(i) 4 pairings involve contracting mq and mg to m’s with the same index, e.g. 


(™a(q) mi(q1))o(7ma(4") mi(a2))0(m; (as) m;(—a1 — a2 — gs))o 
bai 54 555 (27)3454(q + qi) 64(q! + qe) 54(qi + qa) (IV.13) 
7 (t+ Kq?) (t + Kq?) (t+ Kq3) 
where we have used eq.(IV.6). After summing over i and j, and integrating over qj, qo, 
and qg3, such terms make a contribution 
ayo Oey art a) fag 1 avis 
(t+ Kq?)? (2r)4t+ Kq 
(ii) 8 pairings involve contracting mg and mg to m’s with different indices, e.g. 


(™a(q) mi(a1))0(ma(a’) m;(q3))o(mi(a2) m;(—a1 — G2 — Gs))o 
bai 9a; 5g (27)>454(q + G1) 67(q’ + G3) 67(qi + Gs) (IV.15) 
7 (t+ Kq?) (t+ Kq”) (t+ Ka) , 
Summing over all indices, and integrating over the momenta leads to an overall contribution 


bas 27)" 6" (a +’) / d’q. 1 (v.16) 
(t+ Kq?)? (27)4t+ Kas 
Adding up both contributions, we obtain 
bag (27)764 (q+ q’) c or d%k i 
t+ Kq t+ Kq? (Q7)¢4t+ Kk? 


+ O(u?) 
(IV.17) 


(Mo(q) ma(q')) = 
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IV.C Diagrammatic Representation of Perturbation Theory 


The calculations become more involved at higher orders in perturbation theory. A 
diagrammatic representation can be introduced to help keep track of all possible contrac- 
tions. To calculate the ¢-point expectation value ee Me, (qi)), at p* order in u, proceed 


according to the following rules: 


(1) Draw @ external points labelled by (q;,a;) corresponding to the coordinates of the re- 
quired correlation function. Draw p vertices with 4 legs each, labelled by internal momenta 
and indices, e.g. {(kj, 7), (ka, 7), (ks, 7), (ka, 7)}. Since the four legs are not equivalent, the 
four point vertex is indicated by two solid branches joined by a dotted line. (The extension 


to higher order interactions is straightforward.) 


(2) Each point of the graph now corresponds to one factor of mq, (q;), and the unperturbed 
average of the product is computed by Wick’s theorem. This is implemented by joining 
all external and internal points pairwise, by lines connecting one point to another, in all 
topologically distinct ways (see #5 below). 

(3) The algebraic value of each such graph is obtained as follows: (i) A line joining a 
pair of points represents the two point average’; e.g. a connection (q,a) ——~ (q/,{), 
contributes daa (27)254(q + q’)/(t + Kq?); (ii) A vertex is represented by a term such as 


u(27)464(k, + ky + kz + ky) (the delta-function insures that momentum is conserved). 


(4) Integrate over the 4p internal momenta {k;}, and sum over the 2p internal indices. 


Note that each closed loop produces a factor of 6;; = n at this stage. 


(5) There is a numerical factor of 


eu 


p! 


x number of different pairings leading to the same topology. 


The first contribution comes from the expansion of the exponential; the second merely 


states that graphs related by symmetry give the same result, and can be calculated once. 


(6) When calculating cumulants, only fully connected diagrams (without disjoint pieces) 


need to be included. This is a tremendous simplification. 


* Because of its original formulation in quantum field theory, the line joining two points 
is usually called a propagator. In this context, the line represents the world-line of a 
particle in time, while the perturbation UY is an ‘interaction’ between particles. For the 


same reason, the Fourier index is called a ‘momentum.’ 
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IV.D Susceptibility 


It is no accident that the correction term in eq.(IV.17) is similar in form to the 
unperturbed value. This is because the form of the two point correlation function is 


constrained by symmetries, as can be seen from 


bila y= f dx x | dx! lati! (my, (x) ma(x’)). (IV.18) 


The two point correlation function in real space must satisfy translation and rotation 
symmetry (in the high temperature phase), and (ma(x) mg@(x’)) = dag(mi(x—x’) m1(0)). 


Transforming to center of mass and relative coordinates, the above integral becomes, 


(ma(q) me(q')) = 
/ 
dé (* 5 = Jat (x — x’) ellata’)-(x+%')/2gi—*')-(4-4')/25 5 my (x — x’) m1(0)) 
= (20)%0%(q + q')basS(q), 
(IV.19) 
where 
S(a) = (lma(a)l?) = f atxe!4™ (my (x — x) mx(0), (IV.20) 
is the quantity observed in scattering experiments (sec.II.D). 
From eq.(IV.17) we obtain 
1 eet? | d¢k 1 5 
S(q) = ——, {1 - —— ~~ _ | = SS HO . IV.21 
(9) t+ Kq@ t+ Kq (Qnytt+ KR (w) ( ) 
It is useful to examine the expansion of the inverse quantity 
S(q)'=t+K ?+duin+2) f SF + OW) (IV.22) 
Y= (Qn)4t+ KR 


In the high temperature phase, eq.(I[V.20) indicates that the q — 0 limit of S(q) is just 
the magnetic susceptibility x. For this reason, S(q) is sometimes denoted by x(q). From 
eq.(IV.22), the inverse susceptibility is given by 

d¢k 1 


Onttparm t Ow). eee) 


x(t) =t+4uln +2) } 


The susceptibility no longer diverges at t = 0, since 


d¢k 1 A(n+2)u Sq - A(n + 2)u Ae 
—1 = = d—3 _ 
a ae aa / (Ont KR K (Qn / ane Ko \gs9 
(IV.24) 
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is a finite number. This is because in the presence of wu the critical temperature is reduced 
to a negative value. The modified critical point is obtained by requiring y~+(t,) = 0, and 


hence from eq.(IV.23), to order of u, 


d’k 1 Au(n + 2)KqAt? 
t.=—A OY, eae ee Se ee IV.25 
cs iy (Qn)4t. + Kk (d-2)kK ~ ey 
How does the perturbed susceptibility diverge at the shifted critical point? From 
eq.(IV.23), 


Thay te) = t=  ( +2) | d"k Besa : 
- Me EN eis (Qr)¢ \t+ Kk? t,.4+ Kk? 


7 4u(n + 2) d%k 1 5 
See) Dee cam Acacia ices tesa omRae ) 
(IV.26) 


(In going from the first equation to the second, we have changed the position of t. from 
one denominator to another. Since t, = O(u), the corrections to such change only appear 
at O(u*).) The final integral has dimensions of [k7~4]. For d > 4 it is dominated by the 
largest momenta and scales as KgA¢~*. For 2 < d < 4, the integral is convergent at both 
limits. Its magnitude is therefore set by the momentum scale €~! = J(t—t.)/K , which 


can be used to make the integrand dimensionless. Hence, in these dimensions, 


2—d/2 
xt) St = t,) s- .( K ) HOG?) 


IV.2 
re foe (IV.27) 


where c is aconstant. For d < 4, the correction term at the order of u diverges at the phase 
transition, masking the unperturbed singularity of yx with y = 1. Thus the perturbation 
series is inherently inapplicable for describing the divergence of susceptibility in d < 4. 
The same situation applies in calculating any other quantity perturbatively. Although, 
we start by treating u as the perturbation parameter, it is important to realize that it is 
not dimensionless; u/K? has dimensions of (length)?~*. The perturbation series for any 
quantity then takes the form X(t,u) = Xo(t)[1 + f(uat~4/K?, ué4-4/K?)], where f is 
a power series. The two length scales a and € are available to construct dimensionless 
variables. Since € diverges close to the critical point, there is an inherent failure of the 
perturbation series. The effective (dimensionless) perturbation parameter diverges at t, 


and is not small. 
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IV.E Perturbative RG (First Order) 


The last section demonstrated how various expectation values associated with the 
Landau-—Ginzburg Hamiltonian can be calculated perturbatively in powers of u. However, 
the perturbative series is inherently divergent close to the critical point and cannot be used 
to characterize critical behavior in dimensions d < 4. Wilson showed that it is possible to 
combine perturbative and renormalization group approaches into a systematic method for 
calculating critical exponents. Accordingly, we shall extend the RG calculation of Gaussian 
model in sec.III.G to the Landau-Ginzburg Hamiltonian, by treating U = u [ d?xm* asa 


perturbation. 


1. Coarse Grain: This is the most difficult step of the RG procedure. As before, subdivide 


the fluctuations into two components as, 


m(q) for0<q<A/b 
m(q) = : (IV.28) 
a(q) for A/b<q<A 


In the partition function, 


7 A od 2 fe 
L= | pinavata) on} | on oe ) (|rm(q)|? + |o(q)|) - ute). aca 


(IV.29) 
the two sets of modes are mixed by the operator U/. Formally, the result of integrating out 


{G(q)} can be written as 


z= | Dita off" ot (ASS) con x 


exp iy if sid In (t + )| Coase = [pita e™. 


2 A/b (27)4 o 


Here we have defined the partial averages 


a A d 2 
(0), = [ oe eos re ) ton) (IV.31) 


with Z, = { Da(q) exp{—GHo|a]}, being the Gaussian partition function associated with 
the short wavelength fluctuations. From eq.(IV.30), we obtain 


GH[m] = Vo fe + ie a pape oo jin(q)|? — In (e APA) (IV.32) 
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The final expression can be calculated perturbatively as, 


_1)é 
In Cao =— W),.+5 ((u), - Uy?) Ss ( 2 x cumulant of U+---. (IV.33) 


The cumulants can be computed using the rules set in the previous sections. For example, 


at the first order we need to compute 
(u i |) = uf dard aad! asd (o7)454(q + q2 + qs + qa) 
aele. (27)44 ‘ . (IV.34) 
(rian) + #(ax)| + fre(a2) + #(a2)] [rTi(as) + e(as)| - [rma(as) + 4(a4)]) | 


The following types of terms result from expanding the product: 


fl} 1 (1(qu) - (a2) r7(qa) -m(aa)) | 

2] 4 (4(qx) - (a2) r(as) - (a4) 

[3] 2 (a1) - (a2) (as) - F(a) aan 
4) 4 (4(q) (a2) Has) -r7'(qu)) 

[5] 4 (a(a1) &(q2) F(qs3) ria) )_ 

6] 1 (an) - Ha) Fas) H(a4)), 


The second element in each line is the number of terms with the a given ‘symmetry’. 


The total of these coefficients is 2 = 16. Since the averages (Q),, involve only the short 


wavelength fluctuations, only contractions with ¢ appear. The resulting internal momenta 


are integrated from A/b to A. 
Term [1] has no & factors and evaluates to U(r]. The second and fifth terms involve 


an odd number of os and their average is zero. Term [3] has one contraction and evaluates 


to 
d’qi:--d4 55 (2m)*0%(qu + Q2) = a 
1 


A/b d"q ; A d¢k 1 
=) i et ee 
nu | (anata) i. (Qr)¢t+ Kk? 


(IV.36) 


Term [4] also has one contraction but there is no closed loop (the factor 6;;) and hence no 


factor of n. The various contractions of 4 ¢ in term [6] lead to a number of terms with 
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no dependence on m. We shall denote the sum of these terms by uVof,. Summing up all 


terms, the coarse grained Hamiltonian at order of u is given by 


Alb qdg (4 Kq? 


Sri] =v (698 + 8H) + [OG (FSS) ime 


A/b , (IV.37) 
/ d4qid@qzd“q3 aa % ae = 
ae 3d m(qi) + ™(q2)m(qs3) -m(—qi — q2 — qs) 
0 (27) 
where 
t=t+ 4u( +f ane : (IV.38) 
I ES On tt Re 


The coarse grained Hamiltonian is thus described by the same 3 parameters t, K, and uw. 


The other two parameters in the coarse grained Hamiltonian are unchanged, i.e. 


KkK=kK, and t=u. (IV.39) 


2. Rescale by setting q = b~!q’, and 


3. Renormalize, m= zm’, to get 


Ais A dé / 7 rm Kb-2 12 tial 
(amy im'| =V (692 + ws Hh) + [ote (FA) mia’? 


tustyree [Pa ae A (gh) l(a) (la — ah a) 
: (IV.40) 
The renormalized Hamiltonian is characterized by the triplet of interactions (t’, K’,u’), 
such that 


Pa ah rk a ee (IV.41) 


As in the Gaussian model there is a fixed point at t* = u* = 0, provided that we set 
Z= bits, such that kK’ = K. The recursion relations for t and wu in the vicinity of this 


point are given by 


A ad 
d°k 1 
tbe t+ au(n +2) f 


A/b (27)¢t+ Kk? 


(IV.42) 
up, = b4-4u 

While the recursion relation for wu at this order is identical to that obtained by dimensional 
analysis; the one for t is different. It is common to convert the discrete recursion relations 
to continuous differential equations by setting b = e*, such that for an infinitesimal 62, 


t, =t(b) =t(1+60)=t+ wt +O(60?) , u,=u(b)=ut et + O(60?). 
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Expanding eqs.(IV.42) to order of 52, gives 


dt 0 a 
t+ d€— = (14+ 260) (: + 4u(n + 2) agate) 
de (Qny4t+ KA (Iv 43) 


d 
ut i= =(1+(4—d)ie)u 


The differential equations governing the evolution of t and wu under rescaling are then 


d 
dt ae 4du(n + as 
i t+ KA (IV.44) 
U 
cake es 0 eee 
‘7, (4 —d)u 


The recursion relation for u is easily integrated to give u(£) = uge4—%% = upb@-®. 
The recursion relations can be linearized in the vicinity of the fixed point t* = u* = 0, 


by setting t = t* + ot and u = u* + du, as 


d ( ét g  A(n+2)KaA** ot 
may = K TV.45 
dé & (? 4—d ou ( ) 


In the differential form of the recursion relations, the eigenvalues of the matrix determine 
the relevance of operators. Since the above matrix has zero elements on one side, its 
eigenvalues are the diagonal elements, and as in the Gaussian model we can identify y, = 2, 
and y, = 4-—d. The results at this order are identical to those obtained from dimensional 
analysis on the Gaussian model. The only difference is in the eigen—directions. The 
exponent y; = 2 is still associated with u = 0, while y,, = 4 —d is actually associated 
with the direction t = —4u(n + 2)KqgA%?/K. This agrees with the shift in the transition 
temperature calculated to order of wu from the susceptibility. 

For d > 4 the Gaussian fixed point has only one unstable direction associated with y;. 
It thus correctly describes the phase transition. For d < 4 it has two relevant directions 
and is unstable. Unfortunately, the recursion relations have no other fixed point at this 
order and it appears that we have learned little from the perturbative RG. However, since 
we are dealing with an alternating series we can anticipate that the recursion relations at 


the next order are modified to 


d 
a tes ase 2 as — Av? 
U 
— = (4—d)u— Bu? 
7 (4—d)u u 
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with A and B positive. There is now an additional fixed point at u* = (4 — d)/B for 
d <4. For a systematic perturbation theory we need to keep the parameter u small. Thus 
the new fixed point can be explored systematically only for small ¢ = 4 — d; we are led to 
consider an expansion in the dimension of space in the vicinity of d = 4! For a calculation 
valid at O(e) we have to keep track of terms of second order in the recursion relation for 
u, but only to first order in that of t. It is thus unnecessary to calculate the term A in the 


above recursion relation. 
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IV.F Perturbative RG (Second Order) 


The coarse grained Hamiltonian at second order in U is 


ye (AS) m(a)P? + UM), — 5 (U2), - W)2) + OU). 


(IV.47) 


Hm] = Vofe 
To calculate (Ww), — uy?) we need to consider all possible decompositions of two Us 
into m and @ as in eq.(IV.34). Since each YU can be broken up into 6 types of terms as 
in eq.(IV.35), there are 36 such possibilities for two Us which can be arranged in a 6 x 6 
table. Many of the elements of this table are either zero, or can be neglected at this stage, 


due to a number of considerations: 


(i) All the 11 terms involving at least one factor of type [1] are zero because they cannot 
be contracted into a connected piece, and the disconnected elements cancel in calculating 


the cumulant. 


(ii) An additional 12 terms (such as [2] x [3]) involve an odd number of os and are zero 


due to their parity. 


(iii) Two terms, [2] x [5] and [5] x [2], involve a vertex where two Gs are contracted together, 
leaving a 77(q<) and a &(q>). This configuration is not allowed by the d-function which 


ensures momentum conservation for the vertex, as by construction q7 + q< # 0. 


(iv) Terms [3] x [6], [4] x [6] and their partners by exchange have two factors of 7m. They 
involve two loop integrations, and appear as corrections to the coefficient t. We shall denote 
their net effect by A, which as noted earlier does not need to be known precisely at this 


order. 


(v) The term [5] x [5] also involves two factors of 771, while [2] x [2] includes 6 such factors. 
The latter is important as it indicates that the space of parameters is not closed at this 
order. Even if initially zero, a term proportional to m® is generated under RG. In fact, 
considerations of momentum conservation indicate that both these terms are zero for q = 0, 


and are thus contributions to g2m? and q?m®. We shall comment on their effect later on. 


(vi) The contributions resulting from [6] x [6] are constants, and will be collectively denoted 


by u2V6f?. 
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(vii) The terms [3] x [3], [3] x [4], [4] x [8], and [4] x [4] contribute to fi. For example, 


[3] x [3] results in 
A dk, d4kyd%k), dk}, 


o xaxaxe fo Pon dan f 
0 (27)*4 Jaye (27)*4 
x (27)2454(qi + qo + ki + k2)d4(ky + ke + q3 + qa) 
Ovet( 20 q5¢ (ey 4k! ed Og oa( le + k. an ny S, a, 
( ) ( 1 1) ( ) ( = 2) Fai) ; m(qz2)i7i(q3) : m(qa) — 
t+ Kk5 


t+ Kk? 
A/b dé a dé ~ es is = 
anu? [ agar 2m) 5an + q2 + q3 + qa)m(qi) + 7(q2)m(qs) - 7(qa) 
0 
d"k 1 
(IV.48) 


x SS ee eeeeeEeEeEeEEEEE—EEEEEEESSSSSS==—=_—_ 
(27)? (t+ Kk?) (t + K(qi + G2 — k)*) 
The contractions from terms [3] x [4], [4] x [3], and [4] x [4] lead to similar expressions with 


prefactors of 8, 8, and 16 respectively. Apart from the dependence on q, and qp, the final 


result has the form of U/ [77]. In fact the last integral can be expanded as 
dk i 2Kk - (qi + q2) — K(qi + qo)? 
_— — eels 
Man + 4) : (nyt G+ KR (f+ KR) " 
After fourier transforming back to real space we find in addition to m*, such terms as 


(IV.49) 


m?(Vm)?,m?V?m?, ++: 
Putting all contributions together, the coarse grained Hamiltonian at order of u? takes 


the form 
GH=V (ofp + udfp + u7dfZ) 


K@ ae as 


Qnjtt+ Kee 2 


A/b d?q 
i Pe TP 
Lexa) us 
vie d"qi---d"qa = = = z 
[Sr Bila) i(aa)ri(as) F(a) 
u? OS dk af 
x c — —(8n + 64) I, (on)? G+ KR? + O(u?q?)| + O(u?m®q?,---) + O(u?) 
(IV.50) 


IV.G_ The <Expansion 
The parameter space (Kt, u) is no longer closed at this order; several new interactions 


proportional to m?, m*, and m®, all consistent with symmetries of the problem, appear in 
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the coarse grained Hamiltonian at second order in wu. Ignoring these interactions for the 


time being, the coarse grained parameters are given by 
K = K —v?A"(0) 


Faure dGias) ie any. 24 2 A(0) 
= n u —— ——__— — 4 
ayy Om)? t+ KE (IV.51) 


A oad 
k 1 
aau-4in+8ye? f ce 
A/b (20 )* (t + Kk?) 


where A(0) and A”(0) correspond to the first two terms in the expansion of A(t, K,q?) in 
eq.(IV.50) in powers of q. 
After the rescaling q = b~'q’, and renormalization m = zm’, steps of the RG proce- 


dure, we obtain 
K' =)p-*-? 27K, fhe Be tee, w= bo eet. (IV.52) 


As before, the renormalization parameter z is chosen such that AK’ = K, leading to 


Dts be = pet2 O 2 
gris Towa = (1+ O(u*)). (IV.53) 


The value of z does depend on the fixed point position u*. But as u* is of the order of e€, 
Z= git g+O(e?), is not changed at the lowest order. Using this value of z, and following 


the previous steps for constructing differential recursion relations, we obtain 


dt 4u(n+ 2)KqA? 


—=2 — A(t, K, A)u? 
Oe gee (t, K, Aju 
: (IV.54) 
Ae cer _ A(n+ 8)KaA 2 
de ote KR)? 


The fixed points are obtained from dt/dé = du/dé = 0. In addition to the Gaussian 
fixed point at u* = t* = 0, discussed in the previous section, there is now a non-trivial 


fixed point located at 


(et + KA?) K?2 ; 
= a = eee O 
ee intem ae) 
2u*(n + 2)KaA4 (n + 2) 9 9 , 
: 4+ KR ee 


The above expressions have been further simplified by systematically keeping terms to first 


order ine = 4—d. 


66 


Linearizing the recursion relations in the vicinity of the fixed point results in 


(IV.56) 


n+2) K, e A(n+2)KaA4 sg 

d (ot\ _ 2— re u* — Alu? aa ae — 2Au ot 

dé \ ou) — 8(n+8)KaA* , «2 8(n+8)Kah*. « . 
(+R AZ)s © “erHKAzy 


At the Gaussian fixed point, t* = u* = 0, eq.(IV.45) is reproduced. At the new fixed point 
of eqs.(IV.55), 


A(n+2)K4A* Ke 
as 6 ee 1G) pei 
= 8(n+8)KaA* Ke ; : 
eexee O(e*) «ee Tn+8)Ka ne 
We have not explicitly calculated the top element of the second column as it is not necessary 
for calculating the eigenvalues. This is because the lower element of the first column is 


zero to order of ¢«. Hence the eigenvalues are determined by the diagonal elements alone. 


The first eigenvalue is positive, controlling the instability of the fixed point, 


(n + 2) 2 
=2-— O(e*). IV.58 
Yt rey (€*) ( ) 
The second eigenvalue, 
Yu = —€ + O(€?), (IV.59) 


is negative for d < 4. The new fixed point thus has co-dimension of one and can describe the 
phase transition in these dimensions. It is quite satisfying that while various intermediate 
results, such as the position of the fixed point depend on such microscopic parameters as 
K and A, the final eigenvalues are pure numbers, only depending on n and d = 4~—e. 
These eigenvalues characterize the universality classes of rotational symmetry breaking in 
d < 4. (Long range interaction may lead to new universality classes.) 


The divergence of the correlation length, € ~ (6t)~”, is controlled by the exponent 


yodafafi-Oe2 Jy a2 tet? soe) (IV.60) 


The singular part of the free energy scales as, f ~ (dt)?~%, and the heat capacity diverges 


with the exponent 


(4 —€) ln Bae 4—n . 
= — —— —_— 1 — — = ° I 61 
a=2-—dv=2 5 ia ig" ae (IV.61) 


To complete the calculation of critical exponents, we need the eigenvalue associated with 


the (relevant) symmetry breaking field h. This is easily found by adding a term hie 
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fd?xm(x) = —h + m™(q = 0) to the Hamiltonian. This term is not effected by coarse 
graining or rescaling, and after the renormalization step changes to —zh- m/(q’ = 0), 


implying 
W=czh=b0+th, => m=1+S+O(2) =3-£4+0(), (IV.62) 


The vanishing of magnetization as T’ — T, is controlled by the exponent 


— d-Yyn _ fe 1 n+2 5 
b= 7 =( 5 1) x5 (14+ e+ (2) 


IV.63 
See e+ O(e?) | | 
2 2n+8) 
while the susceptibility diverges as y ~ (dt)~7, with 
2yn — d 1 n+ 2 Nn+2 9 
= 7 =2~x ~ (14+ — —~e] =14+ — ~e4+ O(e’). IV.64 
- Yt 2 ( 2(n + 8) ) 2(n + 8) (e) ( ) 


Using the above results, we can estimate various exponents as a function of d and n. 
For example, for n = 1, by setting « = 1 or 2 in eqs.(IV.60) and (IV.63) we obtain the 
values v(1) © 0.58, v(2) © 0.67, and G(1) ¥ 0.33, (2) + 0.17. The best estimates of these 
exponents in d = 3 are v & 0.63, and @ = 0.32. In d = 2 the exact values are known to be 
vy =1and @ = 0.125. The estimates for 9 are quite good, while those for v are less reliable. 
It is important to note that in all cases these estimates are an improvement over the mean 
field (saddle point) values. Since the expansion is around four dimensions, the results are 
more reliable in d = 3 than in d = 2. In any case, they correctly describe the decrease of 
G with lowering dimension, and the increase of v. They also correctly describe the trends 


with varying n at a fixed d as indicated by the following table of exponents a(n). 


ae = 2 = ta 
O(e) ate=1 0.17 0.11 0.06 0 
Experiments in d = 3 0.11 -0.01 -0.12 — 


Although the sign of a is incorrectly predicted at this order for n = 2 and 3, the decrease 


of a with increasing n is correctly described. 
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IV.H_ Irrelevance of other interactions 


The fixed point Hamiltonian at O(e) (from eqs.(IV.55)) has only three terms 


gine d 9 (n+2) ,2 5 A“ Ky 
ae 2 A Scere IV.65 
BH 5 [atx |ovm 8) © m + On 8) reuk : ( ) 


and explicitly depends on the imposed cutoff A ~ 1/a (unlike the exponents). However, as 
described in section III.E, the starting point for RG must be the most general Hamiltonian 
consistent with symmetries. We also discovered that even if some of these terms are left 
out of the original Hamiltonian, they are generated under coarse graining. At second order 
in u, terms proportional to m® were generated; higher powers of m will appear at higher 
orders in u. 

Let us focus on a rotatationally symmetric Hamiltonian for h =0. We can incorporate 
all terms consistent with this symmetry in a perturbative RG by setting GH = GH) +U, 
where 

K 


UL 
BHo = i dx 5m? + £eom)?+ £072m)? +--+), (IV.66) 


includes all quadratic (Gaussian terms), while the remaining higher order terms are placed 


in the perturbation 
w= jax [um* + um?(Vm)? +-+-+ ugm® +---+ugm® +--+]. (IV .67) 
After coarse graining, and steps (ii) and (iii) of RG in real space, x = bx’ and m = Cm’, 
the renormalized weight depends on the parameters 
try b9C7t = bt 
Ke bt?’CK=K 
De bP CL ab AL 


uv b9¢4a = 04% 


IV.68 
vi bt-2¢4% = b27-45 ( ) 


Ug be b°C us =: ame VA 


Ug te bP tig = ba 
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The second set of equalities are obtained by choosing ¢? = b?-¢4K /K = §7-4/1 + 


O(u?, uv, v?,---)], such that K’ = K. By choosing an infinitesimal rescaling, the recursion 


relations take the differential forms 


dt 

dt 
dk 
dt 

dL 
dt 


=2f + O(u, v, U6, Us, 2 -) 
=0 


= —2L4 O(u?, uv, v?,---) 


— =eu— Bu? + O(ur,v”,---) 
: ; (IV.69) 
ee (-2+6€)u+ O(u?, uv, v,- +) 


dl 
dug 8.252 
“de (—2 + 2e)ug + Olu”, ug, - --) 
d 


These recursion relations describe two fixed points: 


(i) The Gaussian fixed point, t* = L* = u* = v* =---=0, and K 40, has eigenvalues 


y? = 2, yf? — =7Z, weeny y? = hE) y? = =2F6, ee ap — —2+2¢, ys = —4+3¢, ae ae (IV.70) 
(ii) Setting eqs.(IV.69) to zero, a non-trivial fixed point is located at 
thru ~ Ole), Le vut ~~ Oe), uve s+ Ole), os. (IV.71) 


The stability of this fixed point is determined by the matrix, 


bt 2—-O(u*) Ole) +++ OI) OW) --- bt 
5L Ole?)  -2 + Oe) 5L 

d | 2 me 

dé| dul O(e*) O(e) du |- 2) 
du O(e?) O 


(€) a 
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Note that as « — 0, the non-trivial fixed part, its eigenvalues and eigendirections 
continuously go over to the Gaussian fixed points. Hence the eigenvalues can only be 


corrected by order of €, and eq.(IV.70) is modified to 


m+2 2 
+ O(2), SPP Olan: 
a i a te) (IV.73) 


Yu = —€+ O(€7), yo = —2+ O(C), +, Y6 = —2+ O(€), yg = —44+ Of), Sas 


V=a2= 


While the eigenvalues are still labelled with the coefficients of the various terms in the 
Landau—Ginzburg expansion, we must remember that the actual eigendirections are now 
rotated away from the axes of this parameter space, although their largest projection is 
still parallel to the corresponding axis. 

Whereas the Gaussian fixed point has two relevant directions in d < 4, the generalized 
O(n) fixed point has only one relevant direction corresponding to y,. At least perturba- 
tively, this fixed point has a basin of attraction of co-dimension one, and thus describes the 
phase transition. The original concept of Kadanoff is thus explicitly realized and the uni- 
versality of exponents is traced to the irrelevance (at least perturbatively) of the multitude 
of other possible interactions. The perturbative approach does not exclude the existence 
of other fixed points at finite values of these parameter. The uniqueness of the critical 
exponents observed so far for each universality class, and their proximity to the values cal- 
culated from the e-expansion, suggests that postulating such non-perturbative fixed points 


is unnecessary. 


IV.I Comments on the --Expansion 


The perturbative implementation of RG for the Landau—Ginzburg Hamiltonian was 
achieved by K.G. Wilson in the early 1970s; the e-expansion was developed jointly with 
M.E. Fisher. This led to a flurry of activity in the topic which still continues. Wilson 
was awarded the Nobel prize in 1982. Historical details can be found in his Nobel lecture 
reprinted in Rev. Mod. Phys. 55, 583 (1983). A few comments on the e-expansion are in 


order at this stage. 


1. Higher orders, and convergence of the series: Calculating the exponents to O(c)? and 
beyond, by going to order of U/° and higher, is quite complicated as we have to keep track 
of many more interactions. It is in fact quite unappealing that the intermediate steps 
of the RG explicitly keep track of the cutoff scale A, while the final exponents must be 


independent of it. In fact there are a number of field theoretical RG schemes (dimensional 
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regularization, summing leading divergences, etc.) that avoid many of these difficulties. 
These methods are harder to visualize and will not be described here. All higher order 
calculations are currently performed using one of these schemes. It is sometimes (but not 
always) possible to prove that these approaches are consistent with each other, and can 
be carried out to all orders. In principle, the problem of evaluating critical exponents in 
d = 3 is now solved: simple computations lead to approximate results, while more refined 
calculations should provide better answers. The situation is somewhat like finding the 
energy levels of a He atom which can not be done exactly, but which may be obtained with 
sufficient accuracy using various approximation methods. 

To estimate how reliable the exponents are, we need some information on the conver- 
gence of the series. The € expansion has been carried out to the fifth order, and the results 
for the exponent y, for n = 1 at d= 3, are 

y = 140.167 + 0.077e? — 0.049e? + 0.180e* — 0.415e° 

: (IV.74) 

1.2385 + 0.0025 = 1.000, 1.167, 1.244, 1.195, 1.375, 0.96 

The second line compares the values obtained at different orders by substituting « = 1, 
with the best estimate of y ~ 1.2385 in d = 3. Note that the elements of the series 
have alternating signs. The truncated series evaluated at € = 1 improves up to third 
order, beyond which it starts to oscillate, and deviates from the left hand side. These are 
characteristics of an asymptotic series. It can be proved that for large p, the coefficients 
in the expansion of most quantities scale as |fp| ~ cp!a~”. As a result, the «expansion 
series is non-convergent, but can be evaluated by the Borel summation method, using the 


identity, tie dx x?e—* = pl, as 


1. ff? aes p 
fle) = Da Sot? = x era f dz z?e~* = / dre” sy — (IV.75) 


The final summation (which is convergent) results in a function of x which can be integrated 
to give f(e). Very good estimates of exponents in d = 3, such as the one for y quoted 
above, are obtained by this summation method. There is no indication of any singularity 


in the exponents up to € = 2, corresponding to the lower critical dimension d = 2. 
2. The 1/n expansion: The fixed point position, 


_ (* + KA?)?(4-d) 
—  A(n+8)KqgAt ’ 


* 


ihe 


vanishes as n — oo. This suggests that a controlled 1/n expansion of the critical exponents 
is also possible. Indeed such an expansion can be developed by a number of methods, such 
as a saddle point expansion that takes advantage of the exponential dependence of the 
Hamiltonian on n, or by an exact resummation of the perturbation series. Eq.(I[V.58) in 
this limit gives, 
2 1 
ye Mink |e ig as =. (IV.76) 
n— oo n- — 
This result is exact in dimensions 4 < d < 2. Above four dimensions the mean field value 


of 1/2 is recovered, while for d < 2 there is no order. 
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V. Position Space Renormalization Group 


V.A_ Lattice Models 


While Wilson’s perturbative RG provides a systematic approach to probing critical 
properties, carrying out the «expansion to high orders is quite cumbersome. Models de- 
fined on a discrete lattice provide a number of alternative computational routes that can 
complement the perturbative RG approach. Because of universality, we expect that all 
models with appropriate microscopic symmetries and range of interactions, no matter how 
simplified, lead to the same critical exponents. Lattice models are convenient for visualiza- 
tion, computer simulation, and series expansion purposes. We shall thus describe models 
in which an appropriate “spin” degree of freedom is placed on each site of a lattice, and 
the spins are subject to simple interaction energies. While such models are formulated in 
terms of explicit ‘microscopic’ degrees of freedom, depending on their degree of complexity, 
they may or may not provide a more accurate description of a specific material than the 
Landau-Ginzburg model. The point is that universality dictates that both descriptions 
describe the same macroscopic physics, and the choice of continuum or discrete models is 
a matter of computational convenience. Some commonly used lattice models are described 
here: 

1. The Ising Model is the simplest and most widely applied paradigm in statistical me- 
chanics. At each site i of a lattice, there is a spin 0; which takes the two values of +1 or 
—1. Each state may correspond to one of two species in a binary mixture, or to empty 
and occupied cells in a lattice approximation to an interacting gas. The simplest possible 


short range interaction involves only neighboring spins, and is described by a Hamiltonian 
H= \" B (Oyo), (V.1) 
<1,j> 


where the notation < 2,7 > is commonly used to indicate the sum over all nearest neighbor 
pairs on the lattice. Since c? = 1, the most general interaction between two spins is 


A 


B(o,o') = -—g- “(0 +o’) — Jao’. (V.2) 


For N spins, there are 2% possible micro-states, and the (Gibbs) partition function is 


Lh yertes Y= exp K S° Baath Dee ; (V.3) 


{oi} {oi} <1,j> 
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where we have set K = BJ,h = Bh, and g = Gg (6 =1/keT, and z is the number of bonds 
per site, i.e. the coordination number of the lattice). For h = 0 at T = 0, the ground 
state has a two fold degeneracy with all spins pointing up or down (K > 0). This order 
is destroyed at a critical K, = J/kgT,. with a phase transition to a disordered state. The 
field h breaks the up-down symmetry and removes the phase transition. The parameter g 
merely shifts the origin of energy, and has no effect on the relative weights of microstates, 
or the macroscopic properties. 


All the following models can be regarded as generalizations of the Ising model. 


2. The O(n) model: Each lattice site is now occupied by an n-component unit vector, i.e 
S; = (Sj, S?,---,S"), with ye ($2)? = (V.4) 


A nearest-neighbor interaction can be written as 


H=-IS\~ §- Se De (V.5) 


<t,j> 


In fact, the most general interaction consistent with spherical symmetry is f (S; . S;) for 
an arbitrary function f. Similarly, the rotational symmetry can be broken by a number of 
“fields” such as yh . §,)P. Specific cases are the Ising model (n = 1), the XY model 
(n = 2), and the Heisenberg model (n = 3). 


3. The Potts model: Each site of the lattice is occupied by a q-valued spin S; = 1,2,---,¢. 


The interactions between the spins are described by the Hamiltonian 


H=-—J S- Os; 8; FD bs. (V.6) 


<i,j> 


The field h now breaks the permutation symmetry amongst the q-states. The Ising model 
is recovered for gq = 2, since 6¢,¢/ = (1+00')/2. The 3 state Potts model can for example 
describe the distortion of a cube along one of its faces. Potts models with gq > 2 represent 
new universality classes not covered by the O(n) model. Actually, the transitions for g > 4 


ind=2, and q>3 ind =3 are discontinuous. 


4. Spin s-models: The spin at each site takes the 2s + 1 values, s; = —s,—-s+1,---,+s. 


A general nearest-neighbor Hamiltonian is 


H = ys (Jisis; + Jo(sis;)? +--+ Jo5(s; i8j)° oe Sj. (V.7) 


<i,j> 
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The Ising model corresponds to s = 1/2, while s = 1 is known as the Blume-Emery-Griffith 
(BEG) model. It describes a mixture of non-magnetic (s = 0), and magnetic (s = +1) 
elements. This model exhibits a tricritical point separating continuous and discontinuous 
transitions. However, since the ordered phase breaks an up-down symmetry, the phase 
transition belongs to the Ising universality class for all values of s. 


Some of the computational tools employed in the study of discrete models are: 


1. Position space renormalizations: These are implementations of Kadanoff’s RG scheme 
on lattice models. Some approximation is usually necessary to keep the space of inter- 
actions tractable. Most such approximations are uncontrolled; a number of them will be 


discussed here. 


2. Series expansions: Low-temperature expansions start with the ordered ground state and 
examine the lowest energy excitations around it. High temperatures expansions begin with 
the collection of non-interacting spins at infinite temperature and include the interactions 
between spins perturbatively. Critical behavior is then extracted from the singularities of 


such series. 


3. Exact solutions can be obtained for a very limited subset of lattice models. These 
include one dimensional chains which will be solved by the transfer matrix method in 


recitations, and the two dimensional Ising model as described later in the course. 


4. Monte-Carlo simulations: The aim of such methods is to generate configurations of spins 
that are distributed with the correct Boltzmann weight exp(—37). There are a number of 
methods, most notably the Metropolis algorithm, for achieving this aim. Various expecta- 
tion values and correlation functions are then directly computed from these configurations. 
With the continuing increase of computer power, numerical simulations have become in- 
creasingly popular. Limitations of the method are due to the size of systems that can be 
studied, and the amount of time needed to ensure that the correctly weighted configura- 
tions are generated. There is an extensive literature on numerical simulations which will 


not be further discussed in this course. 
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V.B Exact treatment in d=1 


An exact RG treatment can be carried out for the Ising model with nearest neighbor 
interactions (eq.(V.1)) in one dimension. The basic idea is to find a transformation that 
reduces the number of degrees of freedom by a factor b, while preserving the partition 
function, i.e. 

Ze So ee Ck eee (V.8) 
{o;|t=1,---,N} {o!, |’ =1,---,N/b} 


There are many mappings {o;} > {o/,}, that satisfy this condition. The choice of the 
transformation is therefore guided by the simplicity of the resulting RG. With b = 2, 
for example, one possible choice is to group pairs of neighboring spins and define the 
renormalized spin as their average. This majority rule, 0) = (a2i-1 + 02;)/2, is in fact not 
very convenient as the new spin has three possible values (0,+1) while the original spins 
are two valued. We can remove the ambiguity by assigning one of the two spins, e.g. o2;_-1, 
the role of tie-breaker whenever the sum is zero. In this case the transformation is simply 
o; = 09:1. Such an RG procedure effectively removes the even numbered spins, s; = 2; 
and is usually called a decimation. Note that since o’ = +1 as in the original model, no 
renormalization factor ¢, is necessary in this case. Since the interaction is over adjacent 


neighbors, the partition function can be written as 


N 
S © B(o;, 0141) 


w=1 


N 
Z= Y= exp 


{oi} 


~ >» do exp S [B(o%, 81) + B(si, %441)| : (V.9) 


| N/2 N/2 N/2 
{oped [it 


Summing over the decimated spins, {s;}, leads to 


N/2 
e- BH lo!) = I] | x POCA) = ons Bo ti41) (V.10) 


i=1 eed. 


where following eq.(V.2), 


h 
B(o1,02)=g+t+ 9 (01 + 92) + Koren, (V.11) 
and 
h! 
B’(o1,05) =9' + ace +05)+ K'o}05, (V.12) 


are the most general interaction forms for Ising spins. 
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Following eq.(V.10), the renormalized interactions are obtained from 


/ 


h 
(oi, 05) = exp [Kook + Flo! +08) +9 


i (V.13) 
= S- exp [Kiet +04) + 5 (7% +04) + hs; + 2g]. 
8y=+1 
To solve for the renormalized interactions it is convenient to set 
g=e*, y=e, Z=e7 
, , , (V.14) 
g=e* , y =e", 2 =e 
The four possible configurations of the bond are 
R(+, +) =a’y’2! = 27y(2?y + 2-7y™) 
R(-, —) =z'y'—1 2! = 2*y\ (xy dic x?y—*) 
(V.15) 
R(+,—) =a! Az! = 27 (yty™) 
R(-,+) =a! 2! = A(yty™) 


The last two equations are identical, resulting in three equations in three unknowns, with 


the solutions, 
24 = (ay ta ?y (a 2y + ?y iy ty YP 


12 __ go y oii tage a 


y ay + ay? : (V.16) 
it ety + ety) (ey + wy") 
way t)? 


Taking the logarithms, we find recursion relations of the form 
g =2g + dg(K, h) 
h' =h+ 6h(K,h) ; (VAT) 
K' =K'(K,h) 

The parameter g is just an additive constant to the Hamiltonian. It does not effect the 


probabilities and hence does not appear in the recursion relations for K and h. dg(K,h) 
is the contribution of the decimated spins to the overall free energy. 

1. Fixed points: The h = 0 subspace is closed by symmetry, and it can be checked that 
for y = 1 eqs.(V.16) imply y’ = 1, and 


aie? (— 2 e 2k 


2 
i 
er = 5 ) ‘ = i= 3 incosh 2K. (V.18) 
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The recursion relation for K has the following fixed points: 

(a) An infinite temperature fixed point at K* = 0, which is the sink for the disordered 
phase. If K is small, K’ ~ In(1 + 4K?/2)/2 ~ K?, is even smaller, indicating that 
this is a stable fixed point with zero correlation length. 

(b) A zero temperature fixed point at K* — oo, describing the ordered phase. For a 
large but finite K, the renormalized interaction K’ ~ In (e?“/2) /2 ~ K —1n2/2, 
is somewhat smaller. This fixed point is thus unstable with an infinite correlation 
length. 

Clearly any finite interaction renormalizes to zero, indicating that the one dimensional 

chain is always disordered at sufficiently long length scales. The absence of any other 

fixed point is apparent by noting that the recursion relation eq.(V.18) can alternatively be 

written as tanh K’ = (tanh K’)?. 


2. Flow diagrams indicate that in the presence of a field h, all flows terminate on a line 
of fixed points with K* = 0 and arbitrary h*. These fixed points describe independent 
spins and all have zero correlation length. The flows originate from the fixed point at 


k* = h* =0 which has two unstable directions in the (K,h) parameter space. 


3. Linearizing the recursion relations around this fixed point (x — co) yields, 


alt mw at /A Fs PY Pa 
, —« 


y? = y* ed 


(V.19) 


We can thus regard e~* and h as scaling fields. Since €’ = €/2, the correlation length in 


the vicinity of the fixed point satisfies the homogeneous form (b = 2) 


é(e* ,h) = 2€ (v2e-*, 2n) 


— 2¢ (ees 2‘) eae 


The second equation is obtained by repeating the RG procedure ¢ times. Choosing @ such 


that 2°/2e-* = 1, we obtain the scaling form 
ia (e*, h) = eK ge (he?*) : (V.21) 


The correlation length diverges on approaching T = 0 for h = 0. However, its divergence 
is not a power law of temperature. There is thus an ambiguity in identifying the exponent 


vy related to the choice of the measure of vicinity to T = 0 (1/K or e~*). 
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The hyperscaling assumption states that the singular part of the free energy in d 


dimensions is proportional to €~¢. Hence we expect 
Arr C.@ec) wes é-1= e °K gy (he?*) : (V.22) 


At zero field, the magnetization is always zero, while the susceptibility behaves as 


Of 2K 


x(k) ~ ah fe ~ (V.23) 


On approaching zero temperature, the divergence of the susceptibility is proportional to 
that of the correlation length. Using the general forms, (s;,5j42) ~ e~*/§/a4-?*", and 
x~f d?*x < 898% >c~ E7-", we conclude that 7 = 1. 

The one dimensional Ising model can be solved exactly by a transfer matrix method 
which will be described in more detail in the recitations. The results of the direct solution 
are in full agreement with those of the RG. The transfer matrix approach can be applied to 
any one dimensional chain with nearest neighbor interactions, since the partition function 


can be written as 


Z= Y= exp 
{si} 


Defining a transfer matriz, 


N N 
> B(si, =) = ys I] ereneary), (V.24) 


i=l {s;}i=1 


(iT ]3j) = eF@+99), (V.25) 


and with periodic boundary conditions, we obtain 
Zit" |e Me (V.26) 


Note that for N — ov, the trace is dominated by the largest eigenvalue Amax. Frobenius’s 
theorem states that for any finite matrix with finite positive elements, the largest eigenvalue 
is always non-degenerate. This implies that Amax and Z are analytic functions of the 
parameters appearing in B, and that one dimensional models can exhibit singularities 
(and hence a phase transition) only at zero temperature (when some matrix elements 
become infinite). The transfer matrix method also provides an alternative RG scheme for 
all such one dimensional chains. For decimation by a factor b, we can use Z = tr (r ye) 


to construct the rescaled bond energy as 


eB(si083) = (si |T"|s5) = (si|T?|s). (V.27) 
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V.C The Niemeijer—van Leeuwen Cumulant Approximation 


Unfortunately, the decimation procedure cannot be performed exactly in higher di- 
mensions. For example, the square lattice can be divided into two sublattices. For an 
RG with b = V2, we can start by decimating the spins on one sublattice. The interac- 
tions between the four spins surrounding each decimated spin are obtained by generalizing 


eq.(V.13). If initially h = g = 0, we obtain 


RG) 65;03;04) = Ss" eKs(aiteztostes) — 2cosh[K(o, +0,+04+04)].  (V.28) 


s=+1 


Clearly the four spins appear symmetrically in the above expression, and hence are subject 
to the same two body interaction. This implies that new interactions along the diagonals 
of the renormalized lattice are also generated, and the nearest neighbor form of the original 


Hamiltonian is not preserved. There is also a four point interaction, and 
R=exp(g' + K'(o}05 + 0405 + 0304 + 040) + 0105 + 0504) + Kyojogo304]. (V.29) 


The number (and range) of new interactions increases with each RG step, and some trun- 
cating approximation is necessary. Two such schemes are described in the following sec- 
tions. 

One of the earliest approaches was developed by Niemeijer and van Leeuwen (NvL) 
for treating the Ising model on a triangular lattice, subject to the usual nearest neighbor 
Hamiltonian —GH = kK ae j) F195: The original lattice sites are grouped into cells of three 
spins (e.g. in alternating up pointing triangles). Labelling the three spins in cell a as 


{o., 02,02}, we can use a majority rule to define the renormalized cell spin as 
ot, = sign [oi +02 +03]. (V.30) 


(There is no ambiguity in the rule for any odd number of sites, and the renormalized spin is 
two-valued.) The renormalized interactions corresponding to the above map are obtained 


from the constrained sum 


e~PH led] — » eSHlel (V.31) 


To truncate the number of interactions in the renormalized Hamiltonian, NvL intro- 


duced a perturbative scheme by setting GH = GH o +U. The unperturbed Hamiltonian 
—BHo = ke) (0402 + 0303 + 030%), (V.32) 


involves only intra—cell interactions. Since the cells are decoupled, this part of the Hamil- 
tonian can be treated exactly. The remaining inter—cell interactions are treated as a 


perturbation 


uU=K Se (of) 0) + ofa). (V.33) 
<a,B> 


The sum is over all neighboring cells, each connected by two bonds. (The actual spins 
involved depend on the relative orientations of the cells.) Eq.(V.31) is now evaluated 


perturbatively as 
6H'(o4] BHoloi) ue 
/ Pal = ~PEFolad | —-—-::-]. V.34 
e 2 | - U + 5 | ( ) 


The renormalized Hamiltonian is given by the cumulant series 
1 
BH'[o4] = In Zoloe] + Mo - 5 (U7) - CG) + OW), (V.35) 


where (), refers to the expectation values with respect to GHo, with the restriction of fixed 
[o/,], and Zp is the corresponding partition function. 
To proceed, we construct a table of all possible configurations of spins within a cell, 


their renormalized value, and contribution to the cell energy: 


+/+ 
+)+]+]P° 


+ 
| 
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The restricted partition function is the product of contributions from the independent cells, 


/ 


ZoloiJ=[] | So ek aratearatoara)| = (eK + 3e-K we (V.36) 
cl Oa cee 
It is independent of |[o/,], thus contributing an additive constant to the Hamiltonian. The 


first cumulant of the interaction is 
o=K SO [(ob) 9 (02 )5 + (08) (08)o] =2K D2 (od)o (03) (V.37) 
<a,B> <a,B> 


where we have taken advantage of the equivalence of the three spins in each cell. Using 


the table, we can evaluate the restricted average of site spins as 


fi =+1 = 
lot) _ aaa ee or Og =+ 7 Bie s i ee) 
Pan =e". -ec* 396% ; ~ @8K 4 3e-K ~e" ; 
Bk gee for o,=- 
Substituting in eq.(V.37) leads to 
N Sys = 3K 4 9K 2 

/ = 

—BH'lon] = 3 In(e + 3e )+2K (SA =) a a3 +O(U 2). (V.39) 


At this order, the renormalized Hamiltonian involves only nearest neighbor interactions, 


with the recursion relation 


eK 4 ¢-K \? 
K'S9K |: |}, A 
(= + —) ey) 


1. Eq.(V.40) has the following fixed points: 

(a) The high temperature sink at K* = 0. If K «1, K' = 2K (2/4)? = K/2 < K, i. 
this fixed point is stable, and has zero correlation length. 

(b) The low temperature sink at K* = oo. If K > 1, then K'’ + 2K > K, i.e. unlike the 
one dimensional case, this fixed point is also stable with zero correlation length. 

(c) Since both of the above fixed points are unstable, there must be at least one stable 
fixed point at finite kK’ = kK = K*. From eq.(V.40), the fixed point position satisfies 

i pak” ag" 


Bee STE ee, ASE 201k" 2 = AK" 43. V.41 
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The fixed point value 


V2-1 


can be compared to the exactly known value of 0.2747 for the triangular lattice. 


1 (8-4/2 
K*=<In (4) ~ 0.3356, (V.42) 


2. Linearizing the recursion relation around the non-trivial fixed point yields, 


OK" 


ae Lay? » (eAK* 41) 
pai = ~ 1.624. V.43 
az |. 72 ( (v.43) 


ae 39 Kk* AK 2 
ae ap c (tk + 3)3 


The fixed point is indeed unstable as required by the continuity of flows. This RG scheme 
removes 1/3 of the degrees of freedom, and corresponds to b = 3. The thermal eigenvalue 
is thus obtained as 


OK! 
OK |x 


_ In(1.624) 
on) 


This can be compared to the exactly known value of y; = 1, for the two dimensional Ising 


bvt = 


~ 0.883. (V.44) 


model. It is certainly better than the mean-field (Gaussian) estimate of y, = 2. From this 


eigenvalue we can estimate the exponents 
y=1/y,21.13 (1), and a=2-—2/y,=—-0.26 (0), 


where the exact values are given in the brackets. 


3. To complete the calculation of exponents, we need the magnetic eigenvalue yp, obtained 


after adding a magnetic field to the Hamiltonian, i.e. from 
BH =BHo +U-hd_ ot, . (V.45) 


Since the fixed point occurs for h* = 0, the added term can also be treated perturbatively, 


and to the lowest order 


BH' = BHo + (U oh A (a +02 +03)),, (V.46) 


where the spins are grouped according to their cells. Using eq.(V.38), 


ee NL 
BH’ =InZ)+K' So o 0-30 (Setar z) (V.47) 


3K 
mass ers + 3e7 
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thus identifying the renormalized magnetic field as 


ert sue 
hiS=3h |. A 
(oe + ir] Ou) 


In the vicinity of the unstable fixed point 


Oh! ek" 44 3 
bye _—_- —_ = ee a V.49 
Ol | es Ak" 43 V2 ( ) 
and Ja 
In (3/2 
BNI 2) ae (V.50) 


This is lower than the exact value of y, = 1.875. (The Gaussian value of y;, = 2 is closer 


to the correct result in this case.) 


4. NvL carried out the approach to the second order in YU. At this order two additional 
interactions over further neighbor spins are generated. The recursion relations in this three 
parameter space have a non-trivial fixed point with one unstable direction. The resulting 
eigenvalue of y; = 1.053, is tantalizingly close to the exact value of 1, but this agreement 


is probably accidental. 


V.D The Migdal—Kadanoff Bond Moving Approximation 


Consider a b = 2 RG for the Ising model on a square lattice, in which every other spin 
along each lattice direction is decimated. As noted earlier, such decimation generates new 
interactions between the remaining spins. One way of overcoming this difficulty is to simply 
remove the bonds not connected to the retained spins. The renormalized spins are then 
connected to their nearest neighbors by two successive bonds. Clearly after decimation, the 
renormalized bond is given by the recursion relation in eq.(V.18), characteristic of a one 
dimensional chain. The approximation of simply removing the unwanted bonds weakens 
the system to the extent that it behaves one dimensionally. This is remedied by using the 
unwanted bonds to strengthen those that are left behind. The spins that are retained are 


now connected by a pair of double bonds (of strength 2A’), and the decimation leads to 
it 
i= 5 Incosh(2 x 2k). (V.51) 


1. Fixed points of this recursion relation are located at 
(a) K* =0: For K <1, K's In(1+8K?)/2=4k? « K, ie. this fixed point is stable. 
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(b) K* — oo: For K > 1, K’ & In(e**/2)/2 = 2K > K, indicating that the low 
temperature sink is also stable. 
(c) The domains of attractions of the above sinks are separated by a third fixed point at 
etk* 4 94K" 


ek" — — ,——, => K*~ 0.305, (V.52) 


which can be compared with the exact value of K, & 0.441. 
2. Linearizing eq.(V.51) near the fixed point gives 


_ OK" 


Ut 
e OK 


= 2tanh4k* = 1.6786, => y% 0.747, (V.53) 
K* 
compared to the exact value of y = 1. 
The bond moving procedure can be extended to higher dimensions. For a hypercubic 
lattice in d-dimensions, the bond moving step strengthens each bond by a factor of 2¢—!. 


After decimation, the recursion relation is 
/ 1 d—1 
KS 5 Incosh [2 x 2°-'K]. (V.54) 
The high and low temperature sinks at K* = 0 and K* — ov, are stable, since 


ta ae ay Sea Sd eas cm Gece 2 (V.55) 


1 
K <i, => KS 5 


and 
1 elk 
KS 1, => K's 5a 27'k > K. (V.56) 
(Note that the above result correctly identifies the lower critical dimension of the Ising 
model, in that the low temperature sink is stable only for d > 1.) The intervening fixed 
point has an eigenvalue 
_ OK! 


ey =? Stank (27K), (V.57) 
K* 


The resulting values of K* ~ 0.065 and y% & 0.934 for d = 3, can be compared with the 
known values of K, & 0.222 and y, © 1.59 on a cubic lattice. Clearly the approximation 
gets worse at higher dimensions. (It fails to identify an upper critical dimension, and as 


d— oo, K* = 220-4 and y — 1.) 
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The Migdal—Kadanoff scheme can also be applied to more general spin systems. For 
a one dimensional model described by the set of interactions {A}, the transfer matrix 


method in eq.(V.27) gives the recursion relations as 
Ty({K"}) = T({K})”. 


For a d-dimensional lattice, the bond moving step strengthens each bond by a factor of 


b¢-1, and the generalized Migdal—Kadanoff recursion relations are 
Ty({K"}) = T({bo* K})’. (V.58) 


The above equations can be used as a quick way of estimating phase diagrams and 
exponents. The procedure is exact in d = 1, and does progressively worse in higher 
dimensions. It thus compliments mean-field (saddle point) approaches that are more 
reliable in higher dimensions. Unfortunately, it is not possible to develop a systematic 
scheme to improve upon its results. The RG procedure also allows evaluation of free 
energies, heat capacities, and other thermodynamic functions. One possible worry is that 
the approximations used to construct RG schemes may result in unphysical behavior, e.g. 
negative values of response functions C' and x. In fact most of these recursion relations 
(e.g. eq.(V.58)) are exact on hierarchical (Berker) lattices. The realizability of such lattices 


ensures that there are no unphysical consequences of the recursion relations. 
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VI. Series Expansions 


VI.A Low-temperature expansions 


Lattice models can also be studied by series expansions. Such expansions start with 
certain exactly solvable limits, and typically represent perturbations around such limits 
by graphs on the lattice. High temperature expansions are described in the next section. 
Here we describe the low temperature expansion for the Ising model with a Hamiltonian 
—-BH = K oe j) %ij, ON a d-dimensional hypercubic lattice. The ground state with 
K = BJ > 0 is ferromagnetic, e.g. with o; = +1 for all spins. A series expansion for the 
partition function is obtained by including low energy excitations around this state. The 
lowest energy excitation is a single overturned spin. Any of N sites can be chosen for this 
excitation, which has an energy cost of 2k x 2d, with respect to the ground state. The 
next lowest energy excitation corresponds to a dimer of negative spins with energy cost of 
2K x (4d—2), and a multiplicity of N x d (there are d possible orientations for the dimer). 
The first few terms in the expansion give 


N(N —2d-1 
( So BaRe ie 5 os 


iio QeNdk 1+ Ne4¢4 4+ dNe t@4-K 4 
2 


(VI.1) 


The fourth term comes from flipping two disjoint single spins. The zeroth order term 
comes from the ground state, which has a two fold degeneracy. The overall factor of 2 is 


insignificant in the N — oo limit, and 


Zw eNdk : e 2kx boundary of droplet (V1.2) 


droplets of © spins 


The free energy per site is obtained from the series, 


InZ 1 N(N — 2d—1 
—Gf = BELA I 4+ —In]1+ Ne 744% 4 dNe4@4-HK 4 N(N ~ 2d — 1) |8ax A ee 
N N 2 
— dK + 0-49 4 de~42d-DK _ Catt) eeu at 
(V1.3) 


Note that there is an explicit cancellation of the term proportional to N? in the expansion 


for Z. Such higher N dependances result from disconnected diagrams with several disjoint 
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droplets. Extensivity of the free energy ensures that these terms are cancelled by products 


of connected graphs. The energy per site is then obtained from 


Bie 0 (\InZ aa 0 (\InZ 
N08 (4) =“, x(a) (V4) 


=—J id — 4de~*"* — 4d(2d — 1)e4C*-Y¥ 4 dd (2d + Le 8* +--+], 


and the heat capacity is proportional to 
CG. TOE Ie Or 
Nkg NkgOT NJOK (VL5) 
—K? [16-44 + 16d(2d — 1)2e~4@4-D¥ _ 394? (2d + 1)e784* 4... 


Can such series be used to extract the critical coupling K,, and more importantly, the 
singularities associated with the disordering transition? Suppose that we have identified 
a number of terms in a series C = )77°,, agu’. From the expected divergence of the heat 


capacity, we expect an asymptotic expansion of the form 


CR et ae Sak ep a og Oe) 
Ue tha Qu? flué 
(V1.6) 


The above singular form is characterized by the three parameters A, u,., and a. We can 
try to extract these parameters from the calculated coefficients in the series by requiring 


that they match at large @, i.e. by fitting the ratio of successive terms to 


ae a+fl-1 wg a—1 
— + | —— ] = 1 : ; 
a ( Om ) Us ( + ? ) (V1.7) 
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Thus a plot of ag/ag_; versus 1/@ should be a straight line with intercept uz *, and slope 
uz'(a—1). Note, however, that adding a finite sum paras deu* does not change the 
asymptotic singular form, but essentially renders the first @,, terms useless in determining 
a and uz!. There is thus no a priori guarantee that such a fitting procedure will succeed 
with a finite number of coefficients. In practice, this procedure works very well and very 
good estimates of critical exponents in d = 3 (e.g. a = 0.105 + 0.007) are obtained in this 
way by including a rather large number of terms. 

The three terms calculated for the heat capacity in eq.(VI.5) have different signs, 
unlike those of eq.(VI.6). As this continues at higher orders in e~*, the ratio fitting 
procedure described above cannot be used directly. The alternation of signs usually signifies 


singularities in the complez z = e~* plane, closer to the origin than the critical point of 
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interest at real z. = ec. If it is possible to construct a mapping u(z) on the complex 
plane such that the spurious singularities are pushed further away than u, = u(z.), then 
the ratio method can be used. In the case of low temperature series, the choice of u = 
tanh K achieves this goal. (As we shall demonstrate shortly, tanh K is also the natural 
variable for the high temperature expansion.) Quite sophisticated methods, such as Padé 
approximants, have been developed for analyzing the singular behavior of series. 

Low temperature expansions can be similarly constructed for other discrete spin sys- 
tems, such as the Potts model. For continuous spins, the low energy excitations are Gold- 
stone modes, and the perturbation series cannot be represented in terms of lattice graphs. 
The low temperature description in this case starts with the Gaussian treatment of the 
Goldstone modes. Further terms in the series involve interactions between such modes, 


and a corresponding calculation will be performed later on. 


VI.B High-temperature Expansions 


High temperature expansions work equally well for discrete and continuous spin sys- 
tems. The basic idea is to start with independent spins, and expand the partition function 


in powers of 3 = (kgT)~?, i.e. 


2 
Z=tr(e at) ta [a ou 2% salted ( (VL8) 
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The averages (), are calculated over non-interacting spins. For the Ising model, it is more 
convenient to organize the expansion in powers of tanh K as follows. Since (oj;0;)? = 1, 


the Boltzmann factor for each bond can be written as 


ek te * K -K 


eX rs — i aaa + = ge cosh K (1+ to;o;), (VI.10) 


where t = tanh K is a good high temperature expansion parameter. Applying this trans- 


formation to each bond of the lattice results in 


Z = So oN Luan = (cosh Kym Ff bonds S~ TT (1 + toics). (VL.11) 
{oi} {oi} (6.3) 
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For N, bonds on the lattice, the above product generates 2% terms, which can be 
represented diagrammatically by drawing a line connecting sites i and 7 for each factor 
of to;a;. Note that there can at most be one such line for each lattice bond, which is 
either empty or occupied. This is a major simplification, resulting from the use of t, rather 


Pi 


than K, as the expansion parameter. Each site now obtains a factor of o;', 


where p; is 
the number of occupied bonds emanating from 7. Summing over the two possible values 
o; = +1, gives a factor of 2 if p; is even, and 0 is p; is odd. Thus, the only graphs that 
survive the sum have an even number of lines passing through each site. The resulting 


graphs are collections of closed paths on the lattice, and 


Are QN x (cosh Kk)Ne , t number of bonds in the graph (VI.12) 
All closed graphs 


For a d-dimensional hypercubic lattice, the smallest closed graph is a square of 4 bonds 


which has d(d — 1)/2 possible orientations. As the next graph has 6 bonds, 


d(d —1)N 


5 did =1)0d =3)e +1 5 (V1.13) 


Z = 2" (cosh K)*” c + 


and 
InZ d(d—1 
M2 —m2+dincosnK + {Dey (VI.14) 
In the following sections, we shall employ high temperature expansions not as a nu- 
merical tool, but to establish the following: (a) Exact solution of the Ising model in d = 1. 


(b) The duality relating models at low and high temperatures. (c) The validity of the 


Gaussian model in high dimensions. (d) Exact solution of the Ising model in d = 2. 


VI.C_ Exact Solution of the One Dimensional Ising Model 


The graphical method provides a rapid way of solving the Ising model at zero field 
in d = 1. We can compare and contrast the solutions on chains with open and closed 


(periodic) boundary conditions. 


1. An open chain of N sites has N, = N —1 bonds. It is impossible to draw any closed 


graphs on such a lattice, and hence 


ESO cost KY Se). Ss She In[2 cosh | — (VI.15) 


The same method can also be used to calculate the correlation function (amon), since 


1 O10; hk)N 
(Oy) — Z a Daeg a "OmOn = a S- OmOn II ( (1 -- toio;). 
{oi} {oi} (2,3) 
(VI1.16) 


The terms in the numerator involve an additional factor of omo,n. To get a finite value 
after summing over o,, = +1 and o, = +1, we have to examine graphs with an odd 
number of bonds emanating from these external sites. The only such graph for the open 


chain is one that directly connects these two sites, and results in 


1 


mon — zlm—n| — o—|m—nl/é ith Soh ee) 
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(VI1.17) 


These results are in agreement with the RG conclusions of sec.V.B; the correlation length 
diverges as e?* for K — oo, and there is no power law modifying the exponential decay 


of correlations. 


2. A closed chain has the same number of sites and bonds, N. It is now possible to draw 


a closed graph that circles the whole chain, and 


Z = (2cosh K)* [i+ t%] =o" (cosh KN +sinh kK), 


InZ In [1 + al (VI.18) 
NV = In(2 cosh kK’) + as | i 


The difference between the free energies of closed and open chains vanishes in the thermo- 
dynamic limit. The correction to the extensive free energy, N In(2 cosh Kk), is of the order 
of 1/N for the open chain, and can be regarded as the surface free energy. There is no such 
correction for the closed chain, which has no boundaries; there is instead an exponential 
term t’. The correlation function can again be calculated from eq.(VI.16). There are two 


paths connecting the points m and n, and 


tlm—n| A tN-|m—n| 


Ll 
144% (Erg) 


(OmOn) = 


Note that the final answer is symmetric with respect to the two ways of measuring the 


distance between the two sites m and n. 
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VI.D_ Self—Duality in the Two Dimensional Ising Model 


Kramers and Wannier discovered a hidden symmetry that relates the properties of the 
Ising model on the square lattice at low and high temperatures. One way of obtaining this 
symmetry is to compare the high and low temperature series expansions of the problem. 


The low temperature expansion has the form 


Vs [1+ Ne*%?* + 2Ne~ OPK 4] 
—p2NK S- e 2K xperimeter of island ; (VI.20) 
Islands of (—) spins 


The high temperature series is 


Z =2% cosh K?% [1 + Ntanh K* +2N tanh K® +.- ‘| 
—9N cosh K2N SS tanh K lensth of graph ; (VI.21) 


graphs with 2 or 4 lines per site 
As the boundary of any island of spins serves as an acceptable graph (and vice versa), 
there is a one to one correspondence between the two series. Defining a function g to 


indicate the logarithm of the above series, the free energy is given by 


ne =2K+g kere =1n2+2IncoshK + g(tanhk). (VI.22) 


The arguments of g in the above equation are related by the duality condition 


e?’* .stanhK, => K = D(K)= -5 Intanh K . (VI.23) 
The function g (which contains the singular part of the free energy) must have a special 
symmetry that relates its values for dual arguments. (For example the function f(x) = 
x/(1+ 27) equals f(a~+), establishing a duality between the arguments at x and x71.) 
Eq.(VI.23) has the following properties: 

1. Low temperatures are mapped to high temperatures, and vice versa. 


2. The mapping connects pairs of points since D(D(k)) = K. This condition is established 


by using trigonometric identities to show that 


sinh2K =2sinh K cosh K = 2tanh K cosh? K 
2tanh k 2e7-2K 2 1 (VI.24) 


~1—tanh?K 1—e-4K eK —e-2K  sinh2K 


94 


Hence, the dual interactions are symmetrically related by 


sinh2K -sinh2k =1 . (VI.25) 


3. If the function g(K) is singular at a point K,, it must also be singular at K,. Since 
the free energy is expected to be analytic everywhere except at the transition, the critical 


model must be self dual. At the self dual point 


qe e2kKe 


oat Ke — — 
e 2 Se iy oT 


which leads to the quadratic equation, 
e *Ke 4 9e-Ke _1 = 0, — e *Ke — 14 /2 
Only the positive solution is acceptable, and 


K.=—5mn(v2~-1) = 5m(v2+1) =0.4a-., (VI.26) 


4. As will be explored in the next section and problem sets, it is possible to obtain dual 
partners for many other spin systems, such as the Potts model, the XY model, etc. While 
such mappings place useful constraints on the shape of phase boundaries, they generally 
provide no information on critical exponents. (The self—duality of many two dimensional 


models does restrict the ratio of critical amplitudes to unity.) 


VI.E Dual of the Three Dimensional Ising model 


We can attempt to follow the same procedure to search for the dual of the Ising model 


on a simple cubic lattice. The low temperature series is 


7 —p3NK [1 obs Ne 2 x6 ae 3Ne72h x10 Jag 7 
—e3NK S- e 2K xarea of island’s boundary (V1.27) 


Islands of (—) spins 
By contrast, the high temperature series takes the form 


Z =2% cosh K°% [1+3N tanh K* + 18N tanh K® +--+] 
—9N cosh K3N S- tanh number of lines ; (VI.28) 


graphs with 2, 4, or 6 lines per site 
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Clearly the two sums are different, and the 3d Ising model is not dual to itself. Is there any 
other Hamiltonian whose high temperature series reproduces the low temperature terms 
for the 3d Ising model? To construct such a model note that: 


~2kxarea and the area is the sum of faces 


1. The low temperature terms have the form e 
covered on the cubic lattice. Thus plaquettes must replace bonds as the building blocks of 


the required high temperature series. 

2. How should these plaquettes be joined together? The number of bonds next to a site 
can be anywhere from 3 to 12 while there are only two or four faces adjacent to each bond. 
This suggests ‘glueing’ the plaquettes together by placing spins o; = +1, on the bonds of 
the lattice. 

3. By analogy to the Ising model, we can construct a partition function, 


z= 3» [] (it tanh Kaba2aba$) oc Sock Le MP8 (VT.29) 


{a4,=+1}  plaquettes P {a} 


where c+, are used to denote the 4 dual spins around each plaquette, and K is given by 
eq.(VI.23). This partition function describes a system dual to the original 3d Ising model, 
in the sense of reproducing its low temperature series. (Some reflection demonstrates 
that the low temperature expansion of the above partition function reproduces the high 
temperature expansion of the Ising model.) 

Eq.(VI.29) describes a Z2 lattice gauge theory. The general rules for constructing 
such theories in all dimensions are: (i) Place Ising spins ¢; = +1, on the bonds of the 
lattice. (ii) The Hamiltonian is —3H = K ¥7,n piaquettes || @p- In addition to the global 
Ising symmetry, 6; — —¢;, the Hamiltonian has a local (gauge) symmetry. To observe this 
symmetry, select any site and change the signs of spins on all bonds emanating from it. As 
in each of the faces adjacent to the chosen site two bond spins change sign, their product, 
and hence the overall energy, is not changed. 

There is a rigorous proof (Elitzur’s theorem) that there can be no spontaneous sym- 
metry breaking for Hamiltonians with a local symmetry. The essence of the proof is that 
even in the presence of a symmetry breaking field h, the energy cost of flipping a spin is 
finite (6h for the gauge theory on the cubic lattice). Hence the expectation value of spin 
changes continuously as h — 0. (By contrast, the energy cost of a spin flip in the Ising 
model grows as Nh.) This theorem presents us with the following paradox: Since the three 


dimensional Ising model undergoes a phase transition, there must be a singularity in its 


96 


partition function, and also that of its dual. How can there be a singularity in the parti- 
tion function of the three dimensional gauge theory if it does not undergo a spontaneous 
symmetry breaking? 

To resolve this contradiction, Wegner suggested the possibility of a phase transition 
without a local order parameter. The two phases are then distinguished by the asymptotic 
behavior of correlation functions. The appropriate correlation function must be invariant 
under the local gauge transformation. For example, the Wilson loop is constructed by 


selecting a closed path of bonds S, on the lattice and examining 


Cg = (Product of ¢ around the loop) = (I a) : (VI1.30) 

ieS 
As any gauge transformation changes the signs of two bonds on the loop, their product is 
unaffected and Cg is gauge invariant. Since the Hamiltonian encourages spins of the same 
sign, this expectation value is always positive. Let us examine the asymptotic dependence 
of C's on the shape of the loop at high and low temperatures. In a high temperature 
expansion, the correlation function is obtained as a sum of all graphs constructed from 


plaquettes with S as a boundary. Each plaquette contributes a factor of tanh K, and 
Coal [Ja ek he thonebet 


= (tanh ro ore E + O (tanh K*)| & exp I-F (tanh A) x Area of s| : 

(VI.31) 

The low temperature expansion starts with the lowest energy configuration. There are 

in fact Ng = 2% such ground states related by gauge transformations. The Np plaquette 
interactions are satisfied in the ground states, and excitations involve creating unsatisfied 
plaquettes. Since C's is gauge independent, it is sufficient to look at one of the ground 
states, e.g. the one with o; = +1 for all 7. Flipping the sign of any of the 3N bonds 
creates an excitation of energy 8K with respect to the ground state. Denoting the number 


of bonds on the perimeter of the Wilson loop by Ps, we obtain 


0, Ne ee E +(3N — Ps)e~?**4(+1) + (-1)Pge2%*4 + --- 
SS OC oO IRS Oe 
Ne cK Ne [1+ 3Ne2h x4 4..,| (VL.32) 
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The asymptotic decay of C's is thus different at high and low temperatures. At high 
temperatures, the decay is controlled by the area of the loop, while at low temperatures it 
depends on its length. The phase transition marks the change from one type of decay to 
the other, and by duality has the same singularities in free energy as the Ising model. 

The prototype of gauge theories in physics is quantum electrodynamics (QED), with 
the action 


S= eax G (-i9 +eX+m)~+ ; (OA Oy Ay"). (VI1.33) 


The spinor ~ is the Dirac field for the electron, and the 4—-vector A describes the electro- 
magnetic potential. The phase of ~ is not observable, and the action is invariant under 
the gauge symmetry 7 + e*°?y, and A, + A, — 0,0, for any o(x). The Z2 lattice gauge 
theory can be regarded as the Ising analog of QED with the bond spins playing the role 
of the electromagnetic field. We can introduce a ‘matter’ field by placing spins s; = +1, 


on the sites of the lattice. The two fields are coupled by the Hamiltonian, 


BH =JS— 8,638; + KY apapeRas, (VI.34) 
(1,9) P 
where o;; is the spin on the bond joining i and 7. The Hamiltonian has the gauge symmetry 
8, +> (—1)s;, and 6;,,, —~ —6;,,, for all bonds emanating from any site 7. 

Regarding one of the lattice directions as time, a Wilson loop is obtained by creating 
two particles at a distance x, propagating them for a time t, and then removing them. The 
probability of such an event is roughly given by C, ~ e~/@)*, where U(z) is the interaction 
between the two particles. In the high temperature phase, Cg decays with the area of the 
loop, suggesting U(x)t = f(tanh K)|z|t. The resulting potential U(x) = f(tanh K)|z|, is 
like a string that connects the particles together. This is also the potential that describes 
the confinement of quarks at large distances in quantum chromodynamics. The decay with 
the length of the loop at low temperatures implies U(x)t = gle) (ar +t). For t > |a| 
the potential is a constant and the force vanishes. (This asymptotic freedom describes the 
behavior of quarks at short distances.) The phase transition implies a change in the nature 


of interactions between particles mediated by the gauge field. 
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VI.F Summing over Phantom Loops 


The high temperature series can be approximately summed so as to reproduce the 
Gaussian model. This correspondence provides a better understanding of why Gaussian 
behavior is applicable in high dimensions, and also prepares the way for the exact summa- 
tion of the series in two dimensions in the next section. The high temperature series for 
the partition function of the Ising model on a d—dimensional hypercubic lattice is obtained 
from 

Z= Ss" po ae, 795 — ON cosh*™ K x S, (VI.35) 

{oi} 
where S is the sum over all allowed graphs on the lattice, each weighted by t = tanh K 
raised to the power of the number of bonds in the graph. The allowed graphs have an 
even number of bonds per site. The simplest graphs have the topology of a single closed 
loop. There are also graphs composed of disconnected closed loops. Keeping in mind the 


cumulant expansion, we can set 
= = sum over contribution of all graphs with one loop, (V1.36) 


and introduce another sum, 


1 1 
S’ = exp (=) =1+5+,() +2) + 


=1+ (1 loop graphs) + (2 loop graphs) + (3 loop graphs) + ---. 


(VI.37) 


Despite their similarities, the sums S and S’ are not identical: There are ambiguities 
associated with loops that intersect at a single site, which will be discussed more fully in 
the next section. More importantly, S’ includes additional graphs where a particular bond 
contributes more than once, while in the original sum S$, each lattice bond contributes a 
factor of 1 or t. This arises because after raising = to a power ¢, a particular bond may 
contribute up to times for a factor of t’. In the spirit of the approximation that includes 
multiple appearances of a bond, we shall allow additional closed paths in ©, in which a 
particular bond is traversed more than once in completing the loop. Qualitatively, S is the 
partition function of a gas of self-avoiding polymer loops with a monomer fugacity of t. 
The self—avoiding constraint is left out in the partition function S’, which thus corresponds 


to a gas of phantom polymer loops, which may pass through each other with impunity. 


99 


Loops of various shapes can be constructed from closed random walks on the lattice, and 


the corresponding free energy of phantom loops is 


lS S- all closed random walks on the lattice x ¢lensth of walk 


té (V1.38) 
=N S- z (number of closed walks of @ steps starting and ending at 0). 
£ 


Note that extensivity is guaranteed since (up to boundary effects) the same loop can be 
started from any point on the lattice. The overall factor of 1/¢ accounts for the @ possible 
starting points for a loop of length @. 

A transfer matriz method can be used to count all possible (phantom) random walks 


on the lattice. Let us introduce a set of N x N matrices, 
(i]W (€)|j) = number of walks from j to i in @ steps, (VI1.39) 


in terms of which eq.(VI.38) becomes 


ae 5 7 (OW (0). (V1.40) 
e 


The additional factor of 2 arises since the same loop can be traversed by two random walks 


moving in opposite directions. Similarly, the spin-spin correlation function 


(o@)o(r)) =F Y= o(O)o(r) |] 1 + toia;), (VI.41) 

{oi} (ij) 
is related to the sum over all paths connecting the points 0 and r on the lattice. In 
addition to the simple paths that directly connect the two points, there are disconnected 
graphs that contain additional closed loops. In the same approximation of allowing all 
intersections between paths, the partition function S’ can be factored out of the numerator 


and denominator of eq.(V1.41), and 
(a(O)o(r) = S° t! (r|W(é)|0). (VI.42) 
e 


The counting of phantom paths on a lattice is easily accomplished by taking advantage 
of their Markovian property. This is the property that each step of a random walk proceeds 
from its last location and is independent of its previous steps. Hence, the number of walks 


can be calculated recursively. First, note that the any walk from O to r in @ steps can be 
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accomplished as a walk from 0 to some other point r’ in @— 1 steps, followed by a single 


step from r’ to r. Summing over all possible locations of the intermediate point leads to 


(r|W(€)|0) = 50 (r|W(1)[r’) x (r’|W(e — 1)|0) 
ry (VI1.43) 
= (r|TW(é— 1)|0), 


where the sum corresponds to the product of two matrices, and we have defined T = W(1). 


The recursion process can be continued and 
W()=TW(E-1)=T?W(e- 2)? =--- = 7. (VI.44) 


Thus all lattice random walks are generated by the transfer matrix T, whose elements are 


(rT Ir’) = { 1 ifr and r’ are nearest neighbors (VL45) 
O otherwise 


(It is also called the adjacency, or connectivity matrix.) For example in d = 2, 
(x, y|T|2',y") = Oy,y! (Oe,a/-+1 iy bx,2'—1) + dz, (Oy,y/41 ae by,y’—1) ) (VI.46) 


and successive actions of T’ on a walker starting at the origin |x, y >= dz,0dy,0, generate 


the patterns 


OO 20 0 
00 0 01 0 02 0 2 0 
Goose FS ae Oa FES a sh oe 
Oo ON So? he BA OP 20 ah, 

0010 0 


The value at each site is the number of walks ending at that point after @ steps. 

Various properties of random walks can be deduced from diagonalizing the matrix 
T. Due to the translational symmetry of the lattice, this is achieved in the Fourier basis 
(r\q) = e’47 /VN. For example in d = 2, starting from eq(VI.46), it can be checked that 


(a, ylT de, dy) = 5 (a, ylT Ia", y') (2, y'ldes Gy) 


tal 
vy 


== jet Cae 4 etde aD) 4 Jee (ctw (vr) 4 emu) 
N 


i og 
ae [2 cos gz + 2cos gy] = T'(de, dy) (2, yldx, Wy) - 
(VL.47) 
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The generalized eigenvalue for a d—dimensional hypercubic lattice is 


Ea) = 2 S- COS Ja - (VI1.48) 


The correlation function in eq.(VI.42) is now evaluated as 


~ te = Yee tT)*|0) 


= 4 ). = ‘ Cl) ay (al (VL49) 


= | d?q ear 1 =) dé eat 
(27)? N 1—2t574_, cos qa (27)4 1—2tS°, cosda 


For t — 0, the shortest path costs least energy and (a(0)o(r)) ~ t!"!. As ¢ increases, larger 
paths dominate the sum because they are more numerous (i.e. entropically favored). 
Eventually, there is a singularity for 1 — t7(0) = 0, ie. at 2d x t. = 1, when arbitrarily 
long paths become important. For t < t,, the partition function is dominated by small 
loops, and a polymer connecting two far away points is stretched by its line tension. When 
the fugacity exceeds t,, the line tension vanishes and loops of arbitrary size are generated. 
Clearly the neglect of intersections (which stabilizes the system at a finite density) is no 
longer justified in this limit. This transition is the manifestation of Ising ordering in the 
language of paths representing the high temperature series. On approaching the transition 
from the high temperature side, the sums are dominated by very long paths. Accordingly, 


the denominator of eq.(VI.49) can be expanded for small q as 


d 
1—tT(q) =1-2t cos qa © (1 — 2dt) + tq + O(q*)  te(E? +? + O(G')), (VL50) 
a=1 
where ij 
1 — 2dt\~ 
£= ; ; ) : (VI.51) 
The resulting correlation functions, (a (0 x f Aa, ei r/(q? + €~7), are identical to 


those obtained from the Gaussian ete i 


a2 for << (n = 0) 
(o(O)o(r)) x) oye (VI.52) 
7(d=ty/2 for 7 > & 
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The correlation length in eq.(VI.51) diverges as € ~ (te — t)71/? 


, le. with the Gaussian 
exponent of vy = 1/2. 


We can also calculate the free energy in eq.(VI.40) as 


= iS 1 el hs 
gop OO jo) = (0 a 0) 
d 
= -} (olin - er)I0) = -% ak oti (O|q)In(1—#7(q)) (ql0)——(V.53) 


1 d4q 
=-5 | ou In (: nay 


In the vicinity of the critical point at t. = 1/(2d), the argument of the logarithm is 
proportional to (q7+£€~7) from eq.(VI.50). This is precisely as in the Gaussian model, and 


as discussed earlier the singular part of the free energy scales as 
joing O86 © Bg tn (VI.54) 


The singular part of the heat capacity, obtained after taking two derivatives, is governed 
by the exponent a = 2 —d/2. Note that in evaluating the sums appearing in eqs.(VI.49) 
and (VI.53), the lower limit for @ is not treated very carefully. The series in eq.(VI.49) is 
assumed to start from @ = 0, and that of eq.(VI.53) from @ = 1. In fact the first few terms 
of both series may be zero because the number of steps is not sufficient to reach 0 from r, 
or to from a closed loop. This is not a serious omission, in that the singular behavior of a 
series is not effected by its first few terms. Treating the first few terms properly can only 
add analytic corrections to the singular forms calculated in eqs.(VI.49) and (VI.53). 

The equivalence of these results to the Gaussian model is a manifestation of field— 
particle duality. In a field theoretical description, (imaginary) time appears as an addi- 
tional dimension, and the two point correlations describe the probability of propagating a 
particle from one point in space-time to another. In a wave description, this probability 
is calculated by evolving the wave function using the Schrédinger equation. Alternatively, 
the probability can be calculated as the sum over all (Feynman) paths connecting the two 
points, each path weighted with the correct action. The second sum is similar to the above 
calculation of (o(r)o(0)). 

This approach provides an interesting geometrical interpretation of the phase tran- 
sition. The establishment of long range order implies that all parts of the system have 


selected the same state. This information is carried by the bonds connecting nearest 
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neighbors, and can be passed from the origin to a point r, through all paths connecting 
these two points. The fugacity t is a measure of the reliability of information transfer 
between neighboring sites. Along a one dimensional chain, unless t = 1, the transferred 
information decays at large distances and it is impossible to establish long range order. 
In higher dimensions there are many more paths, and by accumulating the information 
from all paths it is possible to establish order at t. < 1. Since the number of paths of 
length ¢ grows as (2d) while their information content decays as t*, the transition occurs 
at t. = 1/(2d). (A better approximation is obtained by including some of the constraints 
by noting that the random walk cannot back track. In this case the number of walks grows 
as (2d —1)*.) The total information from paths of length @ is weighted by (2dt)’, and 
decays exponentially for t < t-. The characteristic path length, 2 = —1/In(2dt), diverges 
as (t- — t)~! on approaching the transition. For paths of size  < @ there is very good 
information transfer. Such paths execute random walks on the lattice and cover a distance 
gmt 1/2 
walk nature of the paths. 


. The divergence of v with an exponent of 1/2 is thus a consequence of the random 


Why does the classical picture fail for d < 4? Let us focus on the dominant paths close 
to the phase transition. Is it justified to ignore the intersections of such paths? Random 
walks can be regarded as geometrical entities of fractal (Hausdorf) dimension dy = 2. This 
follows from the general definition of dimension relating the mass and extent of an object 
by M x R4@f, and the observation that the size of a random walk (R « €) is the square 
root of its length (M « @). Two geometrical entities of dimensions d; and d2 will generally 
intersect in d-dimensional space if d; + dz > d. Thus our random walkers are unlikely to 
intersect ind > d, = 2+2 = 4, and the above (Gaussian) results obtained by neglecting the 
intersections are asymptotically valid. Below the upper critical dimension of 4, random 
walks have frequent encounters and their intersections must be treated correctly. The 
diagrams obtained in the perturbative calculation of the propagator with um* correspond 
precisely to taking into account the intersections of paths. (Each factor of u corresponds 
to one intersection.) It is now clear that the constraint of self-avoidance will swell the 
paths beyond their random walk size leading to an increase in the exponent v. Below 
the transition the length of paths grows without bound and the self—avoiding constraint is 
necessary to ensure the stability of the system. 

The loop expansion is easily generalized to n-component spins. The only difference 
is that each closed loop now contributes a factor of n. In the phantom limit, where 


intersections are ignored, the free energy (eq.(VI.53)) is simply multiplied by n, while the 
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correlation function is left unchanged (precisely as in the Gaussian model). Corrections due 
to intersections, which modify critical behavior in d < 4, now depend on n. For example, 
in correcting the two point correlation function, we have to subtract contributions from 
the self-intersection of the random walk, as well as from contacts with loops (which have 
a fugacity of n). The correspondence with the perturbative series of the propagator with 


2 


a nonlinearity u(m -m)°, is again apparent. 
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VI.G_ Exact free energy of the Square Lattice Ising model 


As indicated in eq.(VI.35), the Ising partition function is related to a sum S, over 
collections of paths on the lattice. The allowed graphs for a square lattice have 2 or 4 
bonds per site. Each bond can appear only once in each graph, contributing a factor of 
t = tanh kK. While it is tempting to replace S with the exactly calculable sum S$’, of all 
phantom loops of random walks on the lattice, this leads to an overestimation of S. The 
differences between the two sums arise from intersections of random walks, and can be 
divided into two categories: 

(a) There is an over-counting of graphs which intersect at a site, i.e. with 4 bonds through 
a point. Consider a graph composed of two loops meeting at a site. Since a walker 
entering the intersection has three choices, this graph can be represented by three 
distinct random walks. One choice leads to two disconnected loops; the other two are 
single loops with or without a self-crossing in the walker’s path. 

(b) The independent random walkers in S$’ may go through a particular lattice bond more 
than once. 

Including these constraints amounts to introducing interactions between paths. The 
resulting interacting random walkers are non—Markovian, as each step is no longer inde- 
pendent of previous ones and of other walkers. While such interacting walks are not in 
general amenable to exact treatment, in two dimensions an interesting topological property 


allows us to make the following assertion: 


ie Se collections of loops of random walks with no U turns (VL55) 


x poumber of bonds x ol anes of crossings 


The negative signs for same terms reduce the overestimate and render the exact sum. 


Proof: We shall deal in turn with the two problems mentioned above. 


(a) Consider a graph with many intersections and focus on a particular one. A walker 
must enter and leave such an intersection twice. This can be done in three ways only one 
of which involves the path of the walker crossing itself (when the walker proceeds straight 
through the intersection). This configuration carries an additional factor of (-1) according 
to eq.(VI.55). Thus, independent of other crossings, these three configurations sum up to 
contribute a factor of 1. By repeating this reasoning at each intersection, we see that the 
over-counting problem is removed and the sum over all possible ways of tracing the graph 


leads to the correct factor of one. 
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(b) Consider a bond that is crossed by two walkers (or twice by the same walker). We can 
imagine the bond as an avenue with two sides. For each configuration in which the two 
paths enter and leave on the same side of the avenue, there is another one in which the 
paths go to the opposite side. The latter involves a crossing of paths and hence carries a 
minus sign with respect to the former. The two possibilities thus cancel out! The reasoning 
can be generalized to multiple passes through any bond. The only exception is when the 
doubled bond is created as a result of a U-turn. This is why such backward steps are 
explicitly excluded from eq.(VI.55). 

number of crossings | 


Let us label random walkers with no U-turns, and weighted by (—1) 
as RW*s. Then as in eq.(VI.37) the terms in S can be organized as 


S= 5" (RW*s with 1 loop) + S-(RW*s with 2 loops) + S"(RW*s with 3 loops) + --- 
= exp aw's with 1 loop) | é 
(V1.56) 
The exponentiation of the sum is justified, since the only interaction between RW*s is the 
sign related to their crossings. As two RW* loops always cross an even number of times, 
this is equivalent to no interaction at all. Using eq.(VI.35), the full Ising free energy is 


calculated as 
InZ = NlIn2+2N IncoshK + S- (RW*s with 1 loop x ae Paes) : (VI.57) 


Organizing the sum in terms of the number of bonds, and taking advantage of the trans- 
lational symmetry of the lattice (up to corrections due to boundaries), 
InZ 


2 — tf * 
ar =In (2cosh K) 72 7 Oe (£)|0) , (VI.58) 


where 


(0|W*(£)|0) =number of closed loops of @ steps, with no U turns, from 0 to 0 | 
VI.59 


x (—1)# of crossings 
The absence of U-turns, a local constraint, does not complicate the counting of loops. 
On the other hand, the number of crossings is a function of the complete configuration of 
the loop and is a non—Markovian property. Fortunately, in two dimensions it is possible to 


obtain the parity of the number of crossings from local considerations. The first step is to 
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construct the loops from directed random walks, indicated by placing an arrow along the 


direction that the path is traversed. Since any loop can be traversed in two directions, 


1 
(O|W*(£)|0) = - S| directed RW* loops of @ steps, no U turns, from 0 to 0 x (—1)"*, 
(VI.60) 
where n, is the number of self—crossings of the loop. We can now take advantage of the 


following topological result: 
Whitney’s Theorem: The number of self-crossings of a planar loop is related to the total 
angle O, through which the tangent vector turns in going around the loop by 


‘e) 
(Reve = (1 + =) . (VI.61) 
20 mod 2 


This theorem can be checked by a few examples. A single loop corresponds to 0 = +27, 
while a single intersection results in O = 0. 
Since the total angle 0, is the sum of the angles through which the walker turns at 


each step, the parity of crossings can be obtained using local information alone as 


' 6) i yy 
(-1)"" = e'*"" = exp [i (: - =) = ae? Lia, (V1.62) 


where 6; is the angle through which the walker turns on the j™ step, leading to 


1 
(0|W*(£)|0) = =5 S- directed RW* loops of @ steps, with no U turns, from 0 to 0 


x exp (5 Ss" local change of angle by the tangent vector 
(VI1.63) 
The angle turned can be calculated at each site, if we keep track of the directions of arrival 
and departure of the path. To this end, we introduce a label yu for the 4 directions going 
out of each site, e.g. w = 1 for right, uw = 2 for up, w = 3 for left, and u = 4 for down. We 


next introduce a set of 4N x 4N matrices generalizing eq.(VI1.39) as 


(Z2Y2, H2|W*(£) |r 1y1, 1) = oS directed random walks of @ steps, with no U turns, 


i é : 
departing (71, y1) along j41, proceeding along jz after reaching (2, y2) x e? a, 
(V1.64) 
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Thus jg specifies a direction taken after the walker reaches its destination. It serves to 
exclude some paths (arriving along -j2), and leads to an additional phase. As in eq.(VI.43), 


due to their Markovian property, these matrices can be calculated recursively as 


(t2y2, H2|W*(£)|e1y1, 11) 
= So (woye, polT*|2’y’, w’) (a'y’, w|W* (6 - 1 |arys, a) 


(V1.65) 


aly! pl 
= (L2Ye, Ma|T"W"*(l — 1)|r1yi, pa) = (x2Yyo, poll eign, in), 
where T* = W*(1) describes one step of the walk. The direction of arrival uniquely 
determines the nearest neighbor from which the walker departed, and the angle between 
the two directions fixes the phase of the matrix element. We can thus generalize eq.(VI.46) 


to a 4x 4 matrix that keeps track of both connectivity and phase between pairs of sites. 


The steps taken can be represented diagrammatically as 


= ‘ad ee 


rots 
i ails (VI1.66) 


he. tr | 


and correspond to the matrix 


(2’y'|T*|xy) = 
(x’, y’|a +1, y) (x’,y'|ce +1,y)et 0 (x!,y'|c+1l,ye7 = 
(x',y'|x,y+1j)e"F (c',y'|t,yt1) (a, y'|a,y+ let 0 . 
0 Gay ea layle =.  Aaly a= 19) (z',y'|e —ly)et 
(x',y'|z,y — 1) et 0 (2',y/|e,y—1)e"F (aoa gt 1) 
(VI.67) 


wheré:< 4.90.9" >= 05 Og 
Because of its translational symmetry, the 4N x 4N matrix takes a block diagonal 


form in the Fourier basis, (xy|qrq,) = e'(Iet+dy¥) //N, ie. 


So (ely, wT" acy, u) (eyldedy) = (u!|T* (a)| 4) (2’y'|dedy) - (VL.68) 
xy 
Each 4 x 4 block is labelled by a wavevector q = (dz, dy), and takes the form 
—i(dy+7) —tdy —i(qy-F 0 
* e 4 e e 4 
T"(q) = i oilde—F oie GEE Nec (VI.69) 
e'(W+7) 0) edu -F ety 
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To ensure that a path that starts at the origin completes a loop properly, the final 
arrival direction at the origin must coincide with the original one. Summing over all 4 such 


directions, the total number of such loops is obtained from 


(0|W*(0)|0) = > (00, p|T*|00, p) = = Se (ey, ulT“ley, w) = tr (Tr). (VL70) 


p=1 LYM 


Using eq.(VI.58), the free energy is calculated as 


InZ lpr’ 1 Te 
ee Si (2 cosh? K) — 5) 7 (0|W*(é)|0) = In (2 cosh? K) - an tt S° 7 
1 
= In (2 cosh? K) + oN tr In (1 = tT*) 
1 
— 2 x 
=In (2cosh? K) + => ee ar ey 
(VI.71) 
But for any matrix M with eigenvalues {AQ}, 
trinM = S Indy =n] J Aq = IndetM. 
Converting the sum over q in eq.(VI.71) to an integral leads to 
In Z 
— =In (2 cosh? K) =f 
1—te-4 —te-te-F 0 —teGet+F) 
if d?q ee —te(Wwtt)  1—te% te“ -F 0 
2) Qn? \~ 0 —telie-F) 1 —telte tee + F) 
—tet(W+F) 0 —tet(Gu-F 1 — te?d 
(VI.72) 


Evaluation of the above determinant is straightforward, and the final result is 


nZ _ 


In (2 cosh” K)+ ae: In (a + #2)” — 2 (1- t?) (cos gx + cos | (VI.73) 
N T 2 dx dy . be 


2 J (2m) 


Taking advantage of trigonometric identities, the result can be simplified to 


In Z 1 [* dqedq 2 
= cae In2+ a ie one In [cosh” (2K) — sinh(2K) (cos gz + cos qy)] « (VI.74) 


While it is possible to obtain a closed form expression by performing the integrals exactly, 


the final expression involves a hypergeometric function, and is not any more illuminating. 
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VI.H_ Critical behavior of the two-dimensional Ising model 


To understand the singularity in the free energy of the two dimensional Ising model 
in eq.(VI.73), we start with the simpler expression obtained by the unrestricted sum over 


phantom loops in sec. VI.F. Specializing eq.(VI1.53) to d = 2, 


to =n (2 cosh? K) -;/ 


dqzdqy 
(27)? 


Apart from the argument of the logarithm, this expression is similar to the exact result. 


In [1 — 2t (cos gz + cos qy)]. (VI.75) 


Is it possible that such similar functional forms lead to distinct singular behaviors? The 
singularity results from the vanishing of the argument of the logarithm at t. = 1/4. In the 


vicinity of this point we make an expansion as in eq.(VI.50), 


t 
Agc(t,q) = (1 — 4t) + tq? + O(q*) Gs G + iS) : (VI.76) 


Cc 


where 6t = t. —t. The singular part of eq.(VI.75) can be obtained by focusing of the 
behavior of the integrand as q — 0, and replacing the square Brillouin zone for the range 


of the integral with a circle of radius A & 27, 


1 (* 2nqdq ot 
sing, = — > In(q? +42 
fing. a ze Bla + z 


2 A 
2 [enemy 
81 te e€ 0 


Only the expression evaluated at g = 0 is singular, and 


1 ot ot 
sing. = — |— ]} 1 —— Ise I. 
fan = == (S) im (F) (VL.78) 


The resulting heat capacity, Cg « 0? fg/07t, diverges as 1/dt. Since eq.(VI.75) is not valid 


(V1.7) 


for t > te, we cannot obtain the behavior of heat capacity on the low temperature side. 


For the exact result of eq.(VI.73), the argument of the logarithm is 
A*(t,q) = (1+ Py — 2t (1 — t?) (cos qz + cos qy) . (VI.79) 
The minimum value of this expression, for q = 0, is 


A*(t,0) = (1+ #2)” — 4t, (1 — #2) = (1 — #2)” + 402 — at, (1 — #2) = (1 — #2 — 24)”. 
(VI.80) 


uae 


Since this expression (and hence the argument of the logarithm) is always non-negative, 
the integral exists for all values of t. As required, unlike eq.(VI.75), the exact result is 


valid at all temperatures. There is a singularity when the argument vanishes for 
#4+2.-1=0, = £t-=-1l+Vv2. (V1.81) 


The positive solution describes a ferromagnet, and leads to a value of K, = In (V2 + 1) (2, 
in agreement with the duality arguments of sec. VI.D. Setting 6¢ = t — t., and expanding 
eq.(VI.79) in the vicinity of q — 0 gives 


A*(t,q)  [(—2t, — 2)6t]” + te(1 — t2)q? + + 


wey (*)| | (VI.82) 


The important difference from eq.(VI.76) is that (dt/t.) appears at quadratic order. Fol- 


met? 


lowing the steps in eqs.(VI.77) and (VI.78), the singular part of the free energy is 


Te 
— | naga In |q7 +4 (=) 
0 Cc 


2 Jy 4? 
G +4 (=) ) In an as a] ; (VL83) 


0 


InZ 
N 


sing. 


1 (2) ét 
=——([(—] ln 
T \ te 


Cc 


+ analytic terms. 


The heat capacity is obtained by taking two derivatives and diverges as C'(dt)sing. = 
Ax In |6t|. The logarithmic singularity corresponds to the limit a = 0; the peak is sym- 


metric, characterized by the amplitude ratio A, /A_ = 1. 

The exact partition function of the Ising model on the square lattice was originally 
calculated by Lars Onsager in 1944 (Phys. Rev. 65, page 117). Onsager used a 2” x 2 
transfer matrix to study a lattice of width L. He then identified various symmetries of this 
matrix which allowed him to diagonalize it and obtain the largest eigenvalue as a function 
of L. For any finite L, this eigenvalue is non-degenerate as required by Frobenius’s theorem. 
In the limit DL — oo, the top two eigenvalues become degenerate at K,. The result in this 
limit naturally coincides with eq.(VI.74). Onsager’s paper is quite long and complicated, 
and regarded as a tour de force of mathematical physics. A somewhat streamlined version 
of this solution was developed by B. Kaufman (Phys. Rev. 76, 1232 (1949)) and is 
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more or less reproduced in chapter 15 of Huang. This exact solution, for the first time, 
demonstrated that critical behavior at a phase transition can in fact be quite different from 
predictions of Landau (mean-field) theory. It took almost three decades to reconcile the 
two results by the renormalization group. 

The graphical method presented in this section was originally developed by Kac and 
Ward (Phys. Rev. 88, 1332 (1952)). The main ingredient of the derivation is the result 
that the correct accounting of the paths can be achieved by including a factor of (-1) 
for each intersection. (This conjecture by Feynman is proved in S. Sherman, J. Math. 
Phys. 1, 202 (1960).) The change of sign is reminiscent of the exchange factor between 
fermions. Indeed, Schultz, Mattis, and Lieb (Rev. Mod. Phys. 36, 856 (1964)) describe 
how the Onsager transfer matrix can be regarded as a Hamiltonian for the evolution of non- 
interacting fermions in one dimension. The Pauli exclusion principle prevents intersections 
of the world-lines of these fermions in two dimensional space-time. 

In addition to the partition function, the correlation functions (jj) can also be 
calculated by summing over paths (see Itzykson and Drouffe, Statistical field theory: 1). 
Since the combination q? + 4(dt/t,)? in eq.(VI.82) describes the behavior of these random 
walks, we expect a correlation length € ~ |t./dt|, i.e. diverging with an exponent v = 1 on 
both sides of the phase transition with an amplitude ratio of unity. The exponents a and 
vy are related by the hyperscaling identity a = 2 — 2v. The critical correlations at t. are 
more subtle and decay as (oj0;), ~ 1/|i—j]", with 7 = 1/4. Integrating the correlation 


functions yields the susceptibility, which diverges as y+ ~ Cx|ét|~7, with y = 7/4, and 


C,/C_ =1. A rather simple expression for the zero field magnetization, 


m = (1—sinh~4(2K))'*, (VI.84) 


was presented without proof by Onsager. A rather difficult derivation of this result was 
finally given by C.N. Yang (Phys. Rev. 85, 808 (1952)). The exponent 3 = 1/8 satisfies 


the required exponent identities. 
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VII. Continuous Spins at Low Temperatures 


VII.A The non-linear c-model 


Previously we considered low temperature expansions for discrete spins (Ising, Potts, 
etc.), in which the low energy excitations are droplets of incorrect spin in a uniform back- 
ground selected by broken symmetry. These excitations occur at small scales, and are 
easily described by graphs on the lattice. By contrast, for continuous spins, the lowest 
energy excitations are long-wavelength Goldstone modes, as discussed in section II.C. The 
thermal excitation of these modes destroys the long-range order in dimensions d < 2. For 
d close to 2, the critical temperature must be small, making low temperature expansions 
a viable tool for the study of critical phenomena. As we shall demonstrate next, such an 
approach requires keeping track of the interactions between Goldstone modes. 


Consider unit n-component spins on the sites of a lattice, i.e. 
3 (i) = (81, 82,---,8n), with |s(i)|? =s?+---+52 =1. (VII.1) 


The usual nearest neighbor Hamiltonian can be written as 


se 


—BGH= KY s(i)- 8) =KY_ (1 = a (VIL.2) 
(ii) (ij) 


At low temperatures, the fluctuations between neighboring spins are small and the differ- 


ence in eq.(VII.2) can be replaced by a gradient. Assuming a unit lattice spacing, 
K d 2 9 
pH = =6Eo= = d°x (V8(x))° , (VII.3) 


where the discrete index i has been replaced by a continuous vector x € R¢%. A cutoff 
of A & z is thus implicit in eq.(VII.3). Ignoring the ground state energy, the partition 
function is 


Z= J? [5 (x)6 (s(x)? —1)Je7 2 J x98)”, (VIL.4) 


A possible ground state configuration is §(x) = (0,---,1). There are n — 1 Goldstone 
modes describing the transverse fluctuations. To examine the effects of these fluctuations 


close to zero temperature, set 


§ (x) = (m71(x),---, Tn-1(X), o(X)) = (7 (x), o(x)), (VIL.5) 
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where 7 (x) is an n — 1 component vector. The unit length of the spin fixes a(x) in terms 


of 7 (x). For each degree of freedom 


[esr -1) = | dit doé (n* +o — 1) 


= didod |(o - Vi—m?) (o + 1=7?)] = [ ———-. , 

i sae 2. /1 = 2 
(VIL.6) 

where we have used the identity d(az) = 6(x)/|a|. Using this result, the partition function 


in eq.(VII.4) can be written as 


Ge Di (x) oS f xP +(V via)" 


vinmx? (VII.7) 
= [ Di (exp {— f a's Fr) 45 (vvi-w) + E In( a-7)]}. 


In going from the lattice to the continuum, we have introduced a density p = N/V = 1/a4 
of lattice points. For unit lattice spacing p = 1, but for the purpose of renormalization 
we shall keep an arbitrary p. Whereas the original Hamiltonian was quite simple, the 
one describing the Goldstone modes 7 (x), is rather complicated. In selecting a particular 
ground state, the rotational symmetry was broken. The nonlinear terms in eq.(VII.7) 
ensure that this symmetry is properly reflected when considering only 7. 

We can expand the nonlinear terms for the effective Hamiltonian in powers of 7 (x), 


resulting in a series 


BH[r (x)] = BHo +U1 +U2+---, (VIL8) 
where 


BHo = = / d’x(Vi )?, (VIL.9) 


describes independent Goldstone modes, while 


bys [atx Fe Vi)? — al (VII.10) 


is the first order perturbation when the terms in the series are organized according to 


powers of T = 1/K. Since we expect fluctuations (1?) x T’, GHo is order of one, the two 
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terms in U/, are order of T; remaining terms are order of T? and higher. In the language 


of Fourier modes, 


K d“q 9). 2 
Ho=F | ae l@l, 


dé dt dé 
uy=- > —aaystTa(a) To(Q2) ™3(3) ™a(—G1 — G2 — Gs) (1 - Gs) 


d 
- 5 | pat an. 
(VII.11) 


For the non-interacting (quadratic) theory, the correlation functions of the Goldstone 


modes are as 
_ $a,6(27)°6"(q + Q') 


(ra(a)ma(a'))o a (VIL.12) 
The resulting fluctuations in real space behave as 
Cee ls a4 (1 (a) I"), = — / Ve déq 1 (n—1) Kao? 4-1-4) 
0 J Qn) : yt Qnjég  K (d — 2) 
(VII.13) 


For d > 2 the fluctuations are indeed proportional to 7. However, for d < 2 they diverge 
as L — oo. This is a consequence of the Mermin—Wagner theorem on the absence of long 
range order in d < 2. Polyakov (1975) argued that this implies a critical temperature 
T. ~ O(d — 2) for such systems, and that an RG expansion in powers of T may provide a 
systematic way to explore critical behavior close to two dimensions. 

To construct a perturbative RG, consider a spherical Brillouin zone of radius A, and 
divide the modes as 7 (q) = 7 <(q) +7 7(q). The modes 7 < involve momenta 0 < |q| < 
A/b, while we shall integrate over the short wavelength fluctuations 77 with momenta in 


the shell A/b < |q| < A. To order of T, the coarse-grained Hamiltonian is given by 
BH [#<] = Vf2 + Ho [#<] + Ui [F< + #7])5 + O(T?), (VII.14) 


where ( )% indicates averaging over 7~*. The term proportional to p in eq.(VIL11) results 
in two contributions, one is a constant addition to free energy (from ((77)?)), and the 
other is simply p(7<)?. (The cross terms proportional to #< -#7 vanish by symmetry.) 
The quartic part of /; generates 16 terms. Nontrivial contributions arise from products 
of two 7 < and two 77. There are three types of such contributions; the first has the form 
dq, d4qy d4 
Ut)9 =2x + > /SS a - * (qi - qs) (VIL15) 
(m2 (au)72 (a2))¢ ™5(as)5 (—a1 — a2 — 93). 
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The integral over the shell momentum qj, is odd and this contribution is zero. (Two similar 
vanishing terms arise from contractions with different indices a and (.) The next term is 
a renormalization of p, arising from 


(ut) ur? )=- = qi d"qz d%q3 


(2734 (qi - qs) 


(n2(a)n3(as)). *S(a2)"§(-an - a2 — as) 


ms A/b ool <@ex f d¢k k2 
A/b (27)4 Kk? 


(VII.16) 


ip Mt ye >< 2 - 
ry, alt (al? x 1-0-4), 


(Note that in general p = N/V = Hs d¢q/(27)?.) Finally, a renormalization of K is 


obtained from 


Kf d4q, d4qz d¢q 
c\> 1 2 3 
(Uq)o =- os — myst (a - qs) 
> 
(m2 (a2)03 (—an aa ae as), me (a1) (as) (VIL.17) 
kph dq a2 </\2 _ ta(d) 
-+ PF <(a)P x 2, 
where P ; 2 : . 
d¢k 1 KgA** (1 — b-7 
Iq(b) = >= = VIL18 
O= | oP (T-9) er 


The coarse-grained Hamiltonian in eq.(VII.14) now equals 


Ta(0) se: d'q Qia4< 2 
Oe) tla) 


GH [#<] =Vofp + VOfp + S > ( 


d*qi a d@q3 7 


a o (G1) 72 (a2) m3 (G3) 7 (—G1 — G2 — Qs) (a1 - Gs) 


2 (Qr)\2¢ 
A/b 
= ef oe (a)! x [1- (1—b7%)] + O(7”). 
(VII.19) 
The most important consequence of coarse graining is the change of the elastic coefficient 
K to 
K=K ( am) (VII.20) 
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After rescaling, x’ = x/b, and renormalizing, 7'(x) = 7 <(x)/¢, we obtain the renor- 
malized Hamiltonian in real space as 
Kbt-2¢2 
—BH' =—V6fP —Vof, - aes jax’ 


HZ (VIL.21) 


Kb AC 


5 [ax (i! (x VE! (x!))? + oe fats n'(x')? + O(T?). 


The easiest method for obtaining the rescaling factor ¢, is to take advantage of the ro- 


tational symmetry of spins. After averaging over the short wavelength modes, the spin 


= (ri EE (ERY) (VIL.22) 


2 (: = () +-0«r*)) (mE, = <P) 


We thus identify 


¢=1- — ) + O(T?) =1- wet) i) + O(T), (VII.23) 


as the length of the coarse-grained spin. The renormalized coupling constant in eq.(VII.21) 


is now obtained from 


K' = DOCK 


— pt? c = 1400) K c + zlu()| (VII.24) 


= —¢2K c = 7 1a(b) oO (=) ; 


For infinitesimal rescaling, b = (1 + 62), the shell integral results in 
Iq(b) = KaA* 768. (VII.25) 


The differential recursion relation corresponding to eq.(VII.24) is thus 


x = (d—2)K —(n—2)Kah*2. (VII.26) 
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Alternatively, the scaling of temperature T = K~, is 


dT 1 dK _ 


eee = —(d—2)T + (n—2)K,A% 7? T". (VII.27) 
It may appear that we should also keep track of the evolution of the coefficients of the 
two terms in U/; under RG. In fact, spherical symmetry ensures that the coefficient of the 
quartic term is precisely the same as K at all orders. The apparent difference between the 
two is of order of O(T?), and will vanish when all terms at this order are included. The 
coefficient of the second order term in U/; merely tracks the density of points and also has 
trivial renormalization. 

The behavior of temperature under RG changes drastically at d = 2. For d < 2, the 
linear flow is away from zero, indicating that the ordered phase is unstable and there is no 
broken symmetry. For d > 2, small T’ flows back to zero, indicating that the ordered phase 
is stable. The flows for d = 2 are controlled by the second ordered term which changes 
sign at n = 2. For n > 2 the flow is towards high temperatures, indicating that Heisenberg 
and higher spin models are disordered. The situation for n = 2 is ambiguous, and it can 
in fact be shown that dT/dé is zero to all orders. This special case will be discussed in 
more detail in the next section. For d > 2 and n > 2, there is a phase transition at the 
fixed point, 

" € 2TE 


where € = d — 2 is used as a small parameter. The recursion relation at order of € is 


ita (VII.29) 


Stability of the fixed point is determined by the linearized recursion relation 


=|-+S > 


ane = =r] OT = [—€ + 2e] OT = oT, = Yee: «x. LVIL3O) 
T* By 


de 


The thermal eigenvalue, and the resulting exponents v = 1/e, and a = 2—(2+€)/e x —2/e, 
are independent of n at this order. 

The magnetic eigenvalue can be obtained by adding a term aha if d’x §(x), to the 
Hamiltonian. Under the action of RG, h’ = b¢¢h = b¥h, with 


1400) (VIL31) 


For an infinitesimal rescaling 


—1 
1+ yp,d€ = (14 dd) (1 _ “Sr Kia" 60) ; (VII.32) 
leading to 
n—-1l n—3 
= fe SS eS 2 VII.33 
Vi = dn —2)° POG =a" ), ( ) 


which does depend on n. Using exponent identities, we find 
n=2+d—2y, = —. (VII.34) 
n—2 


The exponent 77 is zero at the lowest order in a 4 —d expansion, but appears at first order 
in the vicinity of two dimensions. The actual values of the exponents calculated at this 


order are not very satisfactory. 
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VII.B Topological Defects in the XY model 


As stated in the previous section, thermal excitation of Goldstone modes destroys 
spontaneous order in two dimensional models with a broken continuous symmetry. The RG 
study of the non-linear g-model confirms that the transition temperature of n-component 
spins vanishes as T* = 27€/(n — 2) for « = (d — 2) — 0. However, the same RG procedure 
appears to suggest a different behavior for n = 2. The first indication of unusual behavior 
for the two dimensional XY model (n = 2), appeared in an analysis of high temperature 
series by Stanley and Kaplan in 1966. The series results strongly suggested the divergence 
of susceptibility at a finite temperature, seemingly in contradiction with the absence of 
symmetry breaking. It was indeed this contradiction that led Wigner to explore the pos- 
sibility of a phase transition without symmetry breaking. The Z2 lattice gauge theory, 
discussed in sec. VI.E as the dual of the three dimensional Ising model, realizes such a 
possibility. The two phases of the Z2 gauge theory are characterized by different func- 
tional forms for the decay of an appropriate correlation function (the Wilson loop). We 
can similarly examine the asymptotic behavior of the spin-spin correlation functions of the 
XY model at high and low temperatures. 

A high temperature expansion for the correlation function for the XY model on a 


lattice is constructed from 


N 27 
= = 1 do; K)~ _._ cos(@;—8;) 
(50° 5+) = (cos(89 — Oy)) = Z II (| _ cos (89 — ,) € <i,j> 
_1n = cos (09 — %,) [] [1+ K cos (6; — 0;) + O (K’)| 
Z 31 Wo 20 ee <i,j> ~~ . 


(VII.35) 
The expansion for the partition function is similar, except for the absence of the factor 
cos(69 — 6,). To the lowest order in kK, each bond on the lattice contributes either a factor 


of one, or K cos(6; — @;). Since, 
27 do, 
i — cos (0; — 02) = 0, (VII.36) 
9 27 


any graph with a single bond emanating from an internal site vanishes. For the numerator 
of eq.(VII.35) to be non-zero, there must be bonds originating from the external points at 


O and r. Integrating over an internal point with two bonds leads to 
27 
dé 1 
| oe cos (6; — 02) cos (82 — 63) = 5 008 (0; — 03). (VII.37) 
0 1 
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The first graph that contributes to eq.(VII.35) is the shortest path (of length r) connecting 


the points 0 and r. Since integrating over the end points gives 


27 
d0od6y 2 1 
ee IL 
i (ony cos (69 — Oy) 5 (VII.38) 
each bond along the path contributes (K/2). (In constructing graphs for the partition 


function, there is an additional factor of 2 for every loop.) Thus to lowest order 


Kk\" 1 
oa ye oe (| (8 ith ———— VIL.39 
Wo m)~(F) er, with expo, (VIL9) 
and the disordered high temperature phase is characterized by an exponential decay of 
correlations. (This is clearly quite generic in all spin systems.) 

At low temperatures, the cost of small fluctuations around the ground state is ob- 
tained by a quadratic expansion, which gives K { d?x(V6)?/2 in the continuum limit. The 


standard rules of Gaussian integration yield, 
. 1 
Ae (eil6o—F.) = Rexp -5 ( (G0 zs a.)") (VIL.40) 
In two dimensions, the Gaussian fluctuations grow as 


; ((60 — 6)?) = — In (=) . (VII.41) 


where a is a short distance cut-off (of the order of the lattice spacing). Hence, at low 


temperatures, 
a 


(50-3) © (=) oe , (VIL.42) 
i.e. the decay of correlations is algebraic rather than exponential. A power law decay of 
correlations implies self-similarity (no correlation length), and is usually associated with 
a critical point. Here it arises from the logarithmic growth of angular fluctuations, which 
is specific to two dimensions. 

The distinct asymptotic decays of correlations at high and low temperatures allows 
for the possibility of a finite temperature phase transition separating the two regimes. 
However, the arguments put forward above are not specific to the XY model. Any contin- 
uous spin model will exhibit exponential decay of correlations at high temperature, and a 
power law decay in a low temperature Gaussian approximation. To show that the Gaus- 


sian behavior persists at finite temperatures, we must prove that it is not modified by the 
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additional terms in the gradient expansion. Quartic terms, such as [ d’x(V0)*, generate 
interactions between the Goldstone modes. The relevance of these interactions was probed 
in the previous section using the non-linear o-model. The zero temperature fixed point in 
d = 2 is unstable for all n > 2, but apparently stable for n = 2. (There is only one branch 
of Goldstone modes for n = 2. It is the interactions between different branches of these 
modes for n > 2 that leads to instability towards high temperature behavior.) The low 
temperature phase of the XY model is said to posses quasi—long range order, as opposed 
to true long range order that accompanies a finite magnetization. 

What is the mechanism for the disordering of the quasi—long range ordered phase? As 
the RG study suggests that higher order terms in the gradient expansion are not relevant, 
we must search for other operators. The gradient expansion describes the energy cost of 
small deformations around the ground state, and applies to configurations that can be 
continuously deformed to the uniformly ordered state. Kosterlitz and Thouless (1973) 
suggested that the disordering is caused by topological defects that can not be regarded 
as simple deformations of the ground state. Since the angle describing the orientation 
of a spin is undefined up to an integer multiple of 27, it is possible to construct spin 
configurations for which in going around a closed path the angle rotates by 27n. The 
integer n is the topological charge enclosed by the path. Because of the discrete nature of 
this charge, it is impossible to continuously deform to the uniformly ordered state in which 
the charge is zero. (More generally, topological defects arise in any model with a compact 
group describing the order parameter.) 

The elementary defect, or vortex, has unit charge. In completing a circle centered 
on the defect the orientation of the spin changes by +27. If the radius r, of the circle 
is sufficiently large, the variations in angle will be small and the lattice structure can be 
ignored. By symmetry, V@ has a uniform magnitude, and points along the circle (i.e. 


perpendicular to the radial vector). The magnitude of the distortion is obtained from 


dé 
$ vo -ds = La = 271, = sae (VII.43) 
ds ds =r 
Since V@ is a radial vector, it can be written as 
Ge ae ee ee ‘ 
Vo =n( = +35,0) n# x 8 = —nV x (ZInr). (VII.44) 


Here, * and Z are unit vectors respectively in the plane and perpendicular to it, and d x b 
indicates the cross product of the two vectors. This (continuum) approximation fails close 


to the center (core) of the vortex, where the lattice structure is important. 
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The energy cost of a single vortex of charge n has contributions from the core region, 
as well as from the relatively uniform distortions away from the center. The distinction 
between regions inside and outside the core is arbitrary, and for simplicity, we shall use a 


circle of radius a to distinguish the two, i.e. 


K 
BEn =BEO(a) + a ‘i d°x(V0)? 


. (VII.45) 


=se%(a) +> f 


a 


(Qnrdr) (ey = B€°(a) + 1Kn?In (=) 


The dominant part of the energy comes from the region outside the core, and diverges 
with the size of the system, L. The large energy cost of defects prevents their spontaneous 
formation close to zero temperature. The partition function for a configuration with a 


single vortex is 


CATA (£) ex [-se%(w = eiene Th ()] 506s a | (VIL46) 


where (L/a)? results from the configurational entropy of possible vortex locations in an 
area of size L?. The entropy and energy of a vortex both grow as In L, and the free energy 
is dominated by one or the other. At low temperatures (large K’) energy dominates and 
Z,, a measure of the weight of configurations with a single vortex, vanishes. At high 
enough temperatures, K < K, = 2/(mn7), the entropy contribution is large enough to 
favor spontaneous formation of vortices. On increasing temperature, the first vortices to 
appear correspond to n = +1 at K. = 2/7. Beyond this point there are many vortices in 
the system, and eq.(VII.46) is no longer useful. 

The coupling K, = 2/7 is a lower bound for the stability of the system to topological 
defects. This is because pairs (dipoles) of defects may appear at larger couplings. Consider 
a pair of charges +1 at a separation d. The distortions far away from the dipole center, 


r > d, can be obtained by superposing those of the individual vortices, and 


V6 = V4 + V6_ wiv (*), (VIL.47) 


decays as d/r?. Integrating this distortion leads to a finite energy, and hence dipoles 
appear with the appropriate Boltzmann weight at any temperature. The low temperature 
phase should thus be visualized as a gas of tightly bound dipoles. The number (and size) of 


dipoles increases with temperature, and the high temperature state is a plasma of unbound 
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vortices. The distinction between the two regimes can be studied by examining a typical 
net topological charge, Q(¢), in a large area of dimension @ >> a. The average charge 
in always zero, while fluctuations in the low temperature phase are due to the dipoles 
straddling the parameter, i.e. (Q(2)?) x . In the high temperature state, charges of 
either sign can appear without restriction and (Q(é)?) « €?. (Note the similarity to the 
distinct behaviors of the Wilson loop in the high and low temperature phases of the Z2 
gauge theory.) 

To describe the transition between the two regimes, we need to properly account 
for the interactions between vortices. The distortion field uv = V@, in the presence of a 
collection of vortices is similar to the velocity of a fluid. In the absence of vorticity, the 
flow is potential, i.c. U = Up = Vo, and V x tio = 0. The topological charge can be related 


to V x u by noting that for any closed path, 
gi ds = [(@xa) ‘Vxu , (VII.48) 


where the second integral is over the area enclosed by the path. Since the left hand side is 


an integer multiple of 27, we can set 


Vx a= 2025 7; 6(x — x), (VII.49) 


describing a collection of vortices of charge {n;} at locations {x;}. The solution to 
eq.(VII.49) can be obtained by setting w = tip — V x (Zw), leading to 


Vxi=2V7p, => Wp = In} 0; 6? — x). (VIL.50) 
Thus w behaves like the potential due to a set of charges {27n,;}. The solution, 


w(x) = dni In (|x — x;|) , (VIL.51) 


is simply a superposition of the potentials as in eq.(VII.44). 


Any two dimensional distortion can thus be written as 
é=tot+t% =Ve-V x (Zy), (VII.52) 
and the corresponding “kinetic energy”, GH = K f d?x |u|?/2, decomposed as 
BH = / d?x (vo) —2V¢6-V x (Sb) + (V x Sw)? } . (VIL.53) 
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The second term vanishes, since following an integration by parts, 
- | @xvo-¥ x (2) = [ @xov -(V x 2y), (VII.54) 


and V-V xu = 0 for any vector. The third term in eq.(VII.53) can be simplified by noting 


that Vw = (O20), Oyw, 0), and V x (2) = (—O,~, Oz, 0), are orthogonal vectors of equal 
length. Hence 


BH, = = > [ex (V x 2)? = sf x (Vy)? = —— 5 | exuvs, (VII.55) 


where the second identity follows an integration by parts. Equations (VII.50) and (VII.51) 


now result in 


K 
ory =- > f ax (Xe. In (|x — x;]) ) Beast (x — x;) 


=_ an? KS nin, C(x; = x; 


a,j 


(VIL.56) 


where C(x) = In(|x|)/27 is the two dimensional Coulomb potential. There is a difficulty 
with the above result for 7 = 7 due to the divergence of the logarithm at small arguments. 
This is a consequence of the breakdown of the continuum treatment at short distances. 


The self-interaction of a vortex is simply its core energy GE°, and 


= DBE —4n?K S- igttj CR = Ky). (VII.57) 

tJ 
The configuration space of the XY model close to zero temperature can thus be parti- 
tioned into different topological segments. The degrees of freedom in each segment are the 
charges {n;}, and locations {x;}, of the vortices, in addition to the field ¢(x) describing 


spin waves. The partition function of the model can thus be approximated as 


Z= n/ s a; K S~ cos(6; — 9;) 


<i,j> 


x [ Peo & f dx (Vo)? 3 fexe® 2 She (+40? KD, ninjC (xi —x5) (VII.58) 
{ni} 
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where Z,.w. is the Gaussian partition function of spin waves, and Zg is the contribution of 
vortices. The latter describes a grand canonical gas of charged particles, interacting via the 
two dimensional Coulomb interaction. In calculating the Hamiltonian GH, in eq.(VII.55), 
we performed an integration by parts. The surface integral that was ignored in the process 
in fact grows with system size as (}°,n;) nL. Thus only configurations that are overall 
neutral are included in calculating Zg. We further simplify the problem by considering 
only the elementary excitations with n; = +1, which are most likely at low temperatures 


due to their lower energy. Setting yo = exp [-BE ails 


lore) N 
ZQ = a WS [Tex exp An? KS giqgC (Xi =) (5 (VII.59) 
N=0 w=1 


t<j 


where g; = +1, and 5°, q@ = 0. 
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VII.C Renormalization Group for the Coulomb Gas 


The two partition functions in eq.(VII.58) are independent and can be calculated 
separately. As the Gaussian partition function is analytic, any phase transitions of the 
XY model must originate in the Coulomb gas. As briefly discussed earlier, in the low 
temperature phase the charges appear only in the small density of tightly bound dipole 
pairs. The dipoles dissociate in the high temperature phase, forming a plasma. The 
two phases can be distinguished by examining the interaction between two external test 
charges at a large separation X. In the absence of any internal charges (for yo = 0) in 
the medium, the two particles interact by the bare Coulomb interaction C(X). A finite 
density of internal charges for small yo partially screens the external charges, and reduces 
the interaction between the test charges to C(X)/e, where « is an effective dielectric 
constant. There is an insulator to metal transition at sufficiently large yo. In the metallic 
(plasma) phase, the external charges are completely screened and their effective interaction 
decays exponentially. 

To quantify the above picture, we shall compute the effective interaction between two 
external charges at x and x’, perturbatively in the fugacity yo. To lowest order, we need 


to include configurations with two internal charges (at y and y’), and 


e BV (x-x’) _ eT Am? KC(x-x’) x 


[1 +08 f d2ydty! ear ROW 9) H40 KCl y)—C(x—y’)—C(y’ —x)+C(x’-y’)] + O(y8)| 
[1 ie ye f d2yd2y! e—47? KC(y—y’) + O(yG)| 
ea eft KC(x-x’) c +48 [ Pydty’ wr Kco-y) Cue ee. as 1) abs O(ys) 
(VII.60) 
where D(x,x’;y,y’) is the interaction between the internal and external dipoles. The 
direct interaction between internal charges tends to keep the separation r = y’ — y small. 


Using the center of mass R = (y + y’)/2, we can change variables to y = R — r/2 and 
y’ =R+4r/2, and expand the dipole-dipole interaction for small r as 


D(x,x'sy,y’) =C (x-R+=) -C(x-R-=)-C(x’'-R+<2)+C(x'-R-=) 


=—r-VRC(x—R)+r-VRC(x’ — R) + O(r’). 
(VIL61) 


To the same order 
ef KD Gx syy') 1 = — 4a? Ky» Vp (C(x — R) — C(x’ — R)) 


(VII.62) 
+8n'K? [r- Vp (C(x — R) — C(x’ — R))]? + O(r°). 


128 


After the change of variables [ d?yd?y’ > [ d?rd?R, the effective interaction becomes 


e BV (x-x’) = aa it +48 f Pr@RetKOCy 


( — 4? Kr- Vp (C(x — R) — C(x! — R)) + 80'K? [r- Vp (C(x — R) — C(x’ — R)))? 


(VII.63) 
Following the angular integrations in d?r, the term linear in r vanishes, while the angular 
average of (r- VRC)? is r2|VRC|?/2. Hence eq.(VII.63) simplifies to 


+ O(r*)) + O(y8) 


e BV (x-x’) = e 4M KC(x-x’) y 


c + y2 [Onrarje KOM ant? [err (C(x — R) — C(x’ — R)) is + or’) , 


(VII.64) 
The remaining integral can be evaluated by parts, 
/ @R. [VR (C(x — R) — C(x’ — R))]? 
= — | d?R (C(x — R) — C(x’ — R)) (V2C(x — R) — V?C(x’ —R 
[OR (Coe B)— OG! — RB) (V’OH—R)- VCH RB) args) 


Z - [er (C(x — R) — C(x’ — R)) (62(x — R) — 6°(x’ — R)) 
= 2C(x — x’) — 2C(0). 


The short distance divergence can again be absorbed into a proper cutoff with C(«) — 
In(a/a)/27, and 


eo BV (x-x’) = e 4m? KC(x—x’) c x 161° Ky2C(x = x’) fae ait tus) ; 
(VII.66) 


The second order term can be exponentiated to give an effective interaction BV(x — x’) = 
An? KegC(x — x’), with 
Keg = K — 4n? K?yoae™ / drr? 7" * 4 O(ys). (VII.67) 


We have thus evaluated the dielectric constant of the medium, ¢ = K/Keg, pertur- 
batively to order of y2. However, the perturbative correction is small only as long as 


the integral in r converges at large r. The breakdown of the perturbation theory for 
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K < K. = 2/z, occurs precisely at the point where the free energy of an isolated vortex 
changes sign. This breakdown of perturbation theory is reminiscent of that encountered 
in the Landau—Ginzburg model for d < 4. Using the experience gained from that problem, 
we shall reorganize the perturbation series into a renormalization group for the parameters 
K and yo. 

To construct an RG for the Coulomb gas, note that the partition function for the 
system in eq.(VII.59), involves two parameters (K, yo), and has an implicit cutoff a, related 
to the minimum separation between vortices. As discussed earlier, the distinction between 
regions inside and outside the core of a vortex are arbitrary. Increasing the core size to 
ba modifies not only the core energy, hence yo, but also the interaction parameter K. 
The latter is a consequence of the change in the dielectric properties of the medium due 
to dipoles of separations between a and ba. The change in fugacity is obtained from 


eq.(VII.46) by changing a to ba as 
Yo(ba) = b?-™*yo(a) (VII.68) 


The modified Coulomb interaction due to dipoles of all size is given in eq.(VII.67). (The 
perturbative calculation at order of y@ incorporates only dipoles.) From dipoles in the size 


range, a to ba, we obtain a contribution 


ba 
K=K c — @ntx) | (27rdr) (ye#r'xO) | r| ; (VII.69) 
where the terms are grouped so as to make the similarity to standard computations of 
the dielectric constant apparent. (The probability of creating a dipole is multiplied by its 
polarizability; the role of @ = (kgT)~' is played by 277K.) 
By choosing an infinitesimal b = e6 = 1 + @, eq.(VII.69) is converted to 
dK 


Ts —4n? K7a*tye + O(yg) - (VII.70) 


Including the fugacity, the recursion relations are 


dk~} 
7 asta’ ys + O(ys) 


d 
- =(2— 7K) yo + O(ys) 


(VIL.71) 


y] 


originally obtained by Kosterlitz in 1975. While dk —!/dé > 0, the recursion relation for yo 


changes sign at K7+ = 7/2. At smaller values of K~! (high temperatures) yo is relevant, 
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while at lower temperatures it is irrelevant. Thus the RG flows separate the parameter 
space into two regions. At low temperatures and small yo, flows terminate on a fixed line 
at Yo = 0 and Keg > 2/7. This is the insulating phase in which only dipoles of finite size 
occur. (Hence the vanishing of yo under coarse graining.) The strength of the effective 
Coulomb interaction is given by the point on the fixed line that the flows terminate on. 
Flows not terminating on the fixed line asymptote to larger values of K~! and yo, where 
perturbation theory breaks down. This is the signal of the high temperature phase with 
an abundance of vortices. 

The critical trajectory that separates the two regions of the phase diagram flows to a 
fixed point at (K 7! = 1/2, yo = 0). To find the critical behavior at the transition, expand 
the recursion relations in the vicinity of this point by setting 2 = K~'—7/2, and y = yoa?. 


To lowest order, eqs.(VII.71) simplify to 


“ =4ry? + O(ry”, y*) 

ee (VII.72) 
— =—ryt+ O(27y,y°) 

dé ot 


The recursion relations are inherently nonlinear in the vicinity of the fixed point. This is 
quite different from the linear recursion relations that we have encountered so far, and the 
resulting critical behavior is non-standard. First note that eqs.(VII.72) imply 

dx? 3.2. pidge d 5 


eee a5 Ns, Me ee Dg 
7 8r°y x2 = 7 Th = 7 my”) 0, => aw—-any’=c.  (VII73) 


The RG flows thus proceed along hyperbolas characterized by different values of c. For 
c > 0, the focus of the hyperbola is along the y-axis, and the flows proceed to (x, y) — oo. 
The hyperbolas with c < 0 have foci along the x-axis, and have two branches in the half 
plane, y > 0: the branches for x < 0 flow to the fixed line, while those in the x > 0 
quadrant flow to infinity. The critical trajectory separating flows to zero and infinite y 
corresponds to c = 0, i.e. x. = —7*y-. Therefore, a small but finite fugacity yo reduces 
the critical temperature to K7 = 1/2 — r?yoa? + O(y@). 

In terms of the original XY model, the low temperature phase is characterized by a line 
of fixed points with Keg = lime... K(@) > 2/a. There is no correlation length at a fixed 
point, and indeed the correlations in this phase decay as a power law, (cos (0; — cos 09)) ~ 
1/r", with n = 1/(27 Keg) < 1/4. Since the parameter c is negative in the low temperature 
phase, and vanishes at the critical point, we can set it to c = —b?(T, — T) close to the 


transition. In other words, the trajectory of initial points tracks a line (xo0(T), yo(Z)). 
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The resulting c = x3 — r*yé x (T. — T) is a linear measure of the vicinity to the phase 
transition. Under renormalization, such trajectories flow to a fixed point at y = 0, and 


x = —b/T, — T. Thus in the vicinity of transition, the effective interaction parameter 
2 A 2 Ab 
Ket = —-- li C)=-+-—3V0T.-T, VIL74 
an OT ee foo Re ( ) 


has a square root singularity. 

The stiffness Keg, can be measured in experiments on superfluid films. In the super- 
fluid phase, the order parameter is the condensate wavefunction 7(x) = Ve’. Variations 
in the phase 0, lead to a superfluid kinetic energy, 


2 2 
i [aise (-* ~) y= oe [atxcvoy, (VII.75) 


2m 2m 


where m is the particle (helium 4) mass. The corresponding XY model has a stiffness 
K =h’p,/(2mkgT), where p, = VW? is the superfluid density. The density p, is measured 
by examining the changes in the inertia of a torsional oscillator; the superfluid fraction, 
Psm, experiences no friction and does not oscillate. Bishop and Reppy (1978) examined 
ps for a variety of superfluid films (of different thickness, helium 3 concentration, etc.) 
wrapped around a torsional cylinder. They constructed the effective stiffness K as a 
function of temperature, and found that for all films it undergoes a universal jump of 
2/7 at the transition. The behavior of K for T < T, was consistent with a square root 
singularity. 

Correlations decay exponentially in the high temperature phase. How does the corre- 
lation length €, diverge at T.? The parameter c = x? — n4y? = b?(T — T.) is now positive 


all along the hyperbolic trajectory. The recursion relation for x, 


a2 andy? = 4 2 POT T.)), (VII.76) 
de 1 


can be integrated to give 


dx 4 1 x 
a RG ots) ee a 
e+e.) mn Qn (=z) 


The contribution of the initial point to the left hand side of the above equation can be left 


4 
=—£.  (VIL.77) 
Tv 


out if ro x (T —T.) < 1. The integration has to be stopped when xz(¢) ~ y(@) ~ 1, since 
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the perturbative calculation is no longer valid beyond this point. This occurs for a value 


of 
7 7 


~ Ab/T — Te 2 
where we have used arctan(1/b/T — T.) & arctan(oo) = 7/2. The resulting correlation 
length is 


et (VII.78) 
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Unlike any of the transitions encountered so far, the divergence of the correlation length 
is not through a power law. This is a consequence of the nonlinear nature of the recursion 
relations in the vicinity of the fixed point. 

Vortices occur in bound pairs for distances smaller than €, while there can be an excess 
of vortices of one sign or the other at larger separations. The interactions between vor- 
tices at large distances can be obtained from the Debye-Hiickel theory of polyelectrolytes. 
According to this theory the free charges screen each other leading to a screened Coulomb 
interaction, exp(—r/€)C(r). On approaching the transition from the high temperature 


side, the singular part of the free energy, 


Ting ce oe (- (VII.80) 


12 
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has only an essential singularity. All derivatives of this function are finite at T.. Thus the 
predicted heat capacity is quite smooth at the transition. Numerical results based on the 
RG equations (Berker and Nelson) indicate a smooth maximum in the heat capacity at a 
temperature higher than T)., corresponding to the point at which the majority of dipoles 
unbind. 

The Kosterlitz-Thouless picture of vortex unbinding has found numerous applications 
in two dimensional systems such as superconducting and superfluid films, thin liquid crys- 
tals, Josephson junction arrays, electrons on the surface of helium films, etc. Perhaps more 
importantly, the general idea of topological defects has had much impact in understanding 
the behavior of many systems. The theory of two dimensional melting developed in the 


next section is one such example. 


133 


VII.D Two Dimensional Solids 


The transition from a liquid to solid involves the breaking of translational symmetry. 
Here we shall examine the types of order that are present in a two dimensional solid. The 


destruction of such order upon melting is discussed next. 


(a) A low temperature treatment of the problem starts with the perfect solid at T = 0. 
The equilibrium configuration of atoms forms a lattice, ro(m,n) = me, + neg, where 
€; and eg are basis vectors, and {m,n} are integers. At finite temperatures, the atoms 


fluctuate away from their equilibrium position, moving to 
r(m,n) = ro(m,n) + u(m, n). (VII.81) 


The low temperature distortions do not vary substantially over nearby atoms, enabling us 
to define a coarse-grained distortion field u(x), where x = (x1, £2) is treated as continuous, 
with an implicit short distance cutoff of the lattice spacing a. The distortion u, is the analog 


of the angle 0, is the XY model. 


(b) Due to translational symmetry, the energy depends only on the strain matriz, 


a 


7 (qiu; + qjus). (VIL82) 


1 
aj (X) = 5 (Ojuj + Ojui), = ij () 


The elastic energy must respect the symmetries of the underlying lattice. For simplicity, 
we shall consider the triangular lattice whose elastic energy is fully isotropic, i.e. invariant 
under all rotations (see Landau and Lifshitz, Theory of Elasticity). In terms of the Lamé 


coefficients A and pL, 


1 
BH == [ex [2uuij wig + AUeU;5] 


2 
5 (VIL83) 
=5 | ae ln + + Naw), 


The rotational invariance of energy is ensured by the implicit sum over the indices (i, 7) 
in the above expression. In the Fourier representation, the energy depends on the quan- 
tities q?, Jul?, and (q-u)? which are clearly independent of rotations. For other lattices, 
there are more elastic coefficients, since the energy should only be invariant under lattice 
rotations. For example, the symmetry of a square lattice permits a term proportional to 


(eva + 9%). 
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(c) The Goldstone modes associated with the broken translational symmetry are phonons, 
the normal modes of vibrations. Equation (VII.83) supports two types of normal modes: 
(i) Transverse modes with q 1 u, have energy 


Le d?q 1 
ae a (Qr)? gup, => — {\ur(a@)l’) = ae (VII.84) 
(ii) Longitudinal modes with q || u, with energy 
2u+2 d?q 1 
BH, = 5 / (Qn)? qruz,, —— (\uz(q)|?) = CDi (VII.85) 


Combining the above results, we find that a general correlation function is given by 


; (27)?5?(q + q’) Utd Gd; 
(ui(q)uj(a’)) = —— ae 0ij — ee z : (VII.86) 
The extent of fluctuations in position space is given by 
*q 2-2 : iQi 
({u(e) - (oy?) = f 9, Pea) |, - LS 2] 
“ e (VIL87) 


3+ In(|x|/a) 
pu( 2p + A) T < 
The unbounded growth of fluctuations with |x| indicates the destruction of long-range 
translational order in two dimensions. A similar calculation in higher dimensions results 
in a finite maximum extent of fluctuations proportional to temperature. The Lindeman 
criterion identifies the melting point heuristically as the temperature for which these ther- 
mal fluctuations reach a fraction (roughly 0.1) of the perfect lattice spacing. According to 
this criterion, the two dimensional solid melts at any finite temperature. This is of coarse 
another manifestation of the general absence of true long-range order in two dimension. Is 
some form of quasi-long-range order, similar to the XY model possible in this case? 
(d) The order parameter describing broken translational symmetry is p(x) = e’FT™, 
where G is any reciprocal lattice vector. Since G-rg is an integer multiple of 27 by definition, 
pa = 1 at zero temperature. Due to fluctuations (pq@(x)) = (ese \, decreases at finite 
temperatures, and its correlations decay as 


(pa(x)pG(0)) = (efS-(ule) (0) | 


ft (cia = 1) (cia'™ = 1) (ua(a)usta)} 


>| d’q 2—2cosq:x (< ptr cat 
ORD tl igi — ee ee er 


(2m)? g wo p2ut+A)  @ 
ay G?(3u+ A) In(|x|/a) | — fa \"9 
SPL QpQp tA) 2m 7 (a) 


(VIL88) 
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The integration over q is performed after replacing (G - q)? with the angular average of 
G?q?/2. This is not strictly correct because of the presence of the term cosq- x which 
breaks the rotational symmetry. Nevertheless, this approximation captures the correct 
asymptotic growth of the integral. Correlations thus decay algebraically with a power 

_ G3 +A) 


ee ee VIL89 
An p(2Qu + A) ( ) 


NG 


The translational correlations are measured in diffraction experiments. The scattering 


amplitude is the Fourier transform of pg, and the scattered intensity at a wave-vector q is 
proportional to the structure factor 
2 
Starr) =Ny~ (ee elen ray) 
m,n m,n (VII.90) 


aN ss eit (Fo(m,n)—F0(0,0)) Cera) 


min 


S(a) = (|A(a)|*) -{ 


where N is the total number of particles. At zero temperature, the structure factor is 
a set of delta-functions (Bragg peaks) at the reciprocal lattice vectors. Even at finite 
temperature the sum is zero due to varying phases, unless q is in the vicinity of a lattice 


vector G. We can then replace the sum with an integral, and 


ie i! 
S(q) x N @xeila-@)= {| @ )\ 2 Nn Yy——__.. 11.91 
(q) | xe Ix| d jq— G)2-76 (VII.91) 


The Bragg peaks are now replaced by power law singularities. The strength of the diver- 
gence decreases with temperature and increasing |G|. The peaks corresponding to suffi- 
ciently large |G| are no longer visible, and gradually more of them disappear on increasing 


temperature. In three dimensions, the structure factor is still a set of delta-functions, but 


_ GP _5yt+2a 
127a p(2u+A) }° 


with magnitudes diminished by the so called Debye-Waller factor of exp ( 


(e) The crystal phase is also characterized by a broken rotational symmetry. We can 
define an orientational order parameter U(x) = e®°), where 0(x) is the angle between 
local lattice bonds and a reference axis. (The factor of 6 results from the equivalence of the 
6 possible directions on the triangular lattice. The appropriate choice for a square lattice 
is e9(*).) The order parameter has unit magnitude at T = 0, and is expected to decrease 
due to fluctuations at finite temperature. The distortion u(x) leads to a change in bond 


angle given by 


_ 1 du, Our i, Ses 
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The decay of orientational fluctuations is now calculated from 
(U(x) U*(0)) = ( ef6(0%e)-0(0)) } 
6 2 
= exp -5 — (Iv x u(x) — V x u(0)| ) 


9 d? d? ‘ jax iq’-x 
=exp {5 f TASH (cle — 1)(€* — tNegeenrnrpaly (uala)uw(a))} 


9 d?q 2 2 2 
sen {3 f S4 2-20080:0 (a? (|u(q)I2) — ((a- a) »)} 


9 d?q 9 
= 000 | a | aye Boe) | wee (—a) 
(VII.93) 


(Note that { d?q/(27)? is the density n = N/L? = 1/a?.) The final result is independent 
of x at large distances, and asymptotes to a constant that decays exponentially with 
temperature (since w « 1/7’). The two dimensional solid is thus characterized by quasi- 


long-range translational order, and true long-range orientational order. 


VII.E Two Dimensional Melting 


There are qualitative similarities between the melting of a two dimensional solid and 
the disordering of the XY spin system. Expanding on the work of Kosterlitz and Thouless, 
Halperin and Nelson, and independently, Young, developed a theory of two dimensional 


melting. The resulting, so called KTHNY theory, is briefly sketched here. 


(a) The topological defects of a solid are dislocations. A single dislocation corresponds 
to an extra half lattice plane, and is characterized by performing a Burger’s circuit: The 
circuit is a series of steps from site to site that on a perfect lattice returns to the starting 
point (e.g. 5 steps to the right, 4 down, followed by 5 left, and 4 up on a square lattice.) If 
the circuit encloses a region with dislocations, it will fail to close. The difference between 
the initial and final sites of the circuit defines b, the Burger’s vector of the dislocation. 
This closure failure is allowed, since the distortion field u(x), is defined only up to a lattice 


vector. In the continuum limit this degeneracy is described by 


¢ Vua ds = b%. (VII.94) 
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Thus each component of the distortion field behaves like the angle (x) is the XY model. 
Using this analogy, the distortion field due to a collection of dislocations {b;}, at locations 


{x;}, has a singular part 


Vi =—-V x (: S- a In |x — «) : (VII.95) 


While w*, correctly describes the singularities due to the presence of vortices, it is not 
a local equilibrium configuration, i.e. not a solution to 2utag + WXeguy, = 0. Allowing 
the particles to equilibrate in the presence of the dislocations leads to a solution Ug, which 
differs from uZ, only by a regular part. The total strain field can now decomposed as Uag = 
dag + Uag, Where Uag results from dislocations, and dag is the contribution of phonons. 
After substituting this form in the elastic energy GH = f d?x [2uuaguag + AUaatggl /2, 
manipulations similar to the XY model lead to GH = GHo + GH, where 


Ho = 5 / x 2ubapbag + boaeal (VIL96) 
and 
BH, =—K Y~ b?b?Caa(xi — xj) — S— Infyo(ba)]. (VIL.97) 
<j a 


The dislocations behave as a grand-canonical gas of particles with vector charges {b;}. 
The fugacity yo comes from the core energy of each dislocation. The charges interact via 
a vectorial generalization of the Coulomb interaction, 


Cag(x) = = oan (=) = — (VIL.98) 


The strength of the interaction in eq.(VII.97) is Ka? = 2u(u+A)/(2u+ A). Implicit in 
the above derivation is an overall neutrality condition, 5), b; = 0. 

The bare interaction between dislocations is screened by other dislocations in the 
medium. The effective strength can be calculated perturbatively and diverges for 
nKa? <2. This signals the unbinding of paired neutral dislocations. Under coarse grain- 


ing, the parameters of the theory evolve as 


dk! 
aa (VIL.99) 
So = (2—nK) y+ Dy’, 
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where A, B, and D are constants. In contrast to the XY model, additional terms appear 
at order of yg. This is because in a scalar Coulomb gas the only neutral configurations 
have an even number of particles. For dislocations on a triangular lattice, it is possible to 
construct a neutral configuration from three dislocations at 120° angles. 

The low temperature phase maps onto a line of fixed points characterized by renor- 
malized Lamé coefficients, zp and AR, and yj = 0. At high temperatures yo and K both 
diverge, indicating the vanishing of the shear modulus tr, and a finite correlation length €. 
On approaching the transition from the low temperature side, the effective shear modulus 
undergoes a discontinuous jump, with a singular behavior, wr & fe +c(T; — T)”. The 
correlation length diverges from the high temperature side as € ~ aexp (c’ /(T — Te)” iF 
Due to the cubic terms appearing in eqs.(VII.99), the value of 7 = 0.36963 --- is different 
from the 1/2 that appears in the XY model, indicating a difference between transitions in 


vector and scalar Coulomb universality classes. 


(b) Does the vanishing of the shear modulus imply that the high temperature phase 
with unbound dislocations is a liquid? As discussed earlier, the crystalline phase has both 
translational and orientational order. A distortion u(x) results in a rotation in bond angles 


according to eq.(VII.92). The net rotation due to a collection of dislocations is 


~ 1 is wk b; - (x — x;) 
Lies a II.100 
a a le 2. er (V ) 
In terms of a continuum dislocation density b(x) = >, b;d?(x — x;), 
j Lf 2, bx’): (x — x’) 
=— eo 11.101 
0(x) = fa x eae (VII.101) 
Alternatively, in Fourier space, 
‘ qx O(x") (x — x") _  b(a)-4 
= 2 2 iq:& 22, 


The angular fluctuations are thus related to correlations in dislocation density via 


(18(a)|?) = “Se (oe(a)o"(a)) (VII.103) 


where 


(b*(q)b?(q)) = pex etx (b*(x)b"(0)). (VII.104) 
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After the dislocations are unbound for T > T,, they interact via a screened Coulomb 
interaction, and 


lim (b™(x)b9(0)) x dage PS, = lim, (b*(q)b°(q)) « dap€?. (VII.105) 


|x|—-00 


Substituting into eq.(VII.103) leads to 


2 
, PON oo G 
lim, (18(a)| ) a 55. (VII.106) 
(In the low temperature phase, the neutrality of the charges at a large scale |q|~! requires 
the vanishing of b(q) as q — 0, and limg-.o (b“(q)b°(q)) Xx dags, leading to a finite 
(|8(a)|?).) 

Eq.(VII.106) implies that after the unbinding of dislocations the orientational fluctu- 


ations are still correlated. In fact such correlations would result from a Hamiltonian 


K 
BH = > / d’x(V0)?, with «= Ky x €?. (VII.107) 
The angular stiffness Kk, is known as the Frank constant. The bond angle order correlations 


now decay as 


(W(x) U*(0)) =e F(()-400)}") = ( f > | (VIL.108) 


Ix] 
with ny = 18/(7K 4). The quasi-long-range decay of orientational fluctuations leads to the 
appearance of a six fold intensity modulation in the diffraction pattern. The dislocations 
are thus not effective in completely destroying order, and their unbinding leads to the 
appearance of an orientationally ordered phase known as a hezxatic. The stiffness of the 


hexatic phase diverges at the transition to a solid according to eq.(VII.107). 


(c) Orientational order disappears at a higher temperature due to the unbinding of a new 
set of topological defects known as disclinations. These are very similar to the vortices in 


the XY model, except that since the bond angle is defined only up to 27/6, they satisfy 
, 2 
ve -ds = ra (VII.109) 


The energy cost of a disclination grows with system size as €, = 1K 4 |In(L/a)/36. Consid- 
ering the entropy of 21In(L/a), we find a disclination unbinding transition for K.4 < 72/7. 


This transition is in the universality class of the scalar Coulomb gas. The resulting high 
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temperature phase has neither orientational or translational order and is a conventional 
liquid. 

This scenario predicts that the melting of a two dimensional solid proceeds through 
an intermediate hexatic phase. In fact, computer simulations suggest that simple two 
dimensional solids, e.g. particles interacting via Lennard-Jones potentials, undergo a direct 
first order melting transition as in three dimensions. More complicated molecular systems 
in three dimensions, e.g. long polymers, are known to have intermediate liquid crystal 
phases. Liquid crystals have order intermediate between solid and fluid and are three 
dimensional analogs of the hexatic phase. Thin films made up of a few mono-layers of 


liquid crystals are the best candidates for examining the KTHNY melting scenario. 
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VIII. Dissipative Dynamics 


VITI.A Brownian Motion of a Particle 


Observations under a microscope indicate that a dust particle in a liquid drop under- 
goes a random jittery motion. This is because of the random impacts of the much smaller 
fluid particles. The theory of such (Brownian) motion was developed by Einstein in 1905 
and starts with the equation of motion for the particle. The displacement Z(t), of a particle 
of mass m is governed by, 

mi = sr; = a +e F inials). (VIIL.1) 
The three forces acting on the particle are: 

(i) A friction force due to the viscosity of the fluid. For a spherical particle of radius R, 
the mobility in the low Reynolds number limit is given by = (677R)~+, where 7 is 
the specific viscosity. 

(ii) The force due to the external potential V(#), e.g. gravity. 

(iii) A random force of zero mean due to the impacts of fluid particles. 

The viscous term usually dominates the inertial one (i.e. the motion is overdamped), 

and we shall henceforth ignore the acceleration term. Eq.(VIII.1) now reduces to a 


Langevin equation, 


# = 0(#) + q(t), (VIII.2) 
where u(#) = —yOV/0zZ is the deterministic velocity. The stochastic velocity, 7(t) = 
[ Faden), has zero mean, 

(A(t)) = 0. (VIIL3) 


It is usually assumed that the probability distribution for the noise in velocity is Gaussian, 


i.e. 


P |ni(t)| x exp - i dr a (VIII.4) 


Note that different components of the noise, and at different times, are independent, and 
the covariance is 


(na(t)ng(t’)) = 2D6a,86(t — t’). (VIII.5) 


141 


The parameter D is related to diffusion of particles in the fluid. In the absence of any 


potential, V(z) = 0, the position of a particle at time t is given by 


X(t) = £(0) +f drn(r). 


Clearly the separation £#(t) — Z(0) which is the sum of random Gaussian variables is itself 


Gaussian distributed with mean zero, and a varaince 


(#0) = z(0))") = / dry dry (71) -A(72)) = 3 x 2Dt. 


For an ensemble of particles released at Z(t) = 0, ie. with P(#,t=0) = 6°(#), the 


particles at time t are distributed according to 


1 oie e 
r= (Teme) ee |-apil) 


which is the solution to the diffusion equation 


OP _ pvp, 
at 


A simple example is provided by a particle connected to a Hookian spring, with 


V(#) = Ka?/2. The deterministic velocity is now v(#) = —ywK2#, and the Langevin 
equation, a = —wKZ + H(t), can be rearranged as 

d wKt > ukt— 

or [era (t)| = eh itt): (VIII.6) 


Integrating the equation from 0 to t¢ yields 
t 
etKt a(t) — #(0) = i dret*T Fr), (VIII.7) 
0 


and 


t 
#(t) = #(0)e "HE + | dre PK (- F(7), (VIIL.8) 
0 


Averaging over the noise indicates that the mean position, 


(E(t) = FO)e*™, (VIIL.9) 


decays with a characteristic relaxation time, T = 1/(4K). Fluctuations around the mean 


behave as 


2D6(71—T2) x3 
t 


((#t8) — @))?) =f drdrme HK Om Ger) Hea) 


t 
os | rae (VIII.10) 
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However, once the dust particle reaches equilibrium with the fluid at a temperature T, its 


probability distribution must satisfy the normalized Boltzmann weight 


, mer Kx? 
Peq.(Z) = (+) exp -- (VIII.11) 


yielding a) = 3kpT/K. Since the dynamics is expected to bring the particle to equilib- 
rium with the fluid at temperature T, eq.(VIII.10) implies the condition 


D=kpTp. (VIII.12) 


This is the Einstein relation connecting the fluctuations of noise to the dissipation in the 
medium. 

Clearly the Langevin equation at long times reproduces the correct mean and variance 
for a particle in equilibrium at a temperature T in the potential V(#) = Ka?/2, provided 
that eq.(VIII.12) is satisfied. Can we show that the whole probability distribution evolves 
to the Boltzmann weight for any potential? Let P(Z,t) = (2#|P(t)|0) denote the probability 
density of finding the particle at 7 at time t, given that it was at 0 at t = 0. This probability 
can be constructed recursively by noting that a particle found at x at time t+ ¢€ must have 


arrived from some other point 7’ at t. Adding up all such probabilities yields 
P(z@,t+e = paez'ra' (#|T.|2"), (VIII.13) 
where (Z|T.|Z") = (Z|P(e)|Z") is the transition probability. For « < 1, 


B=2'+0(# eth , (VIII.14) 
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where 7. = ji drn(r). Clearly, (ij-) = 0, and (n?) = 2De x 3, and following eq.(VIII.4), 


1 3/2 "2 
ay = Se AI. IL 
P(e) (<a) exp | a cv EEE IS) 


The transition rate is simply the probability of finding a noise of the right magnitude 
according to eq.(VII.14), and 


ae g— #! — (#")) 
(@|T(\#") = p(n) = (a) eed 


(i _ az) 2 (VIII.16) 


1 \3/2 
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By subdividing the time interval t, into infinitesimal segments of size €, repeated 


application of the above evolution operator yields 


P(#,t) = (z P(e)" 0) 


: 2 

@) Da(r) 1 (#-a@) (VIIL.17) 

-/ exp -{ dt -——_———— |. 
(0,0) N 0 4D 


The integral is over all paths connecting the initial and final points; each path’s weight 
is related to its deviation from the classical trajectory, z= v(£). The recursion relation 
(eq.(VIII.13)), 


2 PY ae cl (€-#'-ev(#’)) | __, 
P(Z,t) = fe x (=) exp gee P(#’,t—€), (VIII.18) 
can be simplified by the change of variables, 
Y=F'+ 0(z2')-F = 
(VIII.19) 


df =d? 2! (1+€V -0(2’)) = Bz! (1+ eV -o(#) + O(e)). 


Keeping only terms at order of €, we obtain 


‘i 3/2 
P(#,t) =[1—eV - o(@)] [ea ) ei P(#+ G9 —cd(Z),t—©) 


1 3/2 2 
se = aa 38> ee ae 
=[1-eV- v(z)| fe y (=) € x 


iYy — 2eyiv; + Priv; 0 
Pa. t)+ (g—eo(z))- VP + 24 os ee gp = + O(c?) 
=(1-—eV-v(2Z)] [P-eo-V + pvp + O(e€)!.. 


(VIII.20) 
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Equating terms at order of € leads to the Fokker-Planck equation, 


oP +V-J=0, with J=¢P-—DVP . (VIII.21) 


The Fokker-Planck equation is simply the statement of conservation of probability. The 
probability current has a deterministic component vP, and a stochastic part —DVP. A 
stationary distribution, OP/Ot = 0, is obtained if the net current vanishes. It is now 
easy to check that the Boltzmann weight, Peq.(%) « exp[—V(Z)/kpBT], with VPeq. = 
U Peq./(ukeT), leads to a stationary state as long as the fluctuation—dissipation condition 
in eq.(VIII.12) is satisfied. 
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VITI.B Equilibrium Dynamics of a Field 


The next step is to generalize the Langevin formalism to a collection of degrees of 
freedom, most conveniently described by a continuous field. Let us consider the order 
parameter field m(x,t) of a magnet. In equilibrium, the probability to find a coarse- 
grained configuration of the magnetization field is governed by the Boltzmann weight of 


the Landau-Ginzburg Hamiltonian 
~ da | 8 a 2 
Fla = fe x [En +um* + oy (Vm)" +--+]. (VIII.22) 


(To avoid confusion with time, the coefficient of the quadratic term is changed from t to 
r.) Clearly the above energy functional contains no kinetic terms, and should be regarded 
as the analog of the potential energy V(Z) employed in the previous section. To construct 
a Langevin equation governing the dynamics of the field ™(x), we first calculate the analo- 
gous force on each field element from the variations of this potential energy. The functional 
derivative of eq.(VIII.22) yields 

oH|m|] 


Poo)= anes = —rm; — 4um,|m|? + KV?m. (VIII.23) 


The straightforward analog of eq.(VIII.2) is 


Om,(x, t 
omits) = wFi(x) + ni(x,t), (VIII.24) 
with a random velocity, 77, such that 


(n(x, t)) = 0, and (ni(x, t)n;(x’, t’)) = 2D6;;6(x — x’)d(t — t’). (VIII.25) 


The resulting Langevin equation, 


Om(x, t 
ones.) = prim — 4uum mh + pKV7 mn + F(x, t), (VIII.26) 


is known as the time dependent Landau-Ginzburg equation. Because of the nonlinear term 
m’m, it is not possible to integrate this equation exactly. To gain some insight into 
its behavior we start with the disordered phase of the model which is well described by 
the Gaussian weight with u = 0. The resulting linear equation is then easily solved by 


examining the Fourier components, 
m(q,t) = : d¢x e'V*77i(x, t), (VIIL27) 
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which evolve according to 


om(q, t 3 a 
orn) = —p(r + Kq’) m(q,t) + 7(q, t). (VIII.28) 
The Fourier transformed noise, 
n(q,t) = [atx e'T* F(x, t), (VIII.29) 


has zero mean, (7;(q,t)) = 0, and correlations 


2D6;;67(x—x’) 5(t-t’) 
! Pee ee | —————————— 
(ni(a,t)nj(a’,t')) = | xd! eA (ni(x, t)y(x’, t')) 
; (VIII.30) 
=2D6,;6(t — t’) jax eared) 
=2D6;;6(t — t’)(27)764(q + q’). 
Each Fourier mode in eq.(VIII.28) now behaves as an independent particle connected 


to a spring as in eq.(VIII.6). Introducing a decay rate 


= p(r + Kq’), (VII.31) 


the evolution of each mode is similar to eq.(VIII.8), and follows 


t 
m(q,t) = m(q, 0)e7 1H" +f dr e~VDE-7) Fig, 2). (VIII.32) 
0 


Fluctuations in each mode decay with a different relaxation time T(q); (m(q,t)) = 
m(q, 0) exp|—t/7(q)]. When in equilibrium, the order parameter in the Gaussain model 
is correlated over the length scale € = af Kr. In considering relaxation to equilibrium, 
we find that at length scales larger than xi (or q < 1/€), the relaxation time saturates 
Tmax = 1/(pr). On approaching the singular point of the Gaussian model at r = 0, the 
time required to reach equilibrium diverges. This phenomena is know as critical slowing 
down, and is also present for the non-linear equation, albeit with modified exponents. The 
critical point is thus characterized by diverging length and time scales. For the critical 
fluctuations at distances shorter than the correlation length €, the characteristic time scale 
grows with wavelength as t(q) ~ (ukq?)~'. The scaling relation between the critical 
length and time scales is described by a dynamic exponent z, as T «x A*. The value of 


z = 2 for the critical Gaussian model is reminiscent of diffusion processes. 
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Time dependent correlation functions are obtained from 
2.D64;6(71—-T2) (277) 754 (q4+q’) 


t ; er———_. 
(mi(a,tmj(a).= f drydrge "Mima ta) ni (gq, 71) (a, 72) 
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t 
=(2n)457(q + q/) 2D6i; | dre 21a) (t-7) 
0 


/ D —27(q)t 
=(27)454(q + q')5;; te) (1 — e721 i) 


D 
too dsd / 
27)°d°(q + q )d;; ————— >. 
—— ( ) (q q’) J u(r + Kq?) 
(VIII.33) 
However, direct diagonalization of the Hamiltonian in eq.(VIII.22) with u = 0 gives 
K 
a (4K) a q)l?, (VIII.34) 


leading to the equilibrium correlation functions 


(ny(a),(q1)) = 2n)"8%(a + 4!)0,y EE. 


(VIII.35) 
Comparing equations (VIII.33) and (VIII.35) indicates that the long-time dynamics repro- 
duce the correct equilibrium behavior if the fluctuation—dissipation condition, D = kpTy, 
is satisfied. 

In fact, quite generally, the single particle Fokker-Planck equation (VIII.21) can be 


generalized to describe the evolution of the whole probability functional, P([7(x)],t), as 


OP (|7(x)], t) / dy ) OH oP 
ee —p——_— P — D—__ } . VIII.36 
Ot dm; (x) YSmilx) “a dm (x) ( 
For the equilibrium Boltzmann weight 
ACC ae (VIIIL.37) 
kpT 


the functional derivative results in 


bP eq. 1 OH 
= ———_ —— P,,.. VIII.38 
dm; (x) kpT 6m,(x) Pea. ( ) 


The total probability current, 


oH D oH 


J{h(x)| =  Sma(x) + kpT 5mi(x) eq.) (VIII.39) 


vanishes if the fluctuation—dissipation condition, D = wkgT, is satisfied. Once again, the 


Einstein equation ensures that the equilibrium weight indeed describes a steady state. 


148 


VIII.C Dynamics of a Conserved field 


In fact it is possible to obtain the correct equilibrium weight with q dependent mobility 


and noise, as long as the generalized fluctuation—dissipation condition, 
D(a) = kT ula), (VIIT.40) 


holds. This generalized condition is useful in considering the dissipative dynamics of a 
conserved field. The prescription that leads to the Langevin equations (VIII.23)—(VII.25), 
does not conserve the field in the sense that [ d’xm(x, t) can change with time. (Although 
this quantity is on average zero for r > 0, it undergoes stochastic fluctuations.) If we are 
dealing with the a binary mixture (n = 1), the order parameter which measures the 
difference between densities of the two components is conserved. Any concentration that 
is removed from some part of the system must go to a neighboring region in any realistic 


dynamics. Let us then consider a dynamical process constrained such that 


d 5 = Om(x,t) = 
7 dx (x, t) = i dx ar 0. (VIII.41) 


How can we construct a dynamical equation that satisfies eq.(VIII.41)? The integral clearly 


vanishes if the integrand is a total divergence, i.e. 


Om (x, t : 
emit) = —V-j; + n(x, t). (VIII.42) 
The noise itself must be a total divergence, 7; = —V -o;, and hence in Fourier space, 


(ni(a,t)) =0, and (n(q,t)nj(a’,t’)) = 2D6,;q76(t — t')(27)464(q+q’).  (VIII.43) 


We can now take advantage of the generalized Einstein relation in eq.(VIII.40) to ensure 


the correct equilibrium distribution by setting, 


Ji = pV - (- oo] : (VIII.44) 


The standard terminology for such dynamical equations is provided by Hohenberg and 
Halperin: In model A dynamics the field ™ is not conserved, and the mobility and 
diffusion coefficients are constants. In model B dynamics the field ™ is conserved, and 
ji = —pV? and D = —DV?. 
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Let us now reconsider the Gaussian model (u = 0), this time with a conserved order 


parameter, with model B dynamics 


Om(x, t 
ernie) = prV?it — pK Vm + A(x, t). (VIIL.45) 
The evolution of each Fourier mode is given by 


Om(q, t) 
ot 


m(q, t) 
T(q) 


Because of the constraints imposed by the conservation law, the relaxation of the field 


= —pq?(r + Kq?)m(q, t) + 7(q,t) = - + H(q, t). (VIII.46) 


is more difficult, and slower. The relaxation times diverge even away from criticality. 
Depending on wavelength, we find scaling between length and time scales with dynamic 


exponents z, according to 


T(q) = 


1 q- torqgce (z= 2) 
> & F TII.4 
ugq?(r + Kq?) { Ge sora Re" (z =4) wee 


The equilibrium behavior is unchanged, and 


De nD 
lim (|7(q, t)|?) = n——————- = ——.., VIII.48 
ule, ug?(r+ KG?) wlr + Kq?) 
as before. Thus the same static behavior can be achieved by different dynamics. The 
static exponents (e.g. v) are determined by the equilibrium (stationary) state and are 
unchanged, while the dynamic exponents may be different. As a result, dynamical critical 


phenomena involve many more universality classes than the corresponding static ones. 
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VIII.D Generic scale invariance in equilibrium systems 


We live in a world full of complex spatial patterns and structures such as coastlines 
and river networks. There are similarly diverse temporal processes generically exhibiting 
“1 /f”—noise, as in resistance fluctuations, sand flowing through an hour glass, and even in 
traffic and stock market movements. These phenomena lack natural length and time scales 
and exhibit scale invariance and self-similarity. The spacial aspects of scale invariant sys- 
tems can be characterized using fractal geometry|1]. In this section we explore dynamical 
processes that can naturally result in such scale invariant patterns. 

Let us assume that the system of interest is described by a scalar field m(x), distributed 
with a probability P|m]. Scale invariance can be probed by examining the correlation func- 
tions of m(x), such as the two point correlator, C(|x—y|) = (m(x)m(y)) — (m(x)) (m(y)). 
(It is assumed that the system has rotational and translational symmetry.) In a system 
with a characteristic length scale, correlations decay to zero for separations z = |x—y| > €. 
By contrast, if the system possesses scale invariance, correlations are homogeneous at long 
distances, and lim,_,., C(z) ~ 2%. 

As we have seen, in equilibrium statistical mechanics the probability is given by Peg « 
exp (—3H{m]) with 8 = (kgT)~+. Clearly at infinite temperature there are no correlations 
for a finite Hamiltonian. As long as the interactions in H[{m] are short ranged, it can 
be shown by high temperature expansions that correlations at small but finite @ decay 
as C(z) « exp(—z/€), indicating a characteristic length-scalet. The correlation length 
usually increases upon reducing temperature, and may diverge if the system undergoes a 
continuous (critical) phase transition. At a critical transition the system is scale invariant 
and C(z) x z?~4—-", However, such scale invariance is non-generic in the sense that it can 
be obtained only by precise tuning of the system to the critical temperature. Most scale 
invariant processes in nature do not require such precise tuning, and therefore the analogy 
to the critical point is not particularly instructive[2][3]. 

We shall frame our discussion of scale invariance by considering the dynamics of a 
surface, described by its height h(x,t). Specific examples are the distortions of a soap 
film or the fluctuations on the surface of water in a container. In both cases the minimum 


energy configuration is a flat surface (ignoring the small effects of gravity on the soap film). 


* It is of course possible to generate long-range correlations with long ranged interac- 
tions. However, it is most interesting to find out how long-range correlations are generated 


from local, short ranged interactions. 
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The energy cost of small fluctuations for a soap film comes from the increased area and 


surface tension 0. Expanding the area in powers of the slope results in 


He= of ax yi + (Vh)? — 1 pe 5 f ax (Vh)? . (VIIL.49) 


For the surface of water there is an additional gravitational potential energy, obtained by 


adding the contributions from all columns of water as 


h(x) 
H, = fax pgh(x) = O8 fatten? (VIII.50) 
0 


The total (potential) energy of small fluctuations is thus given by 
H= jax E (Vh)? + cen (VIIL.51) 


with the second term absent for the soap film. 


The corresponding Langevin equation, 


Oh(x, t) 
Ot 


= —ppgh + poV7h + n(x, t), (VIII.52) 


is linear, and can be solved by Fourier transforms. Starting with a flat interface, h(x,t = 


0) = h(q,t = 0) = 0, the profile at time t is 


d‘q —iq:x ‘ - +oq7)(t—T 
e.cea / (an)a° i / dre Hester Mt-7) n(q, t), (VIII.53) 
The average height of the surface, H = f d?x (h(x,t)) /L@ is zero, while its overall width 
is defined by 
1 is d‘q 
2 = d DK. 2 
where L is the linear size of the surface. Similar to Eq.(VIII.33), we find that the width 


grOWS as 


weer [2 2 G_e-2nar 
n= | oe eA (1 ). (VIIL55) 


There are a range of time scales in the problem, related to characteristic length scales as 
in Eq.(VIII.31). The shortest time scale, tmin « a?/(wo), is set by an atomic size a. The 
longest time scale is set by either the capillary length (A. = \/a/pg) or the system size (L). 
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For simplicity we shall focus on the soap film where the effects of gravity are negligible and 


tmax « L?/(c). We can now identify three different ranges of behavior in Eq.(VIII.55): 


(a) 


(b) 


For t < tmin, none of the modes has relaxed since 7(q)t < 1 for all gq. Each mode 


grows diffusively, and 


ee ee 
wit,L) = ‘i (On)4 >(q) 2y(q)t = ae (VIII.56) 


For t >> tmax, all modes have relaxed to their equilibrium values since 7(q)t >> 1 for 


all q. The height fluctuations now saturate to a maximum value given by 


d‘q D 
2(t-G -([S= 11.57 
w(t) = | aoa ae (VIII.57) 


The saturated value depends on the dimensionality of the surface, and in a general 


dimension d behaves as 


D ane ford De “(se Se ) 
w GL) 6 —<In(L/a) * ford=2, -(=0T)., (VIII.58) 
aida We as ford<2, (x = 75%) 


where we have defined a roughness exponent x that governs the divergence of the 
width with system size via limp... w(t, L) « LX. (The symbol 0* is used to indicate 
a logarithmic divergence.) The exponent of y = 1/2 in d = 1 indicates that the one 
dimensional interface fluctuates like a random walk. 

For tmin < t < tmax only a fraction of the shorter length scale modes are saturated. 
The integrand in Eq.(VIII.55) (for g = 0) is made dimensionless by setting y = poq?t, 


and 


D 
w(t, L) « — je Go (1 — eure) 
[Lo 


D 1 i t/tmin d—4 
26)" rate eem 
t 


(VIII.59) 


The final integral is convergent for d < 2, and dominated by its upper limit for d > 2. 
The initial growth of the width is described by another exponent (, defined through 


limy.0 w(t, L) « t?, and 


eat ford 2,0 (8 = 0) 
w(t, Lyx 2iIn(t/tmn) ford=2, (6=07) (VIII.60) 
Canara tO? ford<2, (@=(2-—d)/4) 
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The exponents yx and @ also describe the height—height correlation functions which 
assumes the dynamic scaling form 


([atx, t) — h(x’, ty) = |x —x'|/?X gq (=) (VIII.61) 


|x — x’ |? 


Since equilibrium equal time correlations only depend on |x — x’|, lim,—o g(y) should be 
a constant. On the other hand correlations at the same point can only depend on time, 
requiring that limy_... g(y) « y?X/*, and leading to the exponent identity 3 = y/z. 

This scale invariance is broken when the gravitational potential energy, is added to 
the Hamiltonian. The correlations now decay as C(z) « exp(—z/A-) for distances larger 
than the capillary length. What is the underlying difference between these two cases? The 
presence of gravity breaks the translational symmetry, H[h(x) + c] = H[h(x)]. It is this 
continuous symmetry that forbids the presence of a term proportional to [ d’xh(x)? in the 
Hamiltonian and removes the corresponding length scale. (The coefficient of the quadratic 
term is usually referred to as a mass in field theoretical language.) The presence of a 
continuous symmetry is quite a general condition for obtaining generic scale invariance 
(GSI)[4]._ As discussed in previous chapters, there are many low temperature phases of 
matter in which a continuous symmetry is spontaneously broken. The energy cost of small 
fluctuations around such a state must obey the global symmetry. The resulting excitations 
are the “massless” Goldstone modes. We already discussed such modes in connection 
with magnons in ferromagnets (with broken rotational symmetry), and phonons in solids 
(broken translational symmetry). All these cases exhibit scale invariant fluctuations. 

In the realm of dynamics we can ask the more general question of whether temporal 
correlations, e.g. C(|x — x’|,t — t’) = (h(x, t)h(x’, t’))., exhibit a characteristic time scale 
T, or are homogeneous in t — t’. It is natural to expect that scale invariance in the 
spacial and temporal domains are closely interlinked. Establishing correlations at large 
distances requires long times as long as the system follows local dynamical rules (typically 
(t — t') « |x —x’|*). Spacial scale invariance thus implies the lack of a time scale. The 
converse is not true as dynamics provides an additional possibility of removing time scales 
through a conservation law. We already encountered this situation in examining the model 
B dynamics of the surface Hamiltonian in the presence of gravity. Equation (VIII.47) 
indicates that, even though the long wavelength modes are massive, the relaxation time of 


a mode of wavenumber q diverges as 1/q? in the q — 0 limit. 
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Symmetries and conservation laws are intimately linked in equilibrium systems. Con- 
sider a local Hamiltonian that is invariant under the symmetry H[h(x) + c] = H[h(x)]. 
Since H can only depend on VA and higher derivatives, 


me a a (VIII.62) 


= pL = >ft 


Even for model A dynamics, the deterministic part of the velocity is the divergence of a 
current and conserves { d?xh(x,t). The conservation is only statistical and locally broken 
by the non-conserved noise in model A. The above result is a consequence of Noether’s 


theorem. 


VIIIL.E Non-equilibrium dynamics of open systems 


We have to be cautious in applying the methods and lessons of near equilibrium dy- 
namics to the various processes in nature which exhibit generic scale invariance. Many 
such systems, such as a flowing river or a drifting cloud, are very far from equilibrium. 
Furthermore, they are open and extended systems constantly exchanging particles and 
constituents with their environment. It is not clear that there is any Hamiltonian that 
governs the dynamics of such processes and hence the traditional approach presented ear- 
lier is not necessarily appropriate. However, the robust self-similar correlations observed 
in these systems|1] suggests that they can be described by stationary scale invariant prob- 
ability distributions. This section outlines a general approach to the dynamics of open and 
extended systems that is similar in spirit to the construction of effective coarse-grained 
field theories described in Chapter 2. Let us again consider the dynamics of a static field, 
h(x, t): 

1. The starting point in equilibrium statistical mechanics is the Hamiltonian 1{|h]|. Lan- 
dau’s prescription is to include in # all terms consistent with the symmetries of the 
problem. The underlying philosophy is that in a generic situation an allowed term is 
present, and can only vanish by accident. In the case of non-equilibrium dynamics we 
shall assume that the equation of motion is the fundamental object of interest. Over 
sufficiently long time scales, inertial terms (x 07h) are irrelevant in the presence of 


dissipative dynamics, and the evolution of h is governed by 


deterministic stochastic 
_—~__ —— 
Oph(x, t) = v[h(x,t)] + (x,t). (VIII.63) 
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2. If the interactions are short ranged, the velocity at (x,t) depends only on h(x, t) and 


a few derivatives evaluated at (x,t), i.e. 
v(x, t) = v(h(x, t), VA(x,t),---). (VIII.64) 


3. We must next specify the functional form of deterministic velocity, and the correlations 
in noise. Generalizing Landau’s prescription, we assume that all terms consistent with 
the underlying symmetries and conservation laws will generically appear in v. The 
noise, 7(x,t), may be conservative or non-conservative depending on whether there 
are only internal rearrangements, or external inputs and outputs. 

Corollary: Note that with these set of rules there is no reason for the velocity to be deriv- 
able from a potential (v 4 —fdH/dh), and there is no fluctuation-dissipation condition 
(D F jikgT). It is even possible for the deterministic velocity to be conservative, while 
the noise is not. Thus various familiar results of near equilibrium dynamics may no longer 
hold. 

As an example consider the flow of water along a river (or traffic along a highway). 
The deterministic part of the dynamics is conservative (the amount of water, or the number 
of cars is unchanged). Hence the velocity is the divergence of a current, v = “Vi [h]. The 
current, 7, is a vector, and must be constructed out of the other two vectorial quantities 
in the problem: the gradient operator V, and the average transport direction ¢. (The unit 
vector ¢ points along the direction of current or traffic flow.) The lowest order terms in 


the expansion of the current give 
Lia. tk Xr Sa 
—j =t(ah — a ++.) 4yVhAt+ wot(f-V)ht--:. (VIII.65) 


The components of current parallel and perpendicular to the net flow are 


Sagi — Ap? sali 
{ fl a ae + (v1 + ¥2) Oh + (VIII.66) 
Fit pO |: Bi 
The resulting equation of motion is 
Oh(x, t x 
ones) = O\(ah — 3”) + (m+ v2)Oth +07h+--++n(x,t). (VIII.67) 


In the absence of external inputs and outputs (no rain, drainage, or exits), the noise is also 


conservative, with correlations, 


(n(q,t)) =0, and (n(q,t)n(a’,t’)) = 2(Dyaqj + Digi) 6(t — t') (27) 0%(q +’). 
(VIII.68) 
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Note that the symmetries of the problem allow different noise correlations parallel and 
perpendicular to the net flow. 
Equations (VIII.67) and (VIII.68) define a driven diffusion system (DDS)|5],/6]. The 


first term in Eq.(VIII.67) can be eliminated by looking at fluctuations in a moving frame, 
h(ay,x1,t) > A(x —at,x1,4), (VIII.69) 


and shall be ignored henceforth. Neglecting the non-linear terms, these fluctuations satisfy 


the anisotropic noisy diffusion equation 


Oh(x, t) 
Ot 


= yjOpht+viOjh+n(x,t), (VIII.70) 


where 1) = 41 + v2 and vy, =. The Fourier modes now relax with characteristic times, 


1 
7(q) = —,>——>. (VIIL.71) 
i Mig + YL 


Following the steps leading to Eq.(VIII.48), 


Diq + Dia 
lim ([A(q,2)?) = —_ >. (VIIL72) 
=O Yq TUL) 
The stationary correlation functions!’ in real space now behave as(4][6] 
Di qi + Dy, qe 


d—1 
((n(0 = 0(0))") = | SAL (2 — 2eos(qry an x2) 


(> —1) 
x | —— —- — 
Vi Y|| 
(VIII.73) 


Note that these correlations are spatially extended and scale invariant. This is again 


WyGG + VLG 


a consequence of the conservation law. Only in the special case where D,/v_ = Dy/vj 
is the Einstein relation (D(q) « v(q)) satisfied, and the fluctuations become uncorrelated 
(C(x) « 64(x)D/v). The results then correspond to model B dynamics with a Hamiltonian 
H x fd¢xh?. This example illustrates the special nature of near equilibrium dynamics. 


The fluctuation—dissipation condition is needed to ensure approach to the equilibrium state 


* In non-equilibrium circumstances we shall use the term stationary to refer to behavior 


at long times. 
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where there is typically no scale invariance. On the other hand, removing this restriction 
leads to GSI in a conservative system. (D1 /v. #4 D)/™ is like having two different 
temperatures parallel and perpendicular to the flow.) 

Let us now break the conservation law stochastically by adding random inputs and 
outputs to the problem (in the form of rain, drainage, or exits). The properties of the 


noise are now modified to 
(n(x,t)) =0, and (n(x, t)n(x’, t’)) = 2Dd(a) — a) 09" (x1, — x’, )d(t—t’). (VIII.74) 


In the stationary state, 


, D 
Jim, (|h(q, t)|?) = nae (VIII.75) 
in Fourier space, and 
2a 
((h(x) = h(0))”) « D (via? 4 4%) ae (VIII.76) 


in real space. Except for the anisotropy, this is the same result as in Eq.(VIII.61), with 
x= (2-d)/2. 

How are the results modified by the nonlinear term (—A0\h?/2) in Eq.(VIII.67)? We 
first perform a simple dimensional analysis by rescaling lengths and time. Allowing for 
anisotropic scalings, we set 2) — bx), accompanied by t > b*t, 71 — b6#,, and h — bXh. 
Eq.(VIII.67) is now modified to 


_,Oh = as 
era ue Oth + v4 A — 
where Eq.(VIII.74) has been used to determine the scaling of 7. We thus identify the bare 


scalings for these parameters as 


a 
D 


bx ah? +b 2/2—(d—-1)¢/2 V/2y (VIII.77) 


yy > b°*y| 
Vl, 7 be op) 


(VIII.78) 
ee aa. 


= b2—2x-S(d-1)-1 pp 


In the absence of A, the parameters can be made scale invariant (i.e. independent of b), 
by the choice of G9 = 1, 2) = 2, and yo = (2 — d)/2, as encountered before. However, with 


this choice, a small ’ will grow under rescaling as 


A= bY, with yo= ams (VIII.79) 
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Since the non-linearity grows larger under scaling it cannot be ignored in dimensions d < 4. 

Equations (VIII.79) constitute a simple renormalization group (RG) that is valid 
close to the fixed point (scale invariant equation) corresponding to a linearized limit. To 
calculate the RG equations at finite values of nonlinearity in general requires a perturbative 
calculation. Sometimes there are exact non-renormalization conditions that simplify the 
calculation and lead to exponent identities. Fortunately there are enough such identities 
for Eq.(VIII.78) that the three exponents can be obtained exactly. 

1. As the nonlinearity is proportional to qj is Fourier space, it does not generate under 
RG any contributions that can modify v,. The corresponding ‘bare’ scaling of v 
in Eqs.(VII.78) is thus always valid; its fixed point leads to the exponent identity 
g26 = 0: 

2. As the nonlinearity is in the conservative part, it does not renormalize the strength 
of the non-conservative noise. The non—renormalization of D leads to the exponent 
identity z — 2y — ¢(d-—1) —1 = 0. (This condition has a natural counterpart in 
equilibrium model B dynamics, z — 2. —d—2 = 0, leading to the well known relation, 
z=4-7.) 

3. Eq.(VII.67) is invariant under an infinitesimal reparameterization x) — a —dAt, t > 
t, if hh >h+0. Note that the parameter \ appears both as the coefficient of the non- 
linearity in Eq.(VIII.67) and as an invariant factor relating the x) and h reparameteri- 
zations. Hence any renormalization of the driven diffusion equation that preserves this 
symmetry must leave the coefficient \ unchanged, leading to the exponent identity 
i ae eee a8 
The remaining parameter, 1, does indeed follow a non-trivial evolution under RG. 

However, the above three exponent identities are sufficient to give the exact exponents in 


all dimensions d < 4 as[2] 


——— (VIII.80) 


VIII.F Dynamics of a growing surface 


The rapid growth of crystals by deposition, or molecular beam epitaxy, is an important 
technological process. It also provides the simplest example of a non-equilibrium evolution 
process|7]. We would like to understand the dynamic scaling of fluctuations inherent to 
this type of growth. To construct the dependence of the local, deterministic velocity, v, on 


the surface height, h(x,t), note that: 


159 


(1) As long as the rearrangements of particles on the surface can result in holes and 
vacancies, there is no conservation law. 

(2) There is a translation symmetry, v[h(x) + c] = v[h(x)], implying that v depends only 
on gradients of h(x). 

(3) For simplicity, we shall focus on isotropic surfaces, in which all directions in x are 
equivalent [8]. 

(4) There is no up-down symmetry, i.e. vu[h(x)] 4 —v[—h(x)]. The absence of such 
symmetry allows addition of terms of both parities. 

With these conditions, the lowest order terms in the equation of motion give|9], 


Oh(x, t) 


xT +uVPA+ A (Vh)? +---+ (x,t), (VIII.81) 


with the non-conservative noise satisfying the correlations in Eq.(VII.74). 

In Eq.(VHI.81), wu is related to the average growth velocity. In fact, the coefficients 
of all even terms must be proportional to u as they all vanish in the symmetric case with 
no preferred growth direction. The constant u is easily removed by transforming to a 
moving frame, h — h—ut, and will be ignored henceforth. The first non-trivial term is the 
nonlinear contribution, \(Vh)?/2. Geometrically this term can be justified by noting that 
growth by addition of particles proceeds through a parallel transport of the surface gradient 
in the normal direction. (See the inset to Fig. 9.1.) This term cannot be generated from 
the variations of any Hamiltonian, i.e. v 4 —y“dH|h|/dh. Thus, contrary to the equilibrium 
situation (Noether’s theorem), the translational symmetry does not imply a conservation 
law, v # —Vj. 

Further evidence of the relevance of Eq.(VIII.81) to growth phenomena is provided by 
examining deterministic growth. Consider a slow and uniform snowfall, on an initial profile 
which at t = 0 is described by ho(x). The nonlinear equation can in fact be linearized with 


the aid of a “Cole-Hopf” transformation, 
aN 
W (x,t) = exp Puce 0) : (VIII.82) 


The function W(x, t) evolves according to the diffusion equation with multiplicative noise, 


OW (x, t) 


=v’V?W+ = W(x, t). (VIIL83) 
Ot 2v 
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In the absence of noise, 7(x,t) = 0, Eq.(VII.83) can be solved subject to the initial 
condition W(x, t = 0) = exp|Aho(x)/2v], and leads to the growth profile, 
/ 
het) = ~ ln { [atx exp _ + nix 0) \ ; (VIII.84) 
It is instructive to examine the v — 0 limit, which is indeed appropriate to snow falls since 
there is not much rearrangement after deposition. In this limit, the integral in Eq.( VIII.84) 
can be performed by the saddle point method. For each x we have to identify a point x’ 
which maximizes the exponent, leading to a collection of paraboloids described by 
{ An P= xe \ 
h(x, t) = max 4 ho(x') - ———— >? . (VIII.85) 
2At ai 
Such parabolic sequences are quite common in many layer by layer growth processes in 
nature, from biological to geological formations. The patterns for 4 = 1 are identical 
to those obtained by a geometrical method of Huygens, familiar from optics. The growth 
profile (wave front) is constructed from the outer envelop of circles of radius At drawn from 
all points on the initial profile. The nonlinearity in Eq.(VUI.81) algebraically captures this 


process of expanding wave fronts. 


VIII.1. Deterministic growth according to eq,(VIII.81) leads to a pattern of coarsening 
paraboloids. In one dimension, the slope of the interface forms ‘shock fronts.’ Inset depicts 


projection of lateral growth on the vertical direction. 
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As growth proceeds, the surface smoothens by the coarsening of the parabolas. What 
is the typical size of these features at time t? In maximizing the exponent in Eq.(VIII.85), 
we have to balance a reduction |x — x’|?/2At, by a possible gain from ho(x’) in selecting 
a point away from x. The final scaling is controlled by the roughness of the initial profile. 


Let us assume that the original pattern is a self-affine fractal of roughness x, i.e. 
lho(x) — ho(x’)| ~ |x — x'|*. (VIII.86) 


(According to Mandelbrot, . * 0.7 for mountains[1].) Balancing the two terms in 


Eq.(VIII.85) gives 


(dx)? 


(62)* ~ — devt?, with z+x=2. (VIII.87) 


For example, if the initial profile is like a random walk in d= 1, y = 1/2, and z = 3/2. 
This leads to the spreading of information along the profile by a process that is faster than 
diffusion, dx ~ t?/3. 


Note that the slope, @(x,t) = —AVA(z,t), satisfies the equation, 


Dé : “ 
pete!) = a +0- Vi =vV?s —AVn. (VIIL88) 


The above is the Navier-Stokes equation for the velocity of a fluid of viscosity v, which 
is being randomly stirred by a conservative force[10], F = —AV7n. However, the fluid is 
vorticity free since 


Q=Vxtd=-AV x Vh=0. (VIII.89) 


This is the Burgers’ equation|11], which provides a simple example of the formation of 
shock waves in a fluid. The gradient of Eq.(VIII.85) in d = 1 gives a saw tooth pattern of 
shocks which coarsen in time. Further note that in d = 1 Eq.(VIII.88) is also equivalent 
to the driven diffusion equation of (VIII.67), with v playing the role of h. 

To study stochastic roughening in the presence of the nonlinear term, we carry out 
a scaling analysis as in Eq.(VIHI.77). Under the scaling x — bx, t — b*t, and h — bxh, 
Eq.(VIII.81) transforms to 


pore = pbX-2V27h + Spx (Vh)? + 7 (bx, b*t). (VIII.90) 
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The correlations of the transformed noise, 1/(x, t) = b*~*7(bx, b*t), satisfy 


(in! (x, t)n! (x’, t')) =b2*-2* 2D 64 (x — x')b~46(t — t')b-* 


(VIII.91) 
=b?-4-2x 2D 64 (x — x’)6(t — t’) 
Following this scaling the parameters of Eq.(VIII.81) are transformed to 
vy b? 4p 
Ni BS og (VIIL.92) 
D3 $24 p 


For A = 0, the equation is made scale invariant upon the choice of zo = 2, and yo = 
(2—d)/2. Close to this linear fixed point, \ scales to b7°+X0-?) = b2?-4%/?). and is a 
relevant operator for d < 2. In fact a perturbative dynamic renormalization group suggests 
that it is marginally relevant at d = 2, and that in all dimensions a sufficiently large A 
leads to new scaling behavior. (This will be discussed further in the next chapter.) 

Are there any non—renormalization conditions that can help in identifying the expo- 
nents of the full nonlinear stochastic equation? Note that since Eqs.(VIHI.81) and (VIII.88) 
are related by a simple transformation, they must have the same scaling properties. Since 
the Navier-Stokes equation is derivable from Newton’s laws of motion for fluid particles, 
it has the Galilean invariance of changing to a uniformly moving coordinate frame. This 
symmetry is preserved under renormalization to larger scales and requires that the ratio 
of the two terms on the left hand side of Eq.(VIII.88) (Ov and v- Vv) stays at unity. In 
terms of Eq.(VIII.81) this implies the non—-renormalization of the parameter A, and leads 
to the exponent identity 

Yo ea: (VIII.93) 


Unfortunately there is no other non—renormalization condition except in d = 1. Fol- 
lowing Eq.(VIII.36), we can write down a Fokker—Planck equation for the evolution of the 
configurational probability as, 

OP ([h(x)], t) / }6 9 aN 9 oP 
ot = | dk | [V72h + (Va)? |) P- DI]. VIII.94 
at rae k Cae Ocala 5h(x) ( ) 
Since Eq.(VHI.81) was not constructed from a Hamiltonian, in general we do not know 
the stationary solution at long times. In d = 1, we make a guess and try a solution of the 


form 


Pol[h(x)| « exp = / te(2nh) (VIII.95) 
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Since 


= 5 (02h) Po, (VIII.96) 


Eq.(VIII.94) leads to, 


OP 0 dP o 2 r 2 V 92 

ee eS See Ss =p 
FE [eae (vatn + 5 (Och) poh 

: (VIII.97) 

Vv 


_ V a2 Bits 


oh) 25 


We have thus identified the stationary state of the one dimensional equation. (This pro- 
cedure does not work in higher dimensions as it is impossible to write the final result 
as a total derivative.) Surprisingly, the stationary distribution is the same as the one in 
equilibrium at a temperature proportional to D/v. We can thus immediately identify the 
roughness exponent y = 1/2, which together with the exponent identity in Eq.(VIII.93) 
leads to z = 3/2, i.e. super-diffusive behavior. 

The values of the exponents in the strongly non-linear regime are not known exactly 
in higher dimensions. However, extensive numerical simulations of growth have provided 
fairly reliable estimates[7|. In the physically relevant case (d = 2) of a surface grown in 
three dimensions, x ~ 0.39 and z = 1.61/12]. A rather good (but not exact) fit to the 


exponents in a general dimension d is the following estimate by Kim and Kosterlitz/13], 


xe — and 2% ——. (VIII.98) 
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